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The mathematical theory of neutron chain reacting systems bears considerable similarity to 


Ocroser, 1952 


the mathematical theory of quantum mechanics since the Schroedinger equation can be viewed 
as the equation for an infinitely large chain reactor in which the potential function is replaced 
by the distribution of neutron absorption cross sections. The chain reactor theory, not being 
self-adjoint, is somewhat more complicated than the Schroedinger equation theory. Neverthe- 
less, perturbation theory, time dependent theory, and the theory of the stationary, i.e., just crit- 


INTRODUCTION 


HE purpose of this review is to present the 
scientific—or more accurately—the mathe- 
matical principles which underlie the design and 
operation of high-powered nuclear chain re- 
actors. The theory of these devices is a beautiful 
example of the general theory of linear operators; 
because of wartime restrictions the subject has, 
on the whole, not been very widely disseminated. 
Work on the theory of nuclear chain reactors— 
more perhaps than work in any other field of 
theoretical physics—has been almost entirely 
confined to the large national laboratories rather 
than to universities. This situation stems only 
partly from secrecy; nuclear reactor theory is 
now considered generally declassifiable by the 
Atomic Energy Commission. | think a more im- 
portant reason is that nuclear reactors, because 
of their expense, have yet to become part of the 
university laboratory, and so interest in them 
there has been comparatively slight. All of my 
report will be based on work which has gone on 
in the American atomic energy laboratories; 


ical, state can be stated in much the same way as is done in quantum mechanics. 
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much of it will be familiar to those who have 
been associated with one of these laboratories 
during the war. Contributors to the theory in- 
clude many well-known theoretical physicists. 
Because of the wartime circumstances under 
which much of the work took place, the identifi- 
cation of each individual’s contribution has be- 
come very difficult. Suffice to say that contribu- 
tions to the theory discussed here were made by 
E. P. Wigner, F. L. Friedman, E. Fermi, G. 
Placzeck, G. Volkoff, J. A. Wheeler, R. P. 
Feynman, H. Hurwitz, H. Brooks, L. Nordheim, 
H. Soodak, R. F. Christy, and G. Young. 

We shall consider a nuclear chain reaction to 
be characterized by the energy, spatial, and 
temporal distribution ®(£,r,¢) of the neutron 
flux, the neutron flux at energy E being defined 
as the number of neutrons per unit E per cm® 
multiplied by the neutron speed. The angular 
distribution of the flux, i.e., the number of neu- 
trons of particular speed going in a certain direc- 
tion, we do not consider; nor do we consider 
statistical fluctuations. On this account the 
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theory described here may be called ‘‘elemen- 
tary.”’ In effect, the assumption will be made 
that the angular distribution is almost spherically 
symmetric. The calculation of the exact angular 
distribution falls into the province of “‘trans- 
port” theory and is very much more complicated 
than what is treated here. 

It is the purpose of nuclear reactor theory to 
compute ®(£,r,?¢) for various dispositions of 
fissionable material, moderator, and absorbers, 
given the nuclear cross sections of the constitu- 
ents. From the flux distribution the power dis- 
tribution z(r, ¢) is 


a(r, N= fone, r)&(E, r, t)dE, 


where o;(E£,r) is the fission cross section! per 
cm® at energy E and position defined by the 
position vector r. Since the flux at energy E 
multiplied by the fission cross section per cm® 
is the number of fissions per cm’ per second 
caused by neutrons of energy E, the integral 
over energy gives the total number of fissions 
per cm* per second caused by all neutrons, and 
this is proportional to the power output. The 
power distribution is, of course, the major datum 
required for the engineering design of a high- 
powered nuclear reactor. 


THE THERMAL REACTOR; THE FUNDAMENTAL 
THEOREM 


The simplest and most basic situation is the 
steady state in a reactor in which all fissions and 
all absorptions occur only at thermal energy, 
and all cross sections are independent of position. 
An unreflected thermal neutron chain reactor, 
such as is almost exemplified by the huge graphite 
reactor at Oak Ridge, is an example of such a 
system.? The thermal neutron flux distribution, 
®,(r), at thermal energy E, obeys 


Se etctuns f &,(r’)P(E,, 1, r’)dr’=0, (1) 


Ve 


1The cross section per cm* is sometimes called the 
“macroscopic” cross section; it is the cross section per 
atom multiplied by the number of atoms per cm*. The 
mean free path of a neutron is the reciprocal of the macro- 
scopic cross section. 

? The lattice structure of unenriched reactors is ignored 
throughout this discussion. 


where D=diffusion coefficient, o,=thermal ab- 
sorption cross section per cm’, and y=number of 
neutrons per fission. The first term, DA®,, is the 
net increase in neutrons per cm® per second due 
to neutron diffusion; the second term, o,®,, is 
the number of neutrons absorbed per cm? per 
second ; the third term is the number of neutrons 
which become thermalized per cm’ at r as the 
result of fission occurring throughout the reactor. 
Each fission produces v fast neutrons, and these 
fast neutrons, on being slowed by elastic collision 
with the moderator, appear as slow neutrons at 
some distance away from their place of birth. 
The quantity P(E,,r,r’), called the slowing- 
down kernel, is the number of neutrons which 
reach energy E, in the reactor at r due to a fast 
neutron produced from fission at r’ and energy 
Eo. The integration is over the reactor volume 
V. since neutrons produced by fission at any 
point in the reactor have some chance of appear- 
ing as slow neutrons at r. 

The slowing-down kernel P(E,, 1,1’), in gen- 
eral, is not a displacement kernel, i.e., does not 
depend simply on the distance |r—r’| between 
the source and the field points, because if either 
r or r’ is close to a boundary, neutrons can leak 
out of the system. Thus a distinction must be 
made between P(E,,r,r’), the slowing-down 
kernel in the reactor of volume V,, and P.(E,, 
|r—r’|), the slowing-down kernel in a system of 
similar composition but of infinite extent. If the 
reactor is very large in comparison to the root 
mean square slowing down distance 


4 
frre, r)rdr 
Tae 


fr. (E, r)r’dr 


then the central portion of such a reactor will 
hardly ‘‘feel” the effects of the boundary. For 
such points, far from the reactor boundary, the 
so-called ‘“‘asymptotic’’ pile equation will hold; 
viz., 


DA®,—<aaP; 


+0; [ 8.0) Po(Es |r—r’|)dr’=0. (2) 
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Since P..(E., |r—r’|) falls off rapidly as |r—r’| 
exceeds ¥ (r?)w, the extension of the integral to 
all space introduces practically no error. 

The exact functional form of P..(E, 7) can best 
be determined experimentally. Various theoret- 
ical expressions for P..(E, r) have been proposed, 
all of them being approximations in greater or 
lesser detail to the solution of the exact energy 
dependent Boltzmann transport equation. The 
most useful approximation, called the ‘‘age”’ 
theory approximation, is 


—?r/4r(E 
wa 2 fe 
(47) 


where 7(£), called the ‘‘age’’ of the neutron, is 
one-sixth the mean squared distance the neutron 
travels in slowing down to energy E. A formula 
for t as a function of E is given in Eq. (16b) 
although such formulas are only approximations. 

The age-theory kernel is the solution of the 
differential equation 


oP 
to.—(E, r) = DAP(E, ¥}, (3a) 
ou 


where u=InEo/E, called the lethargy of the neu- 
tron; £=average loss in lethargy per collision; 
o,=scattering cross section; and D=diffusion 
coefficient. Equation (3a) can roughly be under- 
stood as follows. From the definition of &, du/é 
is the number of collisions required to cause a 
lethargy change du. The time dé for unit flux to 
suffer this letha:zgy change is 


dt=du/(éos),. 
hence 
0 dua 0 


= éo,—. 
ot dt du Ou 


If the neutrons in the course of their slowing 
down are supposed to diffuse according to the 
usual time dependent diffusion equation, 


dP/dt=DAP, 
then in lethargy they will diffuse according to 
to,(0P/du) = DAP, 


which is Eq. (3a). Of course, it is assumed here 
that P, the probability that a neutron will appear 
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at energy E and point r, is proportional to the 
flux of energy E and r. 

Another formula for P.(E,7r) which has the 
merit of analytic simplicity is the so-called ‘‘one 
group” formula 


P.(E, r) =e—"!4/4arL?, (3b) 


where again L? is one-sixth the mean square dis- 
tance the neutron travels in slowing down to E 
and is to be determined experimentally. An ex- 
ponential formula such as Eq. (3b) is just the 
distribution around a point source of mono- 
energetic neutrons diffusing without energy loss. 
It is, of course, the Green’s function for the time 
independent diffusion equation, 


AP—P/L?=0. 


Thus Eq. (3b) implies a physical picture in which 
the neutrons are supposed to diffuse without 
energy loss until they have wandered so far from 
their place of origin that their mean square 
diffusion distance is the experimentally found 
L’; at that time they suddenly lose energy in a 
single jump down to the value E. 

A natural generalization of Eq. (3b) is to 
consider the slowing down to occur in two stages 
with sudden energy jumps at the end of each 
stage. For this case the slowing-down kernel is 
characterized by two areas, L,;? and L,?, each 
measuring the mean square distance of travel 
while slowing through one energy group. The 
slowing-down kernel in this “‘two group” case is 


P..(E, r) = {L e771! 41 — L2e-7/ 22}, 


4nr(L,4—L-'*) 


Finally, if the slowing down is assumed to occur 
in many infinitely small increments, the distribu- 
tion will reduce to a Gaussian like Eq. (3). The 
physical model leading to the age slowing down 
kernel is therefore sometimes referred to as the 
“continuous” slowing down model. 

Now a solution of the integro-differential Eq. 
(2) is very readily found to be a solution of the 
wave equation 


A®,+ B’6,=0, (4) 

provided B? satisfies the characteristic equation 
kPo(E., B’) 

~—ninnnnnennn iE, (S) 


1+L?B 
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where L=(D/o.)* is called the diffusion length 
and its square is one-sixth the mean square dis- 
tance a thermal neutron travels before capture; 
k=voy;/cq is the multiplication factor, and 


P.(E, B’) = f P.(E, 1)'e®-*dr 


=4r f "P(E, r)(sinBr/Br)r’dr (6) 


is the three-dimensional Fourier transform of the 
slowing-down kernel. Furthermore, the number 
of neutrons passing energy E per cm’ per second 


at r, called the slowing-down density and denoted 
by g.(E, r), is 


Qu(E, t) =ko®,(r)P..(E, B*). (7) 


Equation (7) states that the slowing down 
density at all energies is everywhere proportional 
to the slow neutron flux. This statement can be 
proved from the definition 


qo(E, 1) = hoa f ®.(0))PalE, |r—r’|)dr’ (8) 


oa 


by applying the theorem which is often used in 
quantum-mechanical scattering problems that if 
®, is a solution of the wave equation (4), then 
Eqs. (8) and (7) are equivalent. 

The results embodied in Eqs. (4) through (8) 
are so basic to nuclear reactor theory that they 
have come to be known as the fundamental 
theorem of pile theory. 

The fundamental theorem gives the ‘‘asymp- 
totic,” i.e., far from any boundaries, energy, and 
spatial distribution of neutrons in a bare ther- 
mal neutron reactor. It is a remarkable fact 
that the distribution is separable into a space 
and an energy part (Eq. (7)), i.e., the asymptotic 
spatial distribution of neutrons is the same for 
neutrons of all energies and is simply given by 
solutions of Eq. (4). Thus for a sphere, 


sinBr 
®, ~ dea (E, r) — 
Br 


X function of energy, 


while for a cylinder 


®,~qx(E, r)~Jo(Br) Xfunction of energy. 


ALVIN M. WEINBERG 


The probability that a neutron will not leak 
out of the system from the time it is born at 
energy Eo to the time it reaches energy E is 
evidently the ratio of the total number of neu- 
trons crossing energy E per second to the total 
number being born at energy Ep per second. This 
ratio can be shown from Eq. (7) and from the 
fact that P,,(Eo, B?) =1 to be 


f Qu(E, r)dr 


ce 


=P,(E, B?). (9) 
da (Eo, r)dr 


It is simply the Fourier transform of the slowing- 
down kernel! 

The asymptotic theory developed here actu- 
ally has validity in a bare reactor much closer 
to the boundary of the reactor than is at first 
apparent. Along with the finite pile Eq. (1) 
there are the boundary values 


@=0, g=0, 


which the neutron flux must extrapolate to on 
what is called in transport theory the extrapo- 
lated boundary. Actually, the boundary condition 
at the edge of a chain reactor is simply that no 
neutrons return once they leave the reactor. 
This in effect is a requirement on the angular dis- 
tribution of neutrons, viz., those neutrons travel- 
ing toward the reactor vanish at the boundary. 
Since simple diffusion theory on the whole 
makes no statements regarding the angular dis- 
tribution of the neutrons, transport theory which 
does consider thé neutron angular distribution 
must be invoked to define the boundary condi- 
tion. It is found in transport theory that the 
above-quoted boundary condition is a good 
approximation to the true state of affairs. 

The extrapolated boundary is about 0.71), 
beyond the actual boundary, \, being the “‘trans- 
port” mean free path in the reactor, i.e., the 
scattering mean free path corrected for any 
asymmetry in the angular distribution of the 
scattered neutrons. We shall assume the ex- 
trapolation distance is energy independent. This 
is a fairly good assumption in materials such as 
graphite, but is a poor assumption in water. In 
any case, if the reactor is very large, the extrapo- 











lation di:stance will hardly change the effective 
size of the reactor. Now the finite slowing-down 
kernel, P(E,r,r’) can in all practical cases be 
considered to be a linear superposition of infinite 
system slowing down kernels P..(E, |r—r’|), 
provided both finite and infinite kernels satisfy 
the same linear equation. For example, if the 
moderation is assumed to follow the Gaussian 


exp[— |r—1"|*/47(E)] 
[4xr(E)}! 


P..(E, |r—r’|) = 


where 7(£) is the Fermi age, then P(£,r, r’) 
will be that linear superposition of Gaussian 
functions with varying r’ that vanishes on the 
extrapolated boundary r,. Such linear super- 
positions are exactly analogous to the image 
system used in constructing a solution to a 
partial differential equation which satisfies pre- 
assigned boundary conditions from the solution 
in the infinite system. 

More exactly it can be stated that if ®, is a 
solution of the wave equation which vanishes 
on the extrapolated boundary, and if P(£,r, r’) 


and P.(E, |r—r’|) satisfy the same linear equa- 
tion, then 


f ®,(r’) P(E, r, r’)dr’ 


ce 


= fe), |r—r’|)dr’. (9a) 


o 


Thus the asymptotic pile Eq. (2) is equivalent 
to the finite pile Eq. (1), and the asymptotic 
solution solves the finite equation all the way to 
the boundary except for transport theory devia- 
tions which set in within a mean free path or so 
of the boundary. 

The sense of this theorem is somewhat as 
follows: the neutron distribution in a finite pile 
can be calculated (under the restrictions stated) 
by extending the pile out to infinity and finding 
the asymptotic neutron distribution in the in- 
finite system. This solution oscillates, positive 
neutron densities alternating with negative ones 
ad infinitum. The positive and negative densities 
are so distributed that on the pile boundary their 
superposed effect satisfies the boundary condi- 
tions. 


The alternating ‘positive’ and ‘“‘negative”’ 
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TaBLE I. Slowing down models, Fourier transforms, 
and characteristic equations. 


Characteristic 


Model P.(Es, 7) P.o(Es, B?) equation 
All-thermal 4(r) 1 B?=(k—1)/L? 

2 Group, elk, 1 k 

thermal, Sens aoe ——— | 
nonthermal 477rLi* 1+L::B? (1 +L2B?) (1 +112B?) 


or 





oe go ae B ko tan"B/oe = 
Transport a" B tan = Bat+LBy 1 
Gaussian exp( —r?/4r) k exp( —7rB?) 
or “age” Gt exp( —7B*) 1+L*B? =i 


piles constitute a system of images such as are 
often used to solve ordinary boundary value 
problems. The beauty of the theorem embodied 
in Eq. (9a) lies in the fact that even for an 
arbitrarily-shaped reactor the appropriate image 
system is shown to be that which is generated 
by the analytic continuation of the asymptotic 
neutron distribution outside the pile. 

The recipe for computing the critical size of 
the reactor is now clear: the critical boundary 
simply coincides with the surface defined by 


®,(r., B) =0, (10) 


where ®, satisfies the wave equation (4) and B? 
is determined as the solution of the characteristic 
Eq. (5). For example, for a sphere the critical 
radius r, is determined by Eq. (10) from the 
relation 


sinBr./Br.=0, r.=7/B. (10a) 


The quantity B found as a solution of Eq. (10) 
is a purely geometric property of the critical 
system and is called the geometric buckling of the 
reactor; the B found as a solution of Eq. (5) is 
determined entirely by the intrinsic properties 
of the reactor—its composition and the nuclear 
characteristics of its constituents. It is therefore 
called the material buckling of the reactor. The 
critical condition can therefore be paraphrased 
as follows: a reactor is critical if and only «f its 
material and geometric buckling are equal. 

The explicit form of the characteristic Eq. (5) 
will naturally depend on the model used to de- 
termine P,.(E, |r—r’|). Table I lists some of the 
models that have been used and the correspond- 
ing critical equations. 

In the all-thermal model, all neutrons are 
supposed to be born at thermal energy. The 
slowing-down kernel therefore degenerates into 
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a 6-function since the neutrons become thermal 
exactly at their birthplace. 

In the two-group model, all fast neutrons are 
lumped together into one group, and all ther- 
mals are lumped into another. This model is very 
extensively used for quick, order-of-magnitude 
calculations. 

In the transport model, fast neutrons are sup- 
posed to be degraded into thermal neutrons by a 
single collision characterized by a cross section 
o. The ‘transport’? model is so called because 
oe~°’/4nr? is just the kernel used in the Boltz- 
mann transport equation. It has only theoretical 
significance, although it very roughly approxi- 
mates the situation in HO where a single colli- 
sion with H usually brings the neutron down to 
near thermal energy. 

The Gaussian model, which has already been 
discussed, is the most accurate for reactors 
moderated by elements heavier than hydrogen. 

It is an important fact that the most general 
characteristic equation for a thermal fission 
reactor 

kP..(E., B?) 


=1 (S) 
1+ L?B? 
always reduces to 


k-1 
PP a—__—_, (10b) 
L?+ (r*)n,/6 


provided & is close to unity. This can be seen by 
expanding P..(E,, B?) in Eq. (5) in powers of 
B? and keeping only the first two terms. Hence, 
whenever & is near unity, the critical radius of, 
for example, a sphere is always by Eq. (10a) 


aM 


r-=— ’ (10c) 
v (k—1) 


where M, called the migration length, is 
M =v (L?+(r?)/6). 


The migration length determines the length 
scale for a large chain reactor. The “number of 
M’s”’ in the critical radius is, by Eq. (10c), just 
a/V (k—1). If R>1, a real r, will exist and the 
system will chain react, while if k<1, 7. will be 
imaginary and the system will not be self- 
sustaining, no matter how large it is. 
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RESONANCE REACTORS 


Actually the asymptotic theory of even a 
reactor with resonance neutron as well as ther- 
mal neutron fission can be written down. To 
accomplish this we generalize the continuous 
slowing-down Eq. (3a) to include multiplica- 
tion by thermal and nonthermal fission. Thus 
the most general uniform reactor equation with 
resonance fission, and continuous slowing down, 
is 


D(u) 
Ao .&(r, u) —oab(r, u)+ F(u) 


o,(u 


| fearon, w')du! + (aoa), 


0 
=—io,P(r, u), 
Ou 


(11) 
where u is the lethargy; F(u)du is the distribu- 
tion of fission neutrons in lethargy; uo is the 
lethargy corresponding to thermal energy; k(u) 
=vo;(u)/oa(u); ££ the logarithmic decrement, 
=1+[a?/(1—a?) Ine’; a? =[(M—1)/(M+1)}; 
and M=atomic weight of the moderator. The 
thermal neutrons obey Eq. (11) with €=0 and 
the source term replaced by q(uo). 

The first two terms on the left side of Eq. 
(11) represent diffusion and absorption. The 
next two terms represent production of neutrons, 
first by nonthermal fission, second by thermal 
fission. Finally the term on the right hand repre- 
sents the net rate of change in neutron popula- 
tion at lethargy wu as the result of elastic scatter- 
ing and moderation of higher energy neutrons. 

Now the asymptotic solution of Eq. (11) is 


@(r, u) =y(r)X(u) (12) 
provided 


Ay+ By =0, 


“Ca B’D 
0 oY 


xf Pero] f es 


x js} (14) 


and B? satisfies the following generalized char- 


(13) 
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acteristic equation: 


[Pa 


b(z, w)exp[ — B?r(z, w) ] 
os w) 











R(uo) me 
X F(z)dedw+———— f P(2, Uo) 
1+ L?B*J 


Xexp[ — B’r(z, uo) ]F(z)dz=1, (15) 


where p(z,w) and 7(z,w) are defined in Eggs. 
(16a) and (16b); that is, even for a general re- 
actor in which fissions occur at arbitrary en- 
ergies, the asymptotic flux distribution is sepa- 
rable in space and energy, and the power ‘dis- 
tribution is a solution of the wave equation, 
provided the reactor is uniform and the extrapo- 
lation distance is energy independent. In other 
words, the fundamental theorem holds for reso- 
nance, as well as thermal, reactors. 
The characteristic Eq. (15) can be rewritten 


ky; exp(— B?r;) 


















+k, exp(—B?r,)/(1+L?B?)=1, (16) 
where ky is the so-called “‘fast’’ multiplication 
constant, k, is the “thermal” multiplication 


constant, and 7; and rf, are the ages of neutrons 
which cause nonthermal and thermal fissions, 
respectively. Comparison of Eqs. (15) and (16) 
shows that 


lise ; | —_— ph w)Fle) ded 


_ 204(uo) 
















b (3, ue) F(2)ds 


Oa(Uo) 


sto fo 


Ta(W) 
vote, w)exp[ — B’r(z, w) ]F(z)dzdw 













Eo,(w 


Inn fp ) 
ee . a P(Z, Uo 


Xexp[ — B’r(z, uo) |F(2)dz. 





Herein p(z, w) is the probability that a neutron 
will escape resonance capture in slowing from 
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lethargy z to lethargy w, r(z, w) is the age of a 
neutron slowing from z to w. These functions 
are given by 





” oa(u) 
; = - d 16. 
p(s, w) exo| f Saez uw] (46a) 
» D(u) 
7(%, w) = ——du (16b) 
a (ue 


Equation (16) is a natural generalization of the 
usual characteristic equation (line 4 in Table 1) 


for a thermal fission reactor; it reduces to it 
when k;=0. 


MULTIGROUP METHODS 


The Fourier transform theory described above 
is beautifully elegant since it includes as special 
cases all the usually considered slowing-down 
models. While special forms of the characteristic 
equation were found in project literature from 
the start, it was not until after the war that the 
general formulation, particularly the relation 
between the asymptotic and finite reactor equa- 
tions was worked out. 

In spite of its elegance, the results of the 
Fourier transforms theory are of little use in 
most reactor systems of current importance. 
The reason is that nowadays interest in many 
cases centers on small, enriched reactors which 
are surrounded by reflectors. Thus the nuclear 
characteristics of the system are not independent 
of position, and therefore the asymptotic theory 
is applicable at best at the center of the reactor, 
or, if the reactor is small, not even there. 

In a reflected reactor the neutron energy spec- 
trum evidently varies from point to point: if 
the reflector is a moderator, the spectrum will 
be progressively degraded from reflector-reactor 
interface out to the reflector edge. As the re- 
flector-reactor interface is crossed going towards 
the center of the reactor, the spectrum tends to 
approach the asymptotic spectrum, which in a 
thermal reactor is by Eq. (7) 

qu(E, t) =kou®,(r)P..(E, B?). (17) 
The spectrum will in fact reach this asymptotic 
form only if the reflector-reactor interface is 
sufficiently far from the center of the reactor. 
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Thus, in addition to the asymptotic solutions 
of the reactor equation there are in general non- 
asymptotic solutions which die off as the center 
of the reactor is approached and which, when 
appropriately superposed, describe the neutron 
energy spectrum throughout the reactor. These 
nonasymptotic solutions correspond to the com- 
plex eigenvalues of the basic characteristic Eq. 
(5). The exact form of the nonasymptotic solu- 
tion depends on the form of the slowing-down 
kernel. 

The scheme which has been generally adopted 
for handling such small reactor-reflector systems, 
especially reactors in which considerable ab- 
sorption and fission occurs at several nonthermal 
energies, is the previously discussed multigroup 
method. In this scheme neutrons are assumed to 
pass discontinuously from one energy interval or 
“group” to the next, the neutrons diffusing 
without energy loss until they have suffered a 
number of collisions which on the average are 
sufficient to degrade them out of their energy 
group. The best way to rationalize the group 
method is, as was pointed out by Hurwitz and 
Brooks, to start with the age-diffusion equation 
in the form 

DA®—o.6+Q(r, u) = d(éo.&)/du, (18) 
where Q is the number of neutrons per second 
per cm* being produced by fission at (r, ~). In 
Eq. (18) the first term represents the number of 
neutrons lost per cm® per second per lethargy 
interval by diffusion ; the second term the number 
lost by absorption, and the third the number pro- 
duced by fission. All this is equated to 0(to,@) /du, 
the rate of change of neutrons per cm’ per second 
per lethargy unit. Equation (18) is fairly ac- 
curate if the absorption is weak and the modera- 
tor is not too heavy. The equations of multi- 
group theory are obtained by integrating Eq. 
(18) with respect to u over the logarithmic 
energy range, U;, corresponding to the i-th group. 
Then, upon replacing averages of products by 
products.of averages, there results 


D,A®; —oa,P;= + Sa;+b;* = Saiz - ve Q;(r) ’ 
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where 


1 pet 
ieee f &(u)du 
U.S, 


1 ut 
Sissinwe f D(u)du 
Ud, 


1 ut 
oa; =— oa(u)du 
a. 


bie fo,(ut) ' ss to,(u-) 


U; U; 


@;7= @(ut) ’ 6-= }(u-) 


1 pt 
Qe f QC, wa 


and ut and u~ denote the upper and lower 
logarithmic energy boundary of the z-th group. 

The term Sa;+é;+ represents the neutrons 
which leave the i-th group by slowing down; the 
term Sa;-®;— represents the neutrons which enter 
the i-th group by degradation from a higher 
group. 

The best relation between ®;, @;+, and ®;>- 
has been a matter of much discussion. In the 
earliest literature it had always been assumed 
that 


@+=6;; O-=9,1. 


However, other somewhat more general assump- 
tions, such as 


;=37(O;++;), 
@;+=6;_)- 


have also been proposed. In any case, the group 
equations always form a coupled set of m dif- 
ferential equations essentially of the form 


D,A®; - oa;P; — SaP;+ Saj_1P;_1 


+F; Dd ikjoasb; =0, (19) 


where k; in the last term is the number of neu- 
trons produced per neutron absorbed in energy 
group j; and F; is the fission neutron spectrum 
normalized to unity, i.e., F; is the probability 
that a fission neutron will have its energy in the 
a-th group. It is assumed that the fission spec- 
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trum is independent of the energy at which 
fissions occur. 

The multigroup equations, being essentially 
approximations to the most general reactor Eq. 
(11), must have a separable asymptotic solution, 
and this asymptotic solution must be the true 
solution in case the reactor is bare and uniform 
and the extrapolation distance is energy inde- 
pendent. For in this case the assumption 


A®;+ B46; = 0 


for all 7, leads when substituted into Eq. (19) to 
the linear homogeneous algebraic equations 


— (D,B*+ oa;+ Sa;)®;+ Saj;_1Pj_1 
+F; Dijkjoajb;=0, 
which have a nontrivial solution if and only if 
det | — (D;B?+0;+ Sa;) + Saj_-1 
+ F; > jkjoa;| =0. 


This is the characteristic equation in the multi- 
group method which corresponds to the most 
general characteristic Eq. (16). 

In the simple two-group picture with only 
thermal fissions occurring, the characteristic 
equation is given in Table I as 


k 
secinnineinsinenecapmnpemmenn ian § 
(1+L,°B*)(1+L2?B?) 


L;=D;/(oa;+Sa;). The quantity (1 
+L,?B?)— is of course just the Fourier trans- 
form of e~"/42/4rL.?r, which is the slowing down 
kernel resultin, from the one-group picture. The 
two solutions of Eq. (21) are 


1 1 
toa) 
L? L? 


1 1 
Get 
LL’ L? 


which corresponds to the asymptotic solution 
since B,? is real and 


1 1 1 
Bt=—| —(—4+— 
Z L;? L? 
1 2 


1 
Mitt 
L? L? 


(20) 


(21) 


where 


4(k—1)}! 
+————]}, aa 
LYL? 


4(k—1)7 
+] , (2b) 
L;*L? 
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which is complex and therefore corresponds to a 
complex By and to an eigenfunction which is 
damped; i.e., nonasymptotic. If k is very near 
unity, ByY=(k—1)/(L2+L22), which illustrates 
the fact that for a large reactor the buckling is 
given by Eq. (10b). 

The great convenience of the group equations 
is that the energy does not enter explicitly. 
Thus even a system involving indefinitely many 
reflectors and with arbitrary cross-section varia- 
tion can be solved if the set of Eqs. (19) can be 
handled. These, being coupled wave equations, 
are amenable to numerous computational tech- 
niques. 


PERTURBATION THEORY—THE IMPORTANCE 
FUNCTION 


The time behavior of a chain reactor which is 
not critical, or, being critical, has an original 
neutron distribution which is not the steady- 
state distribution is calculable from the original 
nonsteady-state equations. The basic time de- 
pendent equation in the multigroup picture is 
obtained by equating the left-hand side of Eq. 
(19) to d@;/dt. This equation will then describe 
the time behavior of the reactor to the extent 
that the delayed neutrons are ignored. In opera- 
tor notation, the equations are 

> Hib; = 06; /dt t=1-- “nN, 


(23) 
where Hi; is the matrix differential operator 
Hi; = ||DA eae (oa;+ Sa) oe Sajat F; Dd ik joa;|| 


and ®; is a vector whose components are the 
neutron densities in the ” groups. (The extra 
velocity factors are incorporated in the quanti- 
ties D;, oa;, and Sa;.) The problem is,,given an 
arbitrary initial distribution @;=0;(r), how does 
®; approach its steady-state distribution if the 
reactor is critical, or how does ®; reach its steady 
mode of increase or decrease if the reactor is 
not critical? 

To answer such questions it is necessary to 
develop an expansion theorem for the reactor 
equation (23); i.e., it is necessary to find a set 
of functions which are characteristic functions 
with characteristic values \, for the operator 
Hi; and in terms of which arbitrary functions 
which satisfy the reactor boundary conditions 
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can be expanded. Now because the matrix H is 
not symmetric—i.e., the cross section for a neu- 
tron to go from group 7 to group 7 is not the same 
as the cross section for a neutron to go from group 
j to group 7, the H operator is not in general 
self-adjoint. In this respect the mathematical 
theory of chain reactors is unusual among physi- 
cal theories: most such theories, e.g., Schroed- 
inger’s wave mechanics, being derivable from a 
variation principle involving a quadratic differ- 
ential form, are self-adjoint. 

Mathematically, this means that the functions 
f;® which are characteristic for 


Vi fj=ufi™ 
(24) 
f:;®” =0 on extrapolated boundary 


do not form an orthogonal set; rather they form 
a bi-orthogonal set with their adjoint functions 
gi, the adjoint functions being defined to 
satisfy 


De idij*gi = gi” 
(25) 
g;” =0 on extrapolated boundary, 


and the adjoint operator Hi;* is so chosen that 
Digi fj= Vy fii;*g;, 


where )0i,g:f7i;f; is a notation for the scalar 
product of the vector g; and the vector > ;Hi;f;; 
i.e., 


Cugufi= Le [ele ulnsf(eddr, 


The adjoint operator matrix Hi;* is obtained 
from Hi; simply by transposing. The bi-ortho- 
gonality condition is 


DVifig: =5,,, 


where 6,,=1 if v=y, 6,,=0 if vAp. 

The simplest example, and one which brings 
out the physical significance of the adjoint’ func- 
tion, g:, is a just-critical, steadily running 
reactor in which the neutron distribution is in 
its fundamental mode, ®;= f;(r), and to which 
at time ¢=0, a single neutron is added in the 
j-th energy group at a point r’. Then 9,(r, #) 
must satisfy the time dependent differential 


(26) 


* The term adjoint function, when used without modifier, 
means the lowest eigenfunction of the adjoint equations, 
i.e., 25. 


ALVIN M. WEINBERG 


Eq. (24) and at time ¢=0, the components of 
® must satisfy 


®(r, 0) f(r) iA 


except for the j-th group neutron density which 
must satisfy 


$,(r, 0) f; (r) + 6(r—1’), 


6(r—r’) being a Dirac delta-function. 
To compute each component ®,(r,7¢), we 
assume 


(1, )=DAfiM (re i=1---n (27) 


and substitute into Eq. (23). Then each f;‘”(r) 
is a characteristic function which satisfies Eq. 
(24); and upon applying the bi-orthogonality 
condition, Eq. (26), to the initial distribution 
®,(r, 0) there results 


Bilt, t) =f (0) +g) (0) f() 
+E BC) FO(WMO. (28) 


Now \Ao=0 because the pile is just critical; the 
other eigenvalues A, can all be shown to be 
negative. Consequently the neutron distribution 
eventually approaches 


Bil, «i= fi () +. gj (0) F.0(4) 
= fo) [1 +.95(C) 


Thus a neutron:of energy group j added at r’ to 
a just-critical chain reactor increases the total 
neutron population by a fractional amount 
g;(r’). This ultimate increase in neutron popu- 
lation is the result of successive generations of 
fissions induced by the originally introduced 
neutron. The total number of fissions caused in 
the n-th generation following introduction of a 
j-th energy neutron at r’, as m approaches o, is 
a function of 7 and r’ called the iterated fission 
probability. Since the number of fissions is pro- 
portional to the number of neutrons, the iterated 
fission probability is proportional to the ulti- 
mate increase in neutron population induced by 
introduction of the original neutron—i.e., to the 
adjoint function. In this sense, the adjoint func- 
tion g;(r’) is a measure of how “important” a 
neutron of given energy introduced at r’ is in 
sustaining the chain reaction, and therefore it 


(29) 


















has bec.»me customary to call the adjoint func- 
tion the importance function. 

The importance function depends on the en- 
ergy and on position in the reactor. Thus a neu- 
tron introduced at the outer edges of the reactor 
has a good chance of escaping without causing 
any fissioning at all, and it is therefore not very 
important to the chain reaction. For this reason 
the importance function is smallest at the outer 
edges of a chain-reacting system. 

As is usual in quantum mechanics, so in nu- 
clear reactor theory, the expansion theorem is 
the basis for perturbation theory. Suppose a 
system which is just-critical is subjected to a 
perturbation—say a small change in a cross 
section. The change in cross section will make 
the reactor either sub- or super-critical, and the 
neutron intensity will change in a manner some- 
what like that given in Eq. (27). Eventually 
only a single exponential will survive; i.e., the 
reactor power will rise or fall exponentially. The 
problem solved by perturbation theory is the 
determination of this steady-state period, or 
what amounts to the same thing, the effective 
deviation of the multiplication factor from criti- 
cality caused by this perturbation. This devia- 
tion can be considered equivalent to a change 
in the number of neutrons per fission, if desired. 

In the original system 


Vif; (x) =0 
and the adjoint system satisfies 


Digi 1) At; =0. 


t=1---n, 


(30) 


The perturbed system, characterized by an 
operator H:i;+ 6Hi; will, in general, not be critical. 
It ultimately will rise or fall like e**; hence, the 
perturbed distribution f;’, must satisfy 


Do i(Hij+ 6H) Cf Je =a fi Je. (31) 


Evidently, by application of the bi-orthogonality 
relations, there results in first order 


A= Lig 5M fi) (32) 
or, since f;’ and f; differ only in the order of 
5o0a;, which is assumed to be small, f;’/ can be 
put equal to f; in Eq. (32): 


A= LegiM sy f,. (32a) 
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Now ) is the rate at which the neutron density 
rises per second. The rate at which the neutron 
density rises per generation is \/T, T being the 
time required for a neutron to go through a com- 
plete multiplication cycle or generation. The 
deviation of the multiplication constant from 
criticality is by definition the rise rate per 
generation, and therefore \/T gives the devia- 
tion caused by the introduction of the per- 
turbation. 

As an application of the foregoing general 
theory, we consider a two-group model reactor 
in which the cross sections are allowed to vary 
with position. If the reactor were critical, the 
neutron distributions in the two energy groups, 
fi(r), fer), will satisfy the Eqs. (19) or 
(30), with the D;A operator replaced by V-D.V 
because D; is space dependent, namely 


(V+ DV —oa,— Say+ Fyki00) fy (r) 
+ Fiksoaf2(r) =0 
(Say+ Foki001) fx (tr) + (V-D2V —oa2— Suz 
+ Fokooas) fo (r) =0. 


Since group (1) is at the top of the energy scale, 
neutrons from group (2) cannot slow into group 
(1); for this reason no term in Sa, appears in 
the first equation. 

Now suppose the capture cross section in 
group (1) is changed from oa to oa+6ea, but 
the multiplication constants are unchanged. The 


new distribution will obey a set of equations 
like Eq. (31): 


[V- DV —oa;— Sait Fikioa 
— 00,(1 — Fikx) ]fi’ (x) 
+ Fykooaefe!’ (r) =rfi! Xr) 
(Sai+ Fokioai+ do, F 2k1) fy’ (r) 


a (V > DNV —oa2.— Sao+ Fakeoes fa Sy, 


(33) 


(34) 


The importance functions satisfy the trans- 
posed equation: 


(V- DiV —oa,— Say+ Fyky00;)g1(r) 
+ (Say-+ Fokioa1) g9(r) =0 

F kooa2g(r) + (V+ D2V —oag— Sez 

4+ Fokeoas)g2(r) =0. 


(35) 
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Hence, if Eqs. (34) are multiplied on the left by the vector [g:(r), go(r)], Eqs. (35) are 
multiplied on the left by [f1’(r), fe’ (r)], and the results subtracted, we obtain 


Afr! gy + fe! go) = g,OV- Di V fy — fy’ OV - DV gi + g.V - DV fo! 


- fe! OV -DVg2 —_ 6o0;(1 = Fyki) fy) gs + boa, F oki fy! go, 


(36) 


Upon integrating Eq. (36) over the reactor, the differential terms disappear by Green’s theorem and 
the vanishing of f;’“ and g; on the boundary. Hence 


— f soap @ — Fyky)g1 (r) — Fokige (r) Jdr 


A= 


PCA g1+ fg 


and if the unperturbed and perturbed distribu- 
tions f; and f,;/ are set equal to each other, 
the final result is 


d= = f boosh) LU Fda) 
— Frokige(r) ]dr (38) 


since the denominator of Eq. (37) is normalized 
to unity. 

The sense of Eq. (38) is the following: 
60a;f;°(r) is the additional number of neutrons 
absorbed per second per cm? at r because of the 
increased oa;. The net loss to the chain reaction 
of group (1) neutrons is doa:f;(r)(1— Fk) and 
the net gain of group (2) neutrons resulting from 
the additional fissions implied in an increase of 
ga; is b00;f;("(r) F2k:. Now the importance of a 
group (1) neutron at r is g(r) ; the importance 
of a group (2) neutron at r is g2(r). Conse- 
quently the integral in Eq. (38) is the sum of all 
the neutrons lost to the chain reaction per sec- 
ond, weighted with their importance in sustain- 
ing the reaction, and it is plausible that this 
should be the time constant at which the in- 
tensity of the chain reaction falls off. 

The computation of the importance function 
goes exactly like the computation of the dis- 
tribution itself since Eqs. (33) and (35) are of 
the same form. In the case of a bare, uniform 
reactor the spatial distribution for all neutron 
energies isthe same. In this case g; is propor- 


tional to f;, and by Eq. (32a), or more specifi- 
cally Eq. (38), the effectiveness of an absorber 
placed at r in a chain reactor in depressing the 
reaction is proportional to the square of the 
neutron flux at r. For example, in a cube of edge 
HI the neutron distribution is 


™x TY 72 
sin— sin— sin—. 


An absorber placed at the center of a bare, uni- 
form cubical reactor therefore has the same 
effect on the criticality of the chain reaction as 
8 times as much absorber distributed uniformly 
throughout the reactor. 

The foregoing account of reactor theory is of 
course not intended to be complete. However, 
I hope that I have given enough to enable the 
reader to appreciate some of the elegance which 
underlies the theory of chain reacting systems. 
The analogies with quantum-mechanical mathe- 
matics are extraordinary—in fact, the Schroed- 
inger equation can be thought of as the equation 
for an infinitely large, just-critical system in 
which the neutron absorption is distributed like 
the potential V(r) and in which the slowing-down 
kernel is a delta-function. It may well be that 
instructors in mathematical physics will find the 
theory of multiplying chain reactors an unusu- 
ally timely and interesting addition to the stand- 
ard list of problems which are usually covered in 
such courses. 





Reproductions of Prints, Drawings and Paintings of Interest in the History of Physics 
47. The Heyday of Electrostatic Experimentation 


E. C. Watson 
California Institute of Technology, Pasadena 4, California 


(Received April 21, 1952) 


Two contemporary engravings showing early electrical experiments are reproduced together 
with an entertaining account of such experimentation written by Park Benjamin in 1895. 


ARK BENJAMAN in his The Intellectual Rise 

in Electricity: gives the following entertain- 

ing account of the popular interest in electrical 

experimentation that sprang up about the middle 
of the eighteenth century: 

“The year 1742 probably marks the beginning 
of the singular and sudden interest in things 
electrical which arose in Germany, and which 
swiftly reached a stage of feverish enthusiasm. It 
differed widely from the perfunctory craze which 
had taken possession of the English aristocracy 
at the behest of Charles. It had still less resem- 
blance to the combined onslaught of the French 
philosophers which was designed to take all of 
Nature’s secrets by storm. It was distinctly 
popular. It was the first instance—many times 
since repeated—of the intelligent portion of an 
entire community regarding with absorbing 
wonder the working of electric powers. 

“No unexpected desire for electrical knowledge 
in general had been born. The German natural- 
ists were familiar with progress abroad during 
the last fifty vears, but had shown no emulative 
spirit. The new motive force now came not from 
them, but from the people; and the people, in 
all times and in all ages, have never failed to re- 
spond to an appeal to their sense of the mar- 
velous—to a conviction that something new has 
been found—something at once new and incom- 
prehensible. The masses had cared little for 
Hauksbee’s lights, and less for the vagrant virtue 
on Gray’s lines, assuming that the knowledge of 
either percolated to them ; but when it came to be 
noised about that the strange radiance which the 
English and French philosophers were exhibiting 
was fire,—fire which flamed in jets from the ends 


1. Park {Benjamin, The Intellectual Rise in Electricity 
(London, 1895), pp. 491-495. 


of rods, or, more wondrous still, leaped from the 
tips of men’s fingers—that was a matter for 
everyone’s personal concern. For fire was then 
believed to be a material substance—phlogiston 
—and while perhaps it might exist in iron bars 
and inanimate things of that kind, and be forced 
visibly to come out of them by friction, as well 
as by heating, no one had ever supposed that it 
resided in the human body and could be com- 
pelled to escape, with an accompaniment of 
sparks and crackles, from one’s person. It was 
the idea of a human being becoming such a 
torch that stirred the Teutonic mind to its pro- 
foundest depths. The impetus which electrical 
science had received from the fancy of a dissolute 
king was nearly spent: now progress was re- 
sumed with renewed vigor under that due to the 
astonishment and wonder which the latest elec- 
trical manifestations had created in the now 
thoroughly awakened Germans. 

*“*. . . It is doubtful to whom is due the credit 
of accomplishing the work which began the new 
era; some contemporary writers according it to 
Christian August Hansen, others to George 
Matthias Bose. .. . 

‘‘ . . Both Hansen and Bose found, at about 
the same time, that not only could a practically 
continuous supply of electricity be obtained 
[with a glass sphere mounted lathe fashion and 
rotated rapidly by means of a crank-wheel and 
belt], but one of much greater strength than had 
hitherto been known. Hansen suspended a boy 
with his toes in proximity to the globe, and drew 
sparks from his fingers. Bose disposed twenty 
soldiers in line, with hands touching and ad- 
ministered a shock to all of them at once. . . . No 
one knew better the art of playing to the gallery 
[than Bose ]; in fact, in the great electrical drama 
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Fic. 2. Electrical experimentation in the seventeen forties. 





HEYDAY OF ELECTROSTATIC EXPERIMENTATION 


he creatzd the part of the ‘modern wizard,’ and 
it is doubtful whether anyone since has ever ex- 
celled him in it. He set jets of fire streaming from 
electrified objects and exhibited them to the 
people who flocked to his laboratory. He invited 
guests to an elegant supper-table loaded with 
silver and glass and flowers and viands of every 
description, and, as they were about to regale 
themselves, caused them to stand transfixed with 
wonder at the sight of flames breaking forth from 
the dishes and the food and every object on the 
board. The table was insulated on pitch cakes, 
and received the discharge from the huge glass 
retort which was revolved in another room. He 
introduced his ardent pupils to a young woman of 
transcendent attractions and as they advanced 
to press her fair hand, a spark shot from it ac- 
companied by a shock which made them reel. 
Others, who had the boldness to accept his 
challenge to imprint a chaste salute upon the 
damsel’s lips, received therefrom a discharge 
which Bose says ‘broke their teeth,’ but Bose 
here either exaggerates more than usual, or else 
neglects to explain how the young lady bore her 
share of the injury.” 

By 1746 showmen were traveling even to 
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America with their electrical machines giving 
people shocks for a small fee. It was one of these 
itinerant electricians who aroused BENJAMIN 
FRANKLIN’S interest in the subject. 

The experiment of suspending a boy by silk 
lines (Fig. 1), which demonstrated the conduc- 
tivity of the human body, was first performed by 
STEPHEN GRAY in 1730. It was repeated by 
CHARLES FRANCOIS DE CISTERNAY DuFay, who 
was the first to suspend himself and describe his 
sensations. The account given by PARK BENja- 
MIN in The Intellectual Rise in Electricity of the 
work of Gray and Duray and of their relations 
is well worth reading not only because of the 
scientific importance of the work described but 
as an example of what the personal relationships 
between scientific workers, even of different 
nationalities, can and should be. 

The two plates here reproduced are from the 
French translation (Paris, 1748) of WILLIAM 
Watson’s Experiments and Observations tending to 
illustrate the Nature and Properties of Electricity. 
They show quite well not only the construction 
and operation of the electrostatic machines used, 
but also the types of experiments performed and 
the popular interest in them. 


Physics Today 


As a result of successful efforts to increase outside 
income, the American Institute of Physics will be able in 
1953 to send PHYSICS TODAY to all members. This 
service will be rendered without adding in any way to the 
financial burdens of the Societies or their members. It is 
being undertaken in response to a growing demand, and it 
accords with the recommendation of a Committee appointed 
to advise the Institute on its future course with respect to 
PHYSICS TODAY. 

There are now nearly 15,000 physicists in the Institute. 
A group so large, driven ever to greater specialization by 


the development of their field, can only with the greatest 
difficulty maintain unity and combine strength in advanc- 
ing commonly held objectives. The Institute was originally 
established by physicists to solve this very problem, and 
PHYSICS TODAY has come to be perhaps its most 
effective instrument. One medium of communication go- 
ing to all is vitally needed, if physicists are to gain the 
advantages of community action and meet the demands 
of professional responsibility. It is on these grounds that 
extraordinary measures seem justified to give PHYSICS 
TODA ¥ general circulation. It will begin with the January 
1953 issue 





William Gilbert on the Rainbow 


CarL B. BovER 
Brooklyn College, Brooklyn, New York 


(Received February 8, 1952) 


A qualitatively satisfactory explanation of the rainbow was given shortly before 1311 by 
Theodoric of Freiberg and Qutb al-din al-Shirazi, but this work had disappeared shortly before 
Gilbert was born. The quasi-Aristotelian theory which Gilbert espoused was akin to early medi- 
eval ideas according to which rays from the sun are tinged as they traverse a thin dewy vapor 
and then are reflected to the eye of an observer by a dark cloud or dense object. Gilbert had the 
happy thought that a spherical magnet might serve as a miniature earth, but he missed en- 
tirely the notion of a globe of water as a magnified raindrop. Unaware of the role of the spherical 
drops in the refraction of solar rays, the fuzzy speculations on the rainbow which Gilbert gave 
in De Mundo stand in marked contrast to the experimental philosophy which he advocated in 


De Magnete. 


EW scientists have earned so persistent a 

reputation as did William Gilbert (ca 1544— 
1603), physician to Queen Elizabeth I and King 
James I and father of the science of electricity. 
Dryden correctly predicted that 


Gilbert shall live till loadstones cease to 
draw, 
Or British fleets the boundless ocean awe. 


But a well-launched reputation all too easily 
grows beyond the limits of sober valuation. 
More often than not, in the case of great scien- 
tists, the truth that men write lives after them, 
the error oft is interred with their books. One 
is tempted to add, ‘‘So let it be with Gilbert,” 
were not truth, in the history of science, a 
nobler virtue than adulation. To illuminate 
only those ideas which eventually led on to 
further success, passing over in silence those 
which posterity has judged to be of no value, 
is to foster a totally erroneous impression of the 
way in which science has developed in the past— 
and, presumably, is likely to continue to grow in 
the future. This is not to say that one should 
admire less such a great book as the De Magnete 
Magneticisque Corporibus, Physiologia Nova;' it 
is simply to call attention to the fact that Gilbert 
wrote also the De Mundo Nostro Sublunari, 
Philosophia Nova,? a work of a different char- 


1 Published in 1600 at London. Translated into English 
by P. F. Mottelay (New York, 1893). 

2 Published in 1651 at Amsterdam. A brief account of 
the contents will be found in Lynn Thorndike, A History 
of Magic and Experimental Science (6 vols., New York, 
1923-1941), VI, pp. 61-63, 379-382. For an analysis of its 
influence on Francis Bacon, see Park Benjamin, The 
Intellectual Rise in Electricity (New York, 1895), pp. 
315-331. 


acter. For this reason attention will be focused 
here upon chapters x—x11I in Book IV of the 
latter volume, covering an aspect of Gilbert’s 
thought which appears to be virtually unknown 
today. 

This paper suggested itself when, some years 
ago, I ran across the following passage in a 
history of physics: 

“‘In optics he [Gilbert ] gave an explanation of 
the formation of the rainbow by the refraction 
of light in the drops of water, which recalls that 
of Maurolycus.’”* 

There is an inadvertency here, for Maurolycus 
did not attribute the rainbow to refraction; but 
the thought occurred to me that in Gilbert one 
might find, for the theory of the bow, a worthy 
precursor of Descartes. “‘In an age of empiricism 
and ignorance he rescued the study of the magnet 
from the atmosphere of occult mysticism with 
which it was surrounded.’ Could he not have done 
the same for the rainbow? Turning to the two 
books by Gilbert, one finds that De Magnete has 
nothing on the theory of the rainbow; but more 
than half a dozen pages of De Mundo are devoted 
to it.6 A large part of the latter volume, promis- 
ingly enough, is given over to an attack on the 


3 Paul F. Schurmann, Historia de la Fisica (Anales de 
la Universidad. Entrega no. 139, 1936. Republica Oriental 
del Uruguay, Montevideo), p. 160. 

4Sylvanus P. Thompson, Gilbert, of Colchester; an 
Elizabethan Magnetizer (London, 1891), p. 9. See also 
Charles E. Benham, William Gilbert of Colchester (Col- 
chester, 1902), and A. Herpin, William Gilbert et la Méde- 
cine au X VI® siécle (Paris, 1946). 

5 De Mundo (Amstelodami, 1651), a rare book, is avail- 
able at the libraries of Columbia University and the New 
York Academy of Medicine. Pages 270-276 are on the 
rainbow. 
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meteorclogy of Aristotle (384-322 B. C.); and 
is this not what one should expect from reading 
about Gilbert? “‘Whilst the fantastic philosophies 
of the schoolmen still prevailed, he calmly 
worked out the inductive method of reasoning 
from the known to the unknown, trying all his 
arguments by the touchstone of experiment.”’® 
Aristotle’s elaborate geometrical and physical 
theory of the primary and secondary rainbows, 
based upon the reflection of solar rays by the 
aggregate of drops in a dewy cloud,’ had been 
questioned many times before Gilbert’s day. In 
particular, Avicenna (980-1037) had modified 
the Peripatetic views by suggesting that there 
were two effective agents in the reflection: (1) a 
thin vapor between the sun and a dark cloud or 
hill, and (2) the surface of the cloud or hill. 
These agents he compared to the glass and mer- 
cury components in mirrors.* During the thir- 
teenth century the theory had undergone more 
drastic change through the introduction, espe- 
cially by Robert Grosseteste (1175-1253) and 
Witelo (fl. 1270), of the phenomenon of refrac- 
tion; and early in the following century a quali- 
tatively correct explanation of the rainbow (in- 
cluding both the primary and secondary bows) 
through the refraction and reflection of rays of 
sunlight as they strike individual drops of rain 
had been presented in detail by Theodoric of 
Freiburg (and also, independently, by Arabian 
contemporaries). However, the theory of Theo- 
doric had disappeared by the middle of the six- 
teenth century, and Jerome Cardan (1501-1576) 
was typical of his time in reverting, in the De 
Subtilitate of 4550, to earlier Peripatetic views, 
thinking only in terms of reflection. Franciscus 
Maurolycus (1494-1575) likewise overlooked re- 
fraction, explaining the rainbow geometrically 
through manifold reflections within and without 
the individual raindrops. Johann Fleischer (1539- 
1593) in 1571 published a theory, which he 
modestly regarded as a reconstruction of the 
views of Aristotle and Witelo, according to 
which the sun’s rays were refracted by some of 


* Thompson, reference 4, p. 10. 

7 For references on the theory of Aristotle and some of 
his successors, see my paper on ‘“‘Kepler’s explanation of 
the rainbow,” Am. J. Phys. 18, 360-366 (1950). 

8M. Horten, ‘‘Avicennas Lehre vom Regenbogen nach 
seinem Werk al Schifa,” Meteorologische Zeitschrift, 30, 
533-544 (1913). 
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the drops and reflected by others. The indi- 
vidual-drop theories of Maurolycus and Fleischer 
seem not to have been widely known, and the 
older cloud-as-a-whole views prevailed at the 
time when Gilbert took up the problem. 

The last three of the five books of De Mundo 
are devoted to meteorology,® and it is midway 
in this account that Gilbert presents his views 
on the rainbow, beginning with chapter x of 
Book IV. The illumination of the rainbow, he 
wrote, is generated in a vapor and is checked or 
brought back by a dark object, just as the image 
is not reflected by a mirror unless tin or a similar 
substance is placed behind the glass. Thus the 
rainbow, an image of the sun, is not transmitted 
to the observer unless there is an object nearby, 
such as a cloud, a cliff, a mountain, or a building. 
The primary rainbow is formed by the vapor 
lying nearby and surrounding us, being reflected 
by a dense cloud or a high cliff. Nevertheless, the 
bow is not generated within the cloud, just as, 
in the case of glass mirrors, the light which ap- 
pears to be reflected by the glass, and is in 
reality reflected by the tin or pitch placed behind 
the glass, is not in the tin or pitch. The second- 
ary rainbow is formed by vapor near to, or 
within, the cloud by whose darkness it is re- 
flected; and hence it is appropriate that here 
the brighter colors should arch around nearer to 
the primary bow, whereas in the case of the 
primary rainbow, which is reflected from a dis- 
tant cloud, the order of colors is inverted, with 
the brighter outermost. Where the illumination 
is smaller, due to a greater refraction of the rays, 
there the color is more intense. Thus there is an 
analogy between the generation of the luminous 
bow in the vapor and its colors. However, the 
light will not be reflected to the eye unless there 
is some dense object not too far away for it to 
appear, and unless the proportion is just right. 

The above free translation of chapter x con- 
tains the essentials of Gilbert’s theory. His ideas 
—vague, speculative, and tediously repeated— 
resemble closely those of Avicenna; but this 
does not necessarily mean that they were derived 
directly from Muslim sources, or indirectly from 
somewhat similar views expressed by Albertus 

9 The first two are on “‘physiology,”’ i.e., the constitution 


of the earth and the heavens. Gilbert’s ‘‘physiology” and 
“meteorology” both are directed contra Aristotelem. 
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Fic. 1. Illustration showing the vagueness of Gilbert's 
idea of the relationship between the dewy vapor and the 
dark cloud in the formation of the rainbow. From De 
Mundo, p. 272. 


Magnus (t 1280). Gilbert was predisposed to 
immoderate attacks upon medieval thinkers, 
both Latin and Arabic; and it is possible that 
his explanation was either independently arrived 
at or else was suggested by other scholars of the 
sixteenth century. That he was not alone in his 
views is seen in the fact that Johann Kepler 
(1571-1630), a great admirer of Gilbert, es- 
poused a somewhat similar theory in 1604,!° 
although he renounced it in the following year. 
More than half a century later Cureau de La 
Chambre (1594-1669), physician to Louis XIII 
and Louis XIV and one of the first members of 
the Académie des Sciences, published a long 
treatise on the rainbow which was based upon 
analogous medieval notions." Kepler and La 
Chambre, however, distinctly brought in re- 
fraction, whereas Gilbert seems to be thinking 
solely in terms of reflection. Gilbert’s account 
does indeed at one point include the word re- 
fraction, but this is inconclusive inasmuch as 
the word often had been used, since the time of 
Aristotle, as synonymous with reflection. All in 
all, his explanation strikes one as an artless 
version of medieval philosophizing, and this im- 
pression is not dispelled by the following 
chapters. 

In chapter x1 Gilbert undertakes to refute 
an ancient error which had been maintained 
more or less frequently throughout the sixteenth 
century and which Gilbert mistakenly attributed 


10 In his treatise Ad Vitellionem Paralipomena (Franco- 
furti, 1604). 

11 Nouvelles Observations et Coniectures sur I’Iris (Paris, 
1650, 1662), second edition. 
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to Aristotle—the idea that the secondary rain- 
bow is a mirror-image of the primary. In par- 
ticular, this error had been espoused not only 
by Cardan, but also by his older mathematical 
contemporary, Josse Clichtowe ( fT 1543), against 
whom Gilbert directed his attack. Inasmuch as 
the position and curvature of the secondary are 
the same as for the primary, he pointed out, 
then the colors, too, ought to be similarly or- 
dered, were the one a mirror-image of the other. 
To explain the inverted order of the colors in 
terms of images in a mirror, it would be neces- 
sary to have the tips of the bows bending in 
opposite directions, something which is never 
seen to happen in the atmosphere. Then, too, 
the secondary would have to be brightest when 
the primary is at its brightest, whereas, Gilbert 
thought, the former is more likely to appear 
when the latter begins to disappear. Gilbert here 
labors to demolish a view which had been simi- 
larly refuted in ancient, medieval, and early 
modern times. 

In chapter twelve he considers more closely 
the shape and colors of the rainbow. The bow is 
a sort of reflection, not of the form of the sun, 
but of a representation of the illuminated vapor 
carried to the eye from the air which is opposite. 
If the air is dense and sufficiently smooth and 
placed directly opposite, it will reflect the un- 
distorted form of the sun (as happens in the case 
of glass mirrors). However, the rays do not cling 
together but are diffused in the vaporous air, 
the illumination of which has been increased, 
and concentrated by the force of the light from 
the spherical drops. Then an image, with a 
variety of colors produced when the luminous 


Fic. 2. Diagram showing that the colors (red, yellow, 
green, and purple) in the secondary or outer rainbow are 
ordered inversely to those in the inner “prima iris.”” From 
De Mundo, p. 275. 
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rays ar? weakened in traversing the moist vapor, 
is reflected if something dense, like a cloud or 
mountain is interposed, just as images are not 
reflected by glass unless it is backed by tin or 
asphaltum or some such substance. Moreover, 
it is not truly a reproduction, but a figure arising 
from a mixture of light and vapor. It is a repro- 
duction of light rather than of form and color, 
just as in the calcination of lead by a slow fire 
there is at first a tincturing with a green color, 
then with yellow, and finally, after long and in- 
tense heating, with red as it is converted into 
oxide of lead. Thus in the rainbow the stronger 
light is of a red color, the weaker is yellow, and 
the still feebler is green. 

The distances between the objects and the 
eye, or the great spaces between bodies, do not 
affect the proportions or shape of the bow. If 
the eye is at A (see Fig. 1), and if the interposing 
mountain (against the side of which the rainbow 
is perceived) is at B, a distance of a stadium 
[about 606 English feet] or less, and if a cloud 
is at F, two thousand paces away, the rainbow 
GDC nevertheless will be uniform, and not 
broken or distorted—that is, if between D and 
E there is no cloud seen in the sky but only 
pure, thin, and empty air. 

Chapter xm is devoted to another assault 
upon the already prostrate form of the Clich- 
tovean hypothesis, and the old argument is re- 
peated in superfluous detail. When one views a 
reflection in water, the portion which is higher 
above the water appears lower in the reflection, 
and vice versa. Consequently, if the point of a 
knife is directzd toward the water with the handle 
toward the sky, in the image in the water the 
pointed end of the iron portion will be upper- 
most and the handle will be lowest. And so let 
it be assumed that a log of wood is divided into 
three parts: a, the highest, red; b, the middle, 
green; and c, the lowest, purple; and let it be 
placed above clear water. Then in the image 
appearing in the water, the highest part, c, will 
be purple; the middle part, 5, will be green; and 
the lowest part, a, will be red. This lofty demon- 
stration, says Gilbert, is inappropriate, inept, 
and wholly worthy of Aristotelian subtlety. 
“But what do you say, Clichtovaeus. In the 
general ordering of things should not the entire 
proportion and form be inverted? See how, your 


Fic. 3. Diagram illustrating the essential elements of 
Gilbert’s theory: rays from the sun are transmitted 
through the air, then are tinted upon traversing the 
watery vapor, and finally are reflected to the eye by a dark 
cloud or dense object. From De Mundo, p. 276. 


own demonstration refutes you.’’ And Gilbert 
again points out that under this hypothesis the 
primary and secondary bows should be convex 
toward each other (not concentric), with the 
horns of the secondary pointing upward. Hence 
the secondary is not an image of the primary, 
either in air or in the clouds. 

The infelicity of Gilbert’s account of the rain- 
bow must be admitted by his most sedulous ad- 
mirers; but the purpose of this paper is not 
primarily to evoke such an admission or even to 
call attention to one more of the many erroneous 
theories to be found during the sixteenth cen- 
tury. The object is rather to point out the extent 
to which Gilbert here departs from the method- 
ology and scientific habit of thought which are 
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presumed to be so characteristic of his work— 
from the “new style of philosophizing’’ of which 
he boasts in the preface of De Magnete. The 
above passages, if placed beside comparable 
sections of Grosseteste or Witelo, will be found 
to be no less medieval in spirit and method ; and 
a comparison with the fourteenth-century work 
of Theodoric will be positively devastating to 
Gilbert’s vaunted reputation as “the great 
father of experimental philosophy.” It is not 
without reason that the De Mundo has been 
described as “a nonexperimental ‘philosophia 
contra Aristotelem,’ a typically Renaissance 
treatise which attacked Aristotelian doctrines in 
perceptibly peripatetic style.” A just recogni- 
tion of this fact makes it less difficult for one to 
appreciate the strictures which Francis Bacon 
(1560-1626) passed on the writings of his brilli- 
ant countryman. Nor was such criticism an iso- 
lated case. Priestley candidly reported of the 
De Mundo that it “abounds with scholastic re- 
finement, and contains the most frivolous dis- 
tinctions and definitions. . . . None of Gilbert’s 
opinions on any optical subjects are worth re- 
citing.’”"* How different is the ring of such 
passages from the familiar later judgments of 
De Magnete: a work ‘‘as remarkable for its sci- 
entific spirit as for its content’’;!® ‘‘a sober trea- 
tise based on simple facts in place of the wild and 
speculative treatises which had hitherto passed 
as philosophic.’’!*® Were the De Mundo a juvenal- 
zum, it would not be fair to compare it with the 
product of the author’s mature thought; but the 
evidence indicates that it was prepared after 
1651,!7 ie., probably later than De Magnete.'* 
One might argue further that the posthumous 
publication of De Mundo suggests that the author 
did not take this speculative work seriously ; but 


12 See article, “William Gilbert of Colchester,’”’ Elec. 
World and Eng. 42, 480 (1903); also one by Brother 
Potamian, “Gilbert of Colchester,’’ Pop. Sci. Monthly 59, 
336-350 (1901). 

13 Marie Boas, “Bacon and Gilbert,” J. History Ideas 
12, 466-467 (1951). 

14 Joseph Priestley, The History and Present State of Dis- 
a Relating to Vision, Light, and Colours (London, 
1772 


18 F, Sherwood Taylor, The March of Mind (New York, 
1939), p. 147. 

16 Thompson, reference 4, p. 36. 

17 Edgar Zilsel, ‘“‘The origins of William Gilbert’s sci- 
entific method,” J. History Ideas 2, 1-32 (1941), especially 
1-3 


18 Benjamin, reference 2, pp. 316-320. 
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here, too, evidence points to the conclusion that 
“Gilbert clearly attached great importance to 
his cosmological and physical speculations.’’® It 
appears that Gilbert undertook his work on 
magnetism largely in order to apply it to his 
theories of the heavens.2° The fact that De 
Mundo appeared after the author’s death sug- 
gests also the possibility that the material was 
not in final form; but this factor (which in any 
case would not materially affect the thesis of this 
paper) is mitigated by the circumstance that the 
editor (the brother of the author) collated two 
manuscript copies of the work. How, then, does 
one account for the fact that Gilbert’s study of 
the magnet differs toto caelo from his approach 
to the problem of the rainbow? One is tempted to 
hold that the former lent itself more easily to 
experimental investigation, whereas the latter 
called for quantitative analysis, a field in which 
Gilbert was relatively less competent. Neverthe- 
less, the work of the fourteenth century would 
seem to belie such an explanation. Just as the 
idea of the earth as a huge magnet led Gilbert 
to experiments with a “‘terrella,” so the exactly 
analogous idea of a raindrop as a minute sphere 
led Theodoric (as well as the Arabs, and, later, 
Kepler, De Dominis, and Descartes) to experi- 
ments with a ‘‘globosum’’—a large glass globe 
filled with water. Was the idea of a magnified 
raindrop more elusive than that of a miniature 
earth? And Theodoric seems to have been as 
weak on the quantitative side as was the great 
“electrician’’; his success in the theory of the 
rainbow resulted from the very same factors 
that served Gilbert so effectively—from the 
care with which he made qualitative observa- 
tions. There seems to be little room for doubt 
that what Gilbert did for magnetism and elec- 
tricity he might have accomplished also for the 
rainbow, had he but applied to the latter phe- 
nomenon the same experimental philosophy 
which he had exploited so successfully in con- 
nection with the former. By a strange coinci- 
dence it appears that one of Gilbert’s ac- 
quaintances, the mathematician Thomas Harriot 
(1560-1621), was perhaps even then experiment- 
ing with globes of water in an effort to decipher 


19 Boas, reference 13, p. 467. 
20 Benjamin, reference 2, pp. 321, 327. 
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the rainbow problem. Correspondence between 
Harriot and Kepler in 1606, only three years 
after Gilbert’s death, indicates that these men 
had reached, independently, the idea of the rain- 
bow as resulting from reflections and refractions 
by spherical raindrops. Had Gilbert lived only a 
few years longer, his theory of the rainbow in- 
evitably would have been completely different. 
The obvious inference is that, in some respects 
at least, Gilbert definitely was not “ahead of 
his time,’’ as is frequently maintained. It is not 
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given to any one man to make all the discoveries 
or to transcend all error. If Gilbert’s dexterous 
experiments on electricity strike sparks of mo- 
dernity, let his devious arguments on the rain- 
bow be a reminder that his mental habits were 
those of an earlier day. Taken together, De 
Magnete and De Mundo show Gilbert to have 
been a typically Janus-faced physician of the 
Renaissance, modern in his advocacy of labora- 
tory precision, yet medieval still in his predilec- 
tion for fuzzy speculation. 


Selection Rules for Vibrational Spectra of Linear Molecules* 
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An elementary explanation is given of the group theory method using characters to deter- 
mine the selection rules for the fundamental vibrations of a linear polyatomic molecule. Then 
the method is extended to indicate how the selection rules for the overtones and summation 
and difference tones are obtained. Finally, the procedure for determining the number of funda- 
mental frequencies belonging to a particular vibration type is given. The groups Ca» and Daa 
are used as examples in the treatment since they are the only ones involved in the considera- 


tion of linear molecules. 


N a previous paper,! an explanation was given 

of the application of group-theory methods to 
the determination of selection rules for vibra- 
tional spectra of nonlinear polyatomic molecules. 
Since some modification of the concepts and for- 
mulas given in this earlier paper becomes neces- 
sary in the ccnsideration of linear molecules, it 
seems worth while to give a more elementary 
account of the methods used than can be found 
in the more advanced references.” This treat- 
ment also seems desirable because it is possible 
to understand the techniques used without too 
great a knowledge of the theory of continuous 
groups. 


* Supported in part by Office of Ordnance Research, 
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CONTINUOUS GROUPS 


In determining the selection rules for linear 
molecules one must consider continuous groups. 
This results from the fact that there are an in- 
finite number of symmetry operations which will 
give an equivalent configuration of the molecule. 
For example, if the axis of rotation is chosen as 
the line on which the nuclei lie, there are an 
infinite number of rotations about this axis which 
are symmetry operations and which furnish an 
equivalent configuration of the molecule. Since 
the determination of the selection rules is based 
on the use of a reduction formula which for non- 
linear molecules involves a summation over all 
the elements of a finite group, it is not possible 
to make these calculations in the same way for 
linear molecules because one would have a sum 
to be taken over an infinite number of elements. 

Before proceeding to a discussion of the neces- 
sary modifications of the reduction formula, it is 
convenient to describe the two continuous groups 
Cs, and D., which are necessary for describing 
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the symmetry properties of linear molecules. 
The symmetry operations for the continuous 
groups C,, and D,, are given in Tables I and 
II, respectively. 

In Table I, the covering operations of the 
group C., have been divided into two sets of 
symmetry elements. The first, C(y) corresponds 
to rotations by an angle ¢ about the axis of the 
molecule (0< g¢<27). Similarly o,-C(¢), refers 
to the rotation C(¢g) followed by a reflection oc, 
in a plane passing through the axis of the mole- 
cule. The symmetry operations thus described 
include all the covering operations of the C.» 
group; the identity operation corresponds to 
C(0). In this way, by means of a continuously 
changing parameter gy, one obtains all the ele- 
ments of the continuous group. 

In Table II, a similar division is made for the 
elements of the D., group. In this case, C(¢) 
and o,-C(yg) have the same meaning as before, 
while C.-C(¢) represents the rotation C(¢) fol- 
lowed by a rotation of w-radians about an axis 
through the center of the molecule perpendicular 
to the axis. S(y) describes a rotation C(¢) fol- 
lowed by a reflection in the plane perpendicular 
to the plane of the molecule and passing through 
the center of the molecule. One notices that o;, 
(reflection in the plane through the center of the 
molecule perpendicular to the axis) is simply 
S(0) and the inversion 7 is identical to S(7). 

The vibration types or species (irreducible 
representations) are listed in the first column of 
each table. An A vibration type always refers 
to a nondegenerate vibration, while an E vibra- 
tion type refers to a doubly degenerate vibration. 
In the case of the nondegenerate vibrations, the 
subscripts J and 2 refer to vibrations symmetric 
and antisymmetric, respectively, with respect to 
reflections o,. There are an infinite number of EZ 
types which are designated by the use of the sub- 
script k where k& is a positive integer so that 


TaBLeE I. Characters for the C..» group. 


Cov C(¢) 
A; 1 
As 1 
Ex 2 coskey 
xu(R) 1+2 cose 
Xa(R) 2+2 cosg+2 cos2¢ 
N 


u 
=(R) 


av 'C(¢) 


2(N—2)cose+N-1 
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k=1, 2, ---. The subscripts g and u always refer 
to vibrations symmetric and antisymmetric with 
respect to the inversion. 

A quantity called the character (obtained by 
use of group theory) is necessary for the deter- 
mination of the selection rules and the number of 
fundamental vibrations of each type. For a 
given vibration type, there is a separate char- 
acter for each set of symmetry operations; thus, 
for a molecule of symmetry C.», there are two 
characters for type A1, two for type Az, and two 
for each E vibration type, since there are two 
sets of operations for this group. These char- 
acters for the C,, and D., groups are listed in 
the upper part of Tables I and II, respectively. 
(As will be seen later, the character is the sum 
of the principal diagonal elements of a matrix.) 

The fundamental problem in determining the 
selection rules and the number of vibrations be- 
longing to a given vibration type is to express 
the character of certain physical quantities as a 
sum of the characters of the vibration types 
(irreducible representations) of the group. In the 
case of finite groups, this is accomplished by 
means of the reduction formula 


ai=(1/h)2 x(R)x*(R), 


where a’ is the number of times that x/(R) (the 
character of the jth vibration type for any 
operation R) is contained in x(R) (the character 
of the physical quantity under consideration), 
h is the number of elements of the group, x/(R) 
is the complex conjugate of x/(R), and the sum 
is taken over all the elements of the group. In 
the earlier paper on nonlinear molecules,! this 
formula was applied to a number of examples. 

For continuous groups, it is necessary to re- 
place this summation with an integral and the 
reduction formula becomes‘ 


ai= f x(R)x(R)aR / fx @ 


where the integration is taken over the range of 
the continuously changing parameter R (in this 
case ¢g), to cover all the operations of the group, 
and the other symbols have the same meaning 
as before. 


4 See reference 2, page 110, Eq. (12). 
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INFRARED SPECTRA 


In order to determine which fundamental 
vibrations are allowed in the infrared spectrum 
it is necessary to obtain the character x4(R) of 
the electric dipole moment for the different 
operations of the group. Only those vibration 
types whose characters are contained in the 
character of the dipole moment will be allowed 
in the infrared.5 

Since the dipole moment M transforms as a 
vector under the different covering operations of 
the group, one has 


M;’ cosg sing 0 
M,’' }={ —sing cosg 0 M,}, (2) 
M;’ 0 0 1 M, 


for a proper rotation of angle ¢ about the z axis 
(the axis of the molecule). One recognizes the 
3X3 matrix as the ordinary transformation 
matrix for a rotation. The C(¢) operation of the 
Cx, and D.» groups is of this type. 

Since the character is the trace, or sum of the 
elements of the principal diagonal of the repre- 
sentation matrix, (in this case the transforma- 
tion matrix) one sees that 


xuLC(¢) ]=1+2 cos¢. 


For an improper rotation (a rotation followed 
by a reflection in the plane perpendicular to the 
axis of rotation) about the z axis, such as S(¢), 

M,' cosg sing 0 M,z 

M,' }={ —sing cos¢ 0 M,}, (3) 

M,’ 0 0 —1 M, 
whence 


Mz 


x4LS(¢)]= —1+2 cos¢. 


Now, it is necessary to obtain the character 
of the dipole moment for the operations o,-C(¢) 
and (,-C(g). One might assume that these 
characters could be obtained by multiplying 
the characters of the two representations in- 
volved, but such is not the case. Instead, one 
must determine the representation matrices for 
the two operations, then multiply the two 


9 


WP itincinen 


TABLE II. Characters for the D4 group. 


C(¢) C2°C(¢) S(¢) ov°C(¢) 


1 1 
1 —1 
1 1 
1 
2 coske 0 
ke 2 coske 0 
1 


1+2 cos¢ - 
2+2 — cos2¢ 2 


—1 
2(—1)* coske 
—2(—1)* coske 
—1+2 cose 
2 —2 cose +2 cos2¢ 


xm(R) 
xa(R) 


u mo mo 
=(R) 2(N —2)cosep+N—1 —mo+1 —mot+1+2mocose N—1 


matrices together, and take the sum of the prin- 
cipal diagonal elements of the product. 

For the operation ¢,, which can be considered 
as an improper rotation of angle zero about an 
axis perpendicular to the z axis (say the x axis), 
one has 

M,’ —-1 0 0 
M, j= . i. 4 
M,' 0 0 1 


M, 
M, }. 
M, 
Hence 
will be 


—1 cosy sing 0 
0 ‘| —sing cosy’, 0 
0 0 0 #1 
—cosg -—sing 0 
={ —sing cosy 0 }, 
0 0 1 


the o,-C(¢) transformation matrix 


whence 
xuLor a C(¢) J =1. 


Also, for the operation C2 (a proper rotation 
of x-radians about the x axis) one has 


M,! 1 0 0 
M,'}=[0 -1 0 
M/} \o oO -1 


MM, 
M, }, 
M, 
so that 


xmLC2-C(y)J=—-1. + 


The values of x(R) are given in Tables I and 
II for the operations of the continuous groups 
under consideration. 

Now, it is necessary to reduce xm(R) by using 
Eq. (1). For the C,, group, one has 


f * pal CAC et f xaos Clo) ixLee- Cle) Ke 


0 
2x Qn 
f aot f dg 
0 0 


5 J. E. Rosenthal and G. M. Murphy, Revs. Modern Phys. 8, 317 (1936), section 27. 
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Using 
2x 
f coskg cosk’ ody = Tix: 
0 


where 
1 for k=k’ 
by -| 


0 for R¥k’, 
one obtains 


2x 
ati=(2/4m) f (1+ cosy)dy=1, 
0 





2r 


J 


+ f xuLS(¢) XL S(e) e+ f ker Cin Chae 


2r 2r 2r 2x 
f ao+ [ do+ f ao+ f dg 
0 0 0 0 


and similarly 
a410 = q42u = gA20 = qEko = 0, 
a4uw=1, 
aPku = 5,1; 


xa (R) = X41u(R) +x#1.(R), 
which means that only the A;, and E,;, funda- 
mental vibrations are allowed in the vibrational 
infrared spectrum. 
RAMAN SPECTRA 

In order to determine which fundamental 
vibrations are allowed in the Raman spectrum, 
it is necessary to obtain the character x.(R) of 
the polarizability (a) for the different operations 


y sin’g 


4 sin2¢ 
0 
0 0 


0 
0 
0 1 
0 
0 


where the upper signs (in +cosg, etc.) refer to 
the proper rotations C(¢g) and the lower signs to 
the improper rotations S(¢). 

Hence, the characters for the polarizability are 


Xal C(y) ]=2+2 cose+2 cos2¢, 
and 
XalL.S(¢) ]=2—2 cose+2 cos2y¢. 


The transformation of a under the operation 
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a42=0, 
ak = §,3. 


Hence 
xm(R) = x41(R)+x#1(R), 


which means that only the A; and £, funda- 
mental vibrations are allowed in the vibrational 
infrared spectrum. 

For the D,,, group one has 


xulC(~) xiLC(~) de+ f xa Co Cl xT Cr Cle) Bo 


» (5) 


of the group. Only those vibration types whose 
characters are contained in the character of the 
polarizability will be allowed in the Raman 
spectrum. ® 

The six distinct components of the polariza- 
bility tensor transform under the operations of 
the group as follows: 

ae;;’ ™ Zz Chan U on Un, 

where all the subscripts extend over x, y, and z 
and the Uim are the elements of the transforma- 
tion matrices used in Eqs. (2) and (3) for proper 
and improper rotations, respectively. 

For a proper or improper rotation by an angle 
yg about the z axis, one has 


sin2¢ 0 
—sin2¢ 0 
0 0 
cos2g , 0 | ey |’ 
0 itcose +sing]| jaz, 

Fsing -+tcosg 


Azz 
yy 
Qzz 


Ayz 


and that for C2 is also the same. 
6 Reference 5, section 28. 
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Proceeding in the same way as for the dipole 
moment, one finds that 


XeLo»-C(y) ]=2 
and 


XaLC2-C(y) ]=2. 


In Tables I and II, the values of x.(R) are 
given for the operations of the continuous groups 
under consideration. 

Now, it is necessary to reduce xe(R) using 
Eq. (1). For the C.» group, one has a relation 
the same as Eq. (4) except that the x(R) are 
replaced by the x2(R). Hence one finds that 


a4i= a/anyf (242 cosg+2 cos2¢+2)dy=2, 
a42=0, ; 
aPk = y+ dx2; 
whence 
Xa(R) = 2x41(R)+x21(R)+xz2(R), 


which means that only the A,, £,, and E, funda- 
mental vibrations are allowed in the Raman 
spectrum. 

For the D,,, group, again, one has a relation 
similar to Eq. (5) so that 


a4w=2, 
q4lu = q4% = qEku = 0, 
aPko = §41+ d42; 

thus 


Xa(R) = 2x419(R) + x#19(R) + x#29(R), 


which means that only the Ay, Ey,, and Eos, 
fundamental vibrations are allowed in the 
Raman spectrum. 


COMBINATIONS 


To obtain the selection rules for combination 
frequencies v;-v;, it is necessary to form the 
direct product of the characters of the vibration 
types to which »; and »; belong. This is done by 
multiplying the characters of the vibration types 
for each operation of the group. Thus, for the 
group Cx, the direct product Ai XE; is obtained 
as follows 

C(¢) ao C(¢) 
xA1(R) 1 1 
XEx(R) 2 cosky 


(multiply). 


XA1 XEK(R) 2 cosky 0 


It is seen that 
XA1 XEK(R) = XEK(R), 


which means that if a frequency »; of type E;, 
combines with a frequency »; of type A;, the 
resulting combination frequency »;-+»; will be 
allowed wherever an EF, type frequency is al- 
lowed. For the present case, the combination 
frequency will be allowed in both the Raman and 
the infrared if R=1, only in the Raman if k=2 
and it will be forbidden for both the Raman and 
the infrared for all other values of k. 

Selection rules for other combination fre- 
quencies are obtained in the same way. In 
general, it is necessary to use the reduction for- 
mula, Eq. (1). This may be illustrated by the 
combination of a frequency »; of type FE, with 
another v; of type E,. One has 


C(¢) a»: C(¢) 
2 coskey 0 
2 cospe 


XEKMR) 


E,(& 
—_ (multiply). 


XEx XEp(R) 4 cosky cospey 0 


By means of the reduction formula one finds 


aA1= d,5, 
a42= d,5, 
a¥»=(0, unless n=k—p or n=k+); 


in which case, for k>p 


aEkt+rp=1, 

at-e=1, 
and for k=p 

aE =}, 


The selection rules for combinations for the 
Co» group are summarized in Table III. A similar 
table can be written for the D,, group taking 
into account the facts that two vibration types 
having only g or only u subscripts give vibration 
types with g subscripts; and that for combina- 
tions for which one vibration type has a g sub- 
script and the other a wu subscript, one gets 
vibration types with only wu subscripts. This 
table is not given here since it is similar to 
Table III, except that the g and u subscript rule 
must be taken into account. 


OVERTONES 


To obtain the selection rules for the overtones 
of nondegenerate vibrations, one proceeds in a 
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manner similar to that used for combination fre- 
quencies. Thus, the character of the (w—1)th 
overtone is given by 


xi"(R) = (x:(R)]". 


For example, for the 1st overtone of an Az» 
frequency »v; of the C,, group the character is 
obtained as follows 


C(¢) a» C(¢) 
XA2(R) 1 —1 
R 1 
— (multiply). 
x 472(R) 1 1 


It is seen that x42*(R) = x4:(R) which means 
that the 1st overtone of an A2 frequency will be 
allowed where an A, frequency is allowed, in this 
case in both the Raman and infrared. 

In a similar way, one sees that 


x41"(R) = x41(R), 
and 
x41(R) for n even 


4o"(R) = 
x42"(R) x42(R) for n odd. 


Consequently, for a molecule of symmetry 
Cav, all overtones of A, frequencies but only the 
A; overtones for which 7 is even are allowed in 
both the Raman and infrared. 

For the D., group, one has 


X419"(R) = x419(R), 


xX41,(R) for n even, 
X41u"(R) - 
X41u(R) for 2 odd, 
x41g(R) for n even, 
X42,(R) for 2 odd, 
x41g(R) for n even, 
X4A2u(R) for n odd. 
One sees that the symmetry types of the over- 
tones for the D,, group are the same as those for 


X429"(R) = 


x42u"(R) = | 


TABLE III. Symmetry types for binary combination 
frequencies of a linear molecule with symmetry C @». 





Combinations 


Symmetry type 


AiXAi1 A; 
AiXE; Ex 
A2XAz2 A, 
AoXEx Ex 

Ex XEp Exspt+Ex_p 
(k>p) 

ExrXEx A,+A2tEx 
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the C.. group, if one takes into account the g 
and u subscript rule given in the last section. 
The method of obtaining selection rules for the 
overtones of the degenerate vibrations is some- 
what more involved. One must use the formula’ 


xep"(R) = 3Lxep" "(R)xep(R)+xz,(R")], (7) 


where xz,(R") means the character for the opera- 
tion R performed successive times and xz,""1(R) 
is the character of the (n—2)th overtone for the 
operation R. For the group C,., the values of 
xep"LC(¢) ] are: 


for m even xzp"LC(¢) ]=1+2 cos2pe 
+2 cos4py+ ---+2 cosnpe, 
and (8) 
for n odd xz," C(¢) ]=2 cospe 
+2 cos3py+2 cosSpy+ ---+2 cosnpe. 


For example, by use of Eq. (7) 


xEC(¢) =2 cos?*o+cos2y=1+2 cos2¢, 
xz1°C(¢) =cosg+2 cosy cos2y¢+cos3¢ 
=2 cosg+2 cos3¢, 
xe1'C(¢) =2 cos*o+2 cosy cos3¢ 
+cos4g=1+2 cos2¢+2 cos4y, 


and by induction one obtains the general for- 
mulas given in Eqs. (8). 

For the operation o,-C(¢), one has xz,[Lo, 
-C(¢g) ]=0 from Table I, and from Eq. (7) one 
obtains 


xEp"Lo» ‘ C(¢) J = 3XEp{ Lo, P C(¢) ]"} ’ 


2 for n even; 
. mn 
xEp{Lo» C(¢) ] } 0 for n odd; 
since performing o,-C(¢g) an even number of 
times reproduces the original configuration so 
that it is equivalent to the identity element of 
the group. For ” odd, 


XEp{ [o, I C(¢) }"} = XEpLor ‘ C(¢) J =0. 


Applying now the reduction formula, Eq. (1), 
one obtains: for ” even 





a4i=1, 

a42=0, 

aPk = byop+ deapt - +++ din; 
or 


XEp"(R) =x41(R)+ X#2p(R) 


=i: Spee eal +x#4p(R)+ oat +x£np(R) ; 
7™L. Tisza, Z. Physik. 82, 48 (1933). 
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and for odd, 


g4=¢a4:=0, 
aFk = §py+ disp t - ++ +5kn, 


or 


xEp"(R) = xep(R) + xe3p(R) 
+ xE5p(R)+ - +> +x#np(R). 


For the D,, group the results are the same, but 
one must again take into account the g and u 
subscript rule. Hence, for even 


XEpo"(R) = XEpu" (R)=xa io(R)+ XE2p9 (R) 
+ x#4p9(R)+ ase + xEnpo(R) ’ 


while for ” odd, 


XEpg” (R) _ XEpo(R) + x£3p9 (R) ooo XEnpo(R), 
and 


XEpu"(R) = XEpu(R) + x#3pu(R)+ ee + X£npu(R). 


NUMBER OF FUNDAMENTALS OF EACH TYPE 


In order to find the number of fundamental 
vibrations of each type, one needs a quantity 
=(R), the character of the complete 3N—5-di- 
mensional representation, referred to a particular 
molecule, corresponding to the 3N—S5 degrees 
of vibrational freedom of a linear molecule, 
where JN is the number of atoms.$ 

Since one is interested first in obtaining the 
character of the 3N-dimensional representation, 
=:n(R), corresponding to all possible motions of 
the molecule, only the diagonal elements of the 
transformation matrix need to be obtained. The 
only diagonal elements which will be different 
from zero for an operation R of the group are 
the ones corresponding to nuclei whose equi- 
librium positions in space are not changed by 
the operation R. Thus, for C(¢g) the character of 
the 3N representation will be u(1+2cosg), 
where u is the number of atoms unchanged by 
the operation C(g). Since there are three co- 
ordinates associated with each nucleus which is 
unchanged, the 1+2cosg factor results from a 
transformation matrix like that given in Eq. (2). 
In a similar way, the character of the 3N repre- 
sentation for S(¢) is found to be u(—1+2 cos¢). 


8 Reference 5, section 23. 


TABLE IV. Characters of the complete representation 
for a linear molecule. 


C(¢) 


Esn(R) u(1+2 cose) u —4u 
Zu(R) 1+2 cose 1 —1 
Erot(R) 2 cose 0 0 


=(R) 


av *C(¢) C2°C(¢) S(¢) 


u(—1+2 cos¢) 
—1+2 cose 
—2 cose 


u—1+2(u—2)cosp u—1 —u+1 2ucose—u+l1 


For o,-C(¢g) and C2-C(¢g), one uses the corre- 
sponding matrices given in the section on the 
dipole moment. The values of Z:y(R) are sum- 
marized in Table IV. 

For the translation of the molecule as a whole, 
the components of the position vector of the 
center of mass must be considered, which trans- 
form under any operation R like the components 
of any other vector (e.g., the dipole moment). 
The values of the character =;,(R) are also given 
in Table IV. 

For the rotation of the molecule as a whole, 
one considers the angular momentum compo- 
nents. However, in the case of linear molecules, 
only two components are involved, giving a two- 
dimensional representation whose characters are 
given in Table IV. They are obtained by con- 
sidering the transformation of the x and y com- 
ponents of the angular momentum since the 2 
component is zero. 

Finally =(R) is obtained from the relation 


=(R) = Zsw(R) — Zer(R) — Zrot(R), 


and the values are given in the last row of 
Table IV. 

The values of =(R) for linear molecules of 
Cz» and D., groups are given in Tables | and 
II. For the C., group, u is equal to N the total 
number of atoms in the molecule. However, for 
the D., group u is not the same for all the opera- 
tions; N is used to designate the total number 
of atoms as before and mp means the number of 
nuclei on all the elements of symmetry which in 
this case can take the value of 1 or 0. 

Now, it is necessary to reduce =(R) using the 
reduction formula, Eq. (1). For the C.. group 
one obtains 

a4i=N—1, 
a“2=0, 


aFk=(N—2) 6x1; 
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and similarly for the D., group one has 


a4w=}(N—m)), 

a4u=3(N+m—2), 

a4 = qA2u = 0, 

aEto=3(N—mo—2) 6x1, 
and 


quku= 4(N+ mo — 2) bx1. 


One sees that the total number of fundamental 
frequencies is 


N—1+2(N—2)=3N—5S for Cav, 
and 


3(N—mo)+3(N-+ mo —2)+2[3(N —m— 2) 
+3(N+m—2)]=3N—5 
for Doh. 
When these formulas are applied to specific 
examples one has the following results: 


For the OCS molecule which has the symmetry 
Cs» and for which N=3, 


a41i=?, 

a“2=0, 

a®i=1, 
and 


akk=0 for k¥1. 


Hence there will be two fundamental frequencies 
of the A; type and a doubly degenerate fre- 
quency of the £, type for the OCS molecule. 

For the CO: molecule which has the symmetry 
D., and for which N=3 and m=1, 


a40=1, 
a4uw=1, 
aA% = qAu=q¥w=0, 
qFu={ 


aFt=qFiu=() for kR¥1. 


Consequently there will be one fundamental fre- 
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quency of each one of the types A1,, Aiu, and 
Ey, for the CO2 molecule. 

For the diacetylene molecule (H-—C=C—C 
=C—H) which has the symmetry D,., and for 
which N=6 and m=0, 


a4w=3, g4u=?2, 
a4%7=gqAu=0) 
a¥=2, 
aFu=?, 
a®tw=q¥i=() for 


’ 


kA1. 


Hence there will be 3A1,, 2Aiu, 2E1,, and 2Ey, 
fundamental frequencies for diacetylene. 

One notices that the number of frequencies be- 
longing to a certain type depends on the number 
of atoms in the molecule, but for the selection 
rules this is not’ the case since the selection rules 
were determined by the symmetry of the mole- 
cule alone. In conclusion, it might be well to give 
some molecules as examples of those having the 
symmetry of the two groups, C., and Daa, since 
these groups are the only ones involved in the 
consideration of linear molecules. 

For the group C,, the following molecules 
have this symmetry: HCN, H—C=C-—D, N,O, 
OCS, CICN, (any heteronuclear diatomic mole- 
cule, HCl, HI, etc.). 

For the group D., these molecules have this 
symmetry: H—C=C—H, COs, CS: H—C 
=C—C=C-H, HgCh, C2Ne, (any homonu- 
clear diatomic molecule, He, Oz, Ne, etc.). 
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But the questioning of Socrates on the merits of these speculators sometimes took another form. 
The student of human learning expects, he said, to make something of his studies for the benefit 
of himself or others, as he likes. Do these explorers into the divine operations hope that when they 
have discovered by what forces the various phenomena occur, they will create winds and waters at 
will and fruitful seasons? Will they manipulate these and the like to suit their needs? Or has 
no such notion perhaps ever entered their heads, and will they be content simply to know how such 
things come into existence?—-XENOPHON: Memorabilia, Book 1, Chap. 1, (Translation of H. G. 


Dakyns). 
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A review of the present status of liquid scintillators is given. The theory explaining the 
properties of the counters is discussed. Preparation of counters is described and counter charac- 
teristics such as emission spectra, response to various nuclear radiations, counting efficiencies, 
and emission intensities from various scintillating solutions are given. The applications of 
liquid scintillation counters in various fields of nuclear physics are discussed. 


HE detection of radiation by means of its 

ability to produce luminescence in some 
materials dates back to the beginning of the 
century. In 1903, Crookes,! and also Elster and 
Geitel,?, observed that the light produced by 
alpha-particles on zinc sulfide consisted of indi- 
vidual flashes or scintillations which could be 
seen with a microscope. The importance of a 
detector of this type in the early development 
of nuclear physics may be judged from its use 
in such experiments as those of Geiger and 
Marsden* in 1913 on the scattering of alpha- 
particles, which verified the Rutherford model of 
the atom, as well as in the work of Cockroft and 
Walton‘ in 1932 on the first nuclear transmuta- 
tations produced by artificially accelerated 
particles. 

With the development of electronic circuits in 
conjunction with gas counters for the detection 
and counting of particles, the scintillation tech- 
nique, which depended on laborious visual ob- 
servation, was largely discarded. However, the 
advent of the multiplier phototube and the 
continued development of electronic techniques 
and scintillation materials has brought the scin- 
tillation method back into widespread use in the 
past several years, so much so in fact, that scin- 
tillation counters are tending to supersede gas 
counters in a number of applications. 

The two main properties with respect to which 
scintillation counters can surpass gas counters 


* Research carried out under the auspices of the AEC. 

t Now at Hudson Laboratories, Columbia University, 
Dobbs Ferry, New York. 

1W. Crookes, Proc. Roy. Soc. (London) 71, 405 (1903). 

2 J. E. Elster and H. Geitel, Physik A. 4, 439 (1903). 

3H. Geiger and E. Marsden, Phil. Mag. 25, 605 (1913). 


4J. P. Cockroft and E. T. S. Walton, Proc. R 


oy. Soc. 
(London) 137, 229 (1932). ‘ 


are efficiency and speed of response. The first 
property results from the development of large 
clear pieces of crystalline materials, and later of 
liquids, which give the scintillation counter a 
large stopping power. The duration of the scin- 
tillation, which, as will appear later, largely 
determines the rise time of the pulse, ranges 
from microseconds in some of the inorganic 
crystals to millimicroseconds in some of the 
organic crystals and liquids. 

The growing of large single crystals is a time- 
consuming process requiring methods of chemical 
purification and special ovens beyond the means 
and interests of many laboratories. Although 
good crystals are now available commercially, 
some applications require still larger sizes than 
can be grown. It was recognized that a trans- 
parent liquid scintillator would provide a simple 
way of avoiding the difficulties associated with 
the preparation and handling of large flawless 
crystals while making possible a much larger 
counting volume. A successful search for satis- 
factory solutions was reported by workers in 
Italy and the United States.5~’ It was found that 
small amounts of previously known organic 
scintillators dissolved in solvents, which by 
themselves scintillated only very weakly, gave 
good results. The scintillation properties of the 
solution resemble those of the dissolved organic 
material in regards to the duration of the light 
flash and its spectrum.®-" 


5 Ageno, Chiozzotto, and Querzoli, Atti accad. naz. 
Lincei. classe sci. fis. mat. e nat. 6, 626 (1949). 
ase Harrison, and Salvini, Phys. Rev. 78, 488 

7H. Kallmann, Phys. Rev. 78, 621 (1950). 

8 R. F. Post, Phys. Rev. 79, 735 (1950). 

A. Lundby, Phys. Rev. 80, 477 (1950). 

10F, B. Harrison and G. T. Reynolds, Phys. Rev. 79, 
732 (1950). 

" H. Kallmann and M. Furst, Phys. Rev. 81, 853 (1951). 
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TABLE I.*> Intensities of fluorescent light 
from pure liquids. 








Relative 
Solvent Solutes intensity 
Toluene 0.010-0.012 
Xylene (commercial) 0.012 
Cyclohexane 0.006 
Benzene 0.009-0.012 
Anthracene (0.01g) 0.002 
Terphenyl (0.03g) 0.005 
Xylene 0.013 
Xylene Carotene (0.1g/l) 0.001 
Xylene Fluorenone (0.5¢/1) 0.004 


® See reference 16. 
> The relative intensity is with respect to an anthracene crystal. All 
measurements were made with an 1P28 photomultiplier. 


The properties of crystal scintillation counters 
have been described in a number of articles.!?-13-34 
In this paper, we shall deal primarily with liquid 
scintillation counters and their applications. 


THEORY 

It has been known for quite a while that 
energy can be transported in solids through 
many atomic distances.'*:"® This phenomenon has 
been explained in terms of the motion of the free 
electrons and the migration of excitation energy 
from one molecule to another. Experimentally it 
has been shown" that the light output of some 
fluorescent crystals can be greatly increased by 
the addition of small amounts of impurity. If 
crystals of this type are excited by x-rays, y-rays, 
protons, or other ionizing radiation, they emit 
light which has the spectral characteristics of 
the impurity but has an intensity which is a 
function of the ionization in the crystal itself.'* 
These phenomena, which are explained by energy 
transport, suggested at once the possibility of 
similar processes in liquids which might be very 
useful in the development of radiation detectors. 
A number of investigators’-? found indeed that 
such a liquid scintillation counter was quite 
feasible. A brief review of the basic ideas of the 
theory of the fluorescence of solutions as devel- 
oped by Kallmann"'® is given here in order to 
make it possible to understand some of the 
properties of such counters. 


( 8 5). H. Jordan and P. R. Bell, Nucleonics 5, No. 4, 30 
1949). 

18 R. Hofstadter, Phys. Rev. 75, 796 (1949). 

13a G, T. Reynolds (to be published). 

144 J. Frank and E. Teller, J. Chem. Phys. 6, 861 (1938). 
( 1 E, J. Bowen et al., Proc. Phys. Soc. (London) A62, 26 
1946). 
16H, Kallmann and M. Furst, Phys. Rev. 79, 857 (1950). 
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Pure liquids, when excited, show very small 
emission intensities of fluorescent light®:®'%"” 
(see Table I). However, if solutions are used, 
then some of these solutions will exhibit rather 
strong fluorescence which in all cases has the 
spectral characteristics of the solute molecules 
(see Table II). Experiments'® show that the 
intensity of this fluorescence cannot possibly 
come from the primary excitation of the solute 
molecules of the solution but has to come from 
the excited solvent molecules. This fact suggests 
at once energy transfer from the excited solvent 
molecules to the solute molecules. 

In order to attempt to understand these 
experimental facts one has to consider the 
various aspects of the whole phenomenon: excita- 
tion of the solvent molecules, transport of the 
energy to the solute molecules, quenching of the 
energy of the excited solvent molecules in the 
solvent, retention of the energy by the solute 
molecules, quenching of the solute molecules, and 
finally the light emission from the solute 
molecules. 

The excitation energy of a molecule can be 
quenched two different ways. The molecule can 
interact with a similar unexcited molecule which 
then dissipates the energy by means of non- 
radiative processes and this we will call homo- 
quenching (sometimes called self-quenching in 
the literature). On the other hand, the molecule 
can also interact with dissimilar molecules and 
thus lose the energy again nonradiatively; this 
we will call hetero-quenching (sometimes called 
internal quenching in the literature). As homo- 
quenching can only occur when two similar 
molecules collide, it is, when it occurs between 
solute molecules, a function of the concentration 
of these molecules in the solution; however, 
hetero-quenching is independent of the solute 
concentration. Even though some quenching 
mechanisms have been proposed,!*:!® it must be 
recognized that the present status of quenching 
theory is not very satisfactory. This statement 
applies especially to hetero-quenching. 

When radiation passes through a scintillation 

17 Ageno, Chiozzotto, and Salvini, Phys. Rev. 79, 720 
i C Kortum and B. Finckh, Z. physik. Chem. B52, 263 


42). 
19G. K. Rollefson and H. Boaz, J. Phys. Chem. 52, 518 
(1948). 











solutior. it excites mainly the solvent molecules 
(on account of their greater concentration). As 
the pure solvent is an extremely weak scintillator, 
the excitation energy of the solvent molecules 
has to be transported to the solute molecule in 
order to obtain large emission intensities from the 
solution. This transport has to be fast enough to 
make homo-quenching of the energy in the sol- 
vent very small. Kallmann!* has suggested three 
possible energy transport processes. The plausi- 
bility of these three transport processes can be 
evaluated by comparing their predicted trans- 
port times with the experimentally measured 
decay time® of the light pulses. This has been 
done by Harrison,?° who concludes that the most 
probable mechanism is transport of energy due 
to resonance interactions between the molecules. 

Once the energy is at the solute molecule it has 
to be retained there long enough to be radiated, 
as otherwise it would be returned to the solvent 
where it would be quenched. This energy reten- 
tion®* can be explained in the following way. 
The excited solvent molecule must have an exci- 
tation level which is higher than some levels of 
the solute molecule. When these two molecules 
interact, energy will be transferred to the solute 
molecule exciting it to many levels. The energy 
irom those levels of the solute molecule which are 
higher than levels of the solvent molecules is 
probably immediately retransferred to the sol- 
vent and quenched, while the energy in the 
remaining levels of the solute molecule is re- 
tained. However, the retention of the energy by 
the solute molecules still does not insure the 
occurrence of fluorescence, as the excited solute 
molecules can still be quenched and these quench- 
ing processes are competing with the emission 
process. The homo-quenching effect depends 
again on the concentration of the solute mole- 
cules, while hetero-quenching with the solvent 
molecules, which can also occur, is independent 
of the solute concentration. This theoretical ex- 
planation agrees rather well with the observa- 
tions*-"' that the spectra from scintillating solu- 
tions are always the emission spectra of the 





20F, B. Harrison, thesis, Princeton University, (1951). 

20a The term ‘‘retention” is used instead of “trapping,” 
which is often used in the literature, in order to prevent 
confusion between this process and the process of trapping 
of electrons. 
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solute molecules and are independent of the 
spectra of the solvents. 

It is evident from the above discussion that 
the theory of liquid scintillation counters is still 
rather incomplete and that it does not predict 
what solutions will make good liquid scintillators. 
The only effect which is really predicted is the 
effect of solute concentration on the emission 
intensities. In order to derive this relationship 
and also to sum up in concise form the role played 
by the various phases of the whole phenomenon 
it is constructive to put all the concepts into 
mathematical form.’® Let P,=probability of 
homo-quenching of the excited solvent molecules; 
P,=probability of the energy from the solvent 
transferring to the solute molecule; 2,=number 
of excited solvent molecules produced per unit 
time; and my=number of excited solute molecules 
produced per unit time. Then 


Gea 

nz = {| ——— }n,=( — ; 
P,+P, 1+P,/P, 

However, every excited solute molecule will not 
radiate as its energy can also be lost by means of 
homo- and hetero-quenching. If P.= probability 
of light emission from the solute molecule; 
Pie=probability of hetero-quenching of the 
solute molecule; and P;.=probability of homo- 


quenching of the solute molecule, then the 
emitted light intensity is 





‘ ( F, ) Ny 
=n a = oF gememg 
NPAPwtPred 1+ (ProtPr)/P. 





(14+P./P,)(1+(ProtPre)/Pe) 


As P, and P,,. are proportional to the concentra- 
tion ¢, let 


(P,/P,)=k/c, (Pho/Pe) =k'c, (1+Pre/P.) =k". 


Substituting in the above expression one obtains 


c i 
heat 
c+k/ \k'c+k” 


This expression predicts that the intensity varies 
with the concentration in the following manner. 
At low c the light intensity should be propor- 
tional to the concentration; at higher concentra- 
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tions the function reaches a maximum and at 
still higher concentrations it gradually decreases 
as the denominator of the second factor becomes 
more dominant. This behavior of the light in- 
tensity as a function of concentration is verified 
by experiments'®:!* as can be seen in Fig. 1. 

In some cases small amounts of some other 
second solute decreases the light output of the 
solution.'!! This effect is attributed to absorp- 
tion of the light from the first solute molecules 
by the second solute molecules and consequent 
partial re-emission of the absorbed light. An 
important aspect of this phenomenon, from the 
point of view of making better counters, is the 
following: A second solute might decrease the 
light output of a solution only a little. However, 
if the wavelengths of the emitted light from the 
second solute are longer than the wavelengths of 
the light from the first solute, then the observed 
light from a solution of the two solutes will 
always have the spectrum of that solute which 
has the redder light. This fact is often utilized 
when one tries to match emission spectra to the 
spectral sensitivity of the light detectors and 
will be discussed in more detail in the experi- 
mental section of this paper. 


SOLUTIONS 


From the previous section and with the aid of 
experimental investigations we can now select 
the most favorable solutions for liquid counter 
application. While it is difficult (mainly because 
of the lack of understanding of the quenching 
processes) to predict what solutions would give 


RELATIVE INTENSITY 


CONCENTRATION OF ANTHRACENE IN XYLENE — GM/LITER 
Fic. 1. Relative light output of an anthracene-xylene 


solution as a function of anthracene concentration. The 
fluorescence is observed with a 1P28 photomultiplier. 


21 H, Schneider, thesis, Geissen (1950). 
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maximum light intensities, one can predict that 
one and only one specific concentration of solute 
in a solution will produce maximum intensity. 
This, of course, is due to the homo-quenching 
effect which comes in strongly with large con- 
centrations. By far the most complete investiga- 
tions of the responses of various solutions are 
those of Kallmann."-!6?2 Column 4 in Table II 
gives the relative emission intensities of various 
solutions (anthracene being used as a standard of 
1) as measured with a 1P28 phototube. This 
table is only a condensation of the extensive 
tables published by Kallmann, and only solutions 
which have a relative intensity of 0.20 or greater 
are given here. For special applications it might, 
of course, be desirable to use one of several 
solutions not mentioned in Table II. The relative 
merit of these solutions can probably be found in 
references 11, 16, 20, or 22. 

As the light emitted by the solutions is not 
directly observed by eye, it is necessary to con- 
sider the spectral response of the detector, the 
photomultiplier tube in this case, in order to 
pick a scintillating solution which will produce 
maximum response from the photomultiplier. 
For this reason column 5 of Table II gives the 
range of the emission spectra®'!.”? of the liquid 
phosphors. Looking at Table II one notices 
immediately that the best solvents are: xylene 
and phenylcyclohexane, while the best solutes 
are: terphenyl, anthranilic acid, diphenylhexa- 
triene, and a-naphthylamine. While Table I] 
gives the relative merits of solutions when 
used with a 1P28 (maximum spectral response 
~3400A), liquid scintillators are frequently used 
in conjunction with a 5819 phototube which has 
a large photosurface on the end of the tube. As 
the spectral response of the 5819 photomultiplier 
has its maximum at about 4800A and a half- 
width of about 1000A, it is desirable to increase 
the wavelengths of the emission spectra of the 
solutions if they are to be used with this tube. 
If one adds a small amount of a second solute, 
whose spectral range lies in a longer wavelength 
region than that of the first solute, it is possible 
to shift the spectral range of the fluorescence to 
that of the second solute while the relative light 
intensity of the solution behaves as if only the 
primary solute were present. For example di- 
~ 9H, Kallmann, Nucleonics 8, 32 (1951). 
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TABLE II.* Intensities of fluorescent light from solutions. 


Solute Solvent 





Phenyl-a-naphthylamine 
Terphenyl 

Terphenyl 

Terphenyl 
Diphenylhexatriene 
Phenyl-a-naphthylamine 
Terphenyl 

Anthranilic acid 
a-naphthylamine 
Phenyl-a-naphthylamine 
Terphenyl 

Anthranilic acid 
m-Diphenylbenzene 
Diphenylbutadiene 
a-naphthylamine 
3-naphthylamine 
Phenyl-a-naphthylamine 
lerphenyl 

Carbazole 
Diphenylhexatriene 
Terphenyl 

Anthranilic acid 
rerphenyl 

Anthranilic acid 
m-Diphenylbenzene 
Diphenylhexatriene 
a-naphthylamine 
Phenyl-a-naphthylamine 
lerphenyl 


benzene 

benzene 

p-cymene 
p-dioxane 
diphenyl oxide 
diphenyl oxide 
diphenyl oxide 
phenetole 
phenetole 
phenetole 
phenetole 
phenylcyclohexane 
phenylcyclohexane 
phenylcyclohexane 
phenylcyclohexane 
phenylcyclohexane 
phenylcyclohexane 
phenylcyclohexane 
phenylcyclohexane 
phenylcyclohexane 
tetrahydronaphthalene 
toluene 

toluene 

xylene 

xylene 

xylene 

xylene 

xylene 

xylene 


(gms/liter) 


Number of 
hydrogen 
nuclei 
X1022/cem? 

of solvent 


Concentration Relative Spectral 


intensity range (A) 


0.28 
0.37 
0.24 
0.23 
0.37 
0.48 
0.52 
0.22 
0.22 
0.28 
0.32 
0.28 
0.24 
0.24 
0.33 
0.28 
0.42 
0.52 
0.20 
0.36 
0.23 
0.20 
0.42 
0.23 
0.20 
0.24 
0.24 
0.31 
0.46 


N 
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3300-4000 
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3600-4200 


3700-4500 
3200-4100 
4300-3800 
3900-5500 
3500-4000 


3200-4000 
4000-5500 
3600-4400 
3200-4500 
3250-4000 
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See reference 19. 


phenylhexatriene and fluoranthene have both 
emission spectral ranges between 4000 and 
5000A" and are thus good solutes for solutions 
which are intended to be used with 5819 photo- 
multipliers. Considering all the available informa- 
tion, one can see for example that a terphenyl- 
phenylcyclohexane solution with a slight con- 
tamination of diphenylhexatriene should give 
the best response to ionizing radiation when used 
with a 5819 tube. (Relative intensity ~0.70.?%) 
These organic liquid scintillators can be used 
as fast neutron detectors if one utilizes the proton 
recoils produced by the neutrons. A number of 
investigators have used this principle.**~°* In 
this case another factor has to be considered in 
determining the suitability of a solution and 
that is the number of hydrogen nuclei per cm’. 
Column 6 in Table II gives this information and 
it can be seen that the various solutions do not 
vary vastly in their concentration of hydrogen. 


*8H. Kallmann (private communication). 

* Falk, Poss, and Yuan, Phys. Rev. 83, 176 (1951). 

*G. R. Keepin and R. H. Lovberg, private com- 
munication. 

% Jastram, Benade, Cleland, and Hughes, Phys. Rev. 
81, 327 (1951). : 


It should be pointed out that xylene, toluene, 
or benzene are very toxic. Poisoning can develop 
through excess absorption of the liquids through 
the skin or through the breathing of the vapors. 
For this reason caution should be used when 
handling these liquids. 


CONSTRUCTION AND ELECTRONICS 


Before discussing the response of liquid scin- 
tillation counters to various types of radiation, 
we shall mention some aspects of the construction 
of such counters and the circuitry involved in 
recording the signals.* Of the many different 
phototubes available in'this country the following 
are most frequently used : the RCA types 931-A, 
1P21 (name given to best 10 percent of 931-A), 
1P28 (same as 1P21 but made with an envelope 
of different glass) and 5819. All of these tubes 
except the 5819 have internal photocathodes and 
the 5819 cathode is by far the largest one of all 
the tubes. In the following discussion we shall 
have in mind the RCA 5819 photomultiplier 
which is often used for general scintillation 
counting. It has a 14-inch diameter photocathode 
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Fic. 2. One method of mounting liquid counting cells 
on RCA 5819 photomultipliers. 


on the inner glass surface at the end of the bulb 
and 10 electrostatically focused multiplying 
stages. The glass end of the bulb is slightly 
convex so that a flat crystal mounted on it does 
not touch the surface at all points. 

The actual construction of a liquid scintillation 
cell depends on the use to which it will be put. 
If one is interested in counting without energy 
resolution one can actually use very large con- 
tainers® (such as Dewar flasks) and just sub- 
merge the photomultiplier surface in the liquid. 
This is not the most efficient use of the solution 
as a considerable amount of the fluorescence is 
lost by absorption in the liquid and in the re- 
flecting surfaces. In a great many applications, 
however, it is desirable to use the scintillator in 
such a way as to utilize maximum counting 
efficiency as well as maximum energy resolution. 
For these cases it is advisable to have liquid 
cells which can be attached to the end or sides of 
the phototube. One way of doing this is by 
putting the liquid in a glass cell and coupling the 
latter optically to the photomultiplier surface 
through a layer of Canada balsam. However, 
the glass surface of the cell between the liquid 
and the photocathode not only absorbs a con- 
siderable fraction of the light but it reflects light 
too. To overcome this deficiency it is advisable 
to place the scintillating liquid directly in contact 
with the photomultiplier surface. The one prob- 
lem in this procedure is to find a suitable cement, 
that is a cement which will hold the container 
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securely to the photomultiplier surface and yet 
will not dissolve in the organic fluid of the cell. 
As a phototube cannot be heated to more than 
50°C only cold-setting cements can be used. Two 
cements of this type which have been used fre- 
quently with liquid counters are Sauereisen and 
Araldite Cement (CIBA). The disadvantage of 
Sauereisen is that contact with moist air will 
slowly soften it. However, there are several ways 
to overcome this difficulty: for example Falk 
et al.* used a technique where the glass cell was 
first cemented on with Sauereisen and after this 
cement had dried the joint was covered with a 
layer of sodium silicate (water glass). The whole 
process was then repeated. This procedure pro- 
duces a strong joint which will stand up for 
months. Kallmann® uses a technique where the 
glass cell is cemented to a metal ring which then 
in turn is cemented to the phototube, the cement 
also being Sauereisen. This method has the 
advantage that only a very small surface area is 
exposed to air. Figure 2 shows a sketch of one 
of these two mounting methods. It should be 
mentioned that o-ring type seals or gaskets can 
be used. However, they do not always make a 
very steady support and great care has to be 
taken that the o-ring, which will usually be in 
contact with the scintillating solution, does not 


- STRAY 
a CAPACITANCE 
' 


CATHODE 


Fic. 3. Photomultiplier wiring diagram with positive 
end of power supply grounded. R=1M, C=0.01 uf. Cathode 
follower stage using 6J6 twin-triode, both units in parallel, 
also is shown. 


27H. Kallmann, U. S. Signal Corps, Progress Report 2, 
New York University (1951). 












poison the counter (gaskets made of Resistoflex 
have been used very successfully). 

As it is desirable to collect as much of the 
fluorescent light as possible on the photomulti- 
plier surface, it is important to have good re- 
flecting surfaces in the liquid cell. In general it is 
advisable to use metallic reflectors instead of 
diffuse ones (such as magnesia or titanium oxide 
powder). As for the metallic reflectors, aluminum 
is the best all-around metal. Its reflectivity below 
4100A is far superior to that of silver or any 
other metal. Above 4100A silver’s reflectivity is 
5-10 percent better than that of aluminum.”® 
However, this comparison is only valid with 
fresh metallic surfaces, and while the aluminum 
remains untarnished for long periods the silver 
reflectors will tarnish quite rapidly. The authors 
have found that glass cells coated on the inside 
with an evaporated layer of aluminum make 
excellent reflecting containers. 

Since the gain of a photomultiplier varies 
markedly with voltage, it is advisable to obtain 
the dynode voltages from a regulated power 
supply. The voltages can be conveniently taken 
from taps on a voltage divider connected across 
the supply. The higher the resistance of the 
voltage divider, the less current it will draw, and 
the better will be the regulation of the power 
supply. However, if the signal current passing 
through the tube becomes comparable to the 
voltage divider current, the dynode potentials 
will fluctuate and hence the gain will be a func- 
tion of signal size. In practice, we have used 
l-megohm dynode resistors as shown in Fig. 3 
bypassing the last three dynodes with 0.01-yf 
condensers to maintain their potentials during 
large pulses. We have operated the photomulti- 
pliers at an over-all potential of 1000 volts, or 
about 90 volts per stage. 
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Fic. 4, Block diagram of the electronic equipment used 
in conjunction with a scintillation detector. 
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28J. Strong, Procedures in Experimental 
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(Prentice-Hall, Inc., New York, 1938), p. 375. : 
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Fic. 5. Equivalent output circuit of photomultiplier. 


It is convenient to have a power supply in 
which the positive terminal can be grounded. 
The anode of the photomultiplier can then be 
operated at ground potential as shown in Fig. 3. 
If the photomultiplier is to be operated with a 
negatively grounded power supply, the anode 
will be a high potential and a good quality 
coupling condenser is then required to pass the 
signal and not introduce any spurious pulses. 
A low loss condenser across the power supply is 
desirable in this case so that all of the signal 
voltage will appear across the load resistor. 

Unlike the output of a Geiger counter, which 
consists of more or less uniform pulses, the output 
pulses of a scintillation counter form a con- 
tinuous distribution in size. The counting rate 
then has significance only if it refers to the 
number of pulses above a certain level, or be- 
tween two levels, other conditions remaining 
fixed. The pulses from the photomultiplier are of 
the order of several tenths of a volt and have 
to be amplified a few hundred times to bring 
them up to the range where amplitude discrim- 
inator circuits operate well. The discriminator 
output pulses can then be counted with a scaler. 

In some applications, such as in accelerator 
experiments, where the scintillation counter may 
be separated by a large distance from its ampli- 
fier and scaler, a cathode follower circuit, such as 
shown in Fig. 3, provides a low impedance out- 
put and can drive cables 100 feet long without 
much loss in signal. A block diagram of a typical 
arrangement is shown in Fig. 4. 

A further point that deserves mention is the 
time constant of the output of the photomulti- 
plier circuit. If we neglect the time required for 
the light to build up and for the photomultiplier 
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to respond, we can assume that the current 
reaching the anode of the photomultiplier due 
to a scintillation is given by Joe~“/7, where Io 
would be proportional to the number of photons 
incident on the photocathode and 7 is a decay 
constant characteristic of the phosphor. This 
current pulse will be fed into the anode resistance 
R and stray capacitance C as indicated in Fig. 5. 
Under these conditions, the output voltage will be 


IyR 
~ (RC/t)—-1 


V (e~#/RC — ett), (1) 


For the special case RC =r, 
V=RI(t/r)e—"". (2) 
Equations (1) and (2) have a maximum value 
Vinax = RI ge—*!*, (3) 
where for Eq. (2), 


to=T, (4) 
and for Eq. (1), 


. 
to = ————— log. RC/r. (5) 
C 


1—7/R 


For a given 7, if we vary RC from 0.27 to 
10007, the rise time f) only varies from 0.47 to 
6.97. Ameasurement of 7 for terpheny] in toluene® 
gives r=(2.2+0.3) X10-* sec so that even for 
large values of RC, the rise times will be faster 
than those of most conventional amplifiers. As 
can be seen from Eq. (1), the decay of the pulse is 
determined by RC when it is large compared with 
t, which will usually be the case for a liquid 
scintillator. In the circuit of Fig. 3 with R=10K, 
RC is in the range of 0.1 to 0.2 usec which permits 
a high counting rate without pile-up. In applica- 
cations where one is dealing with a pulsed accel- 
erator beam of short duration, or time-of-flight 
coincidence measurements, a smaller value of 
RC is usually necessary. If larger time constants 
can be tolerated, some increase in signal can be 


realized by increasing R. As can be seen from. 


Eqs. (3) and (5), this increase is not proportional 
to R, since a larger R results in a smaller value 
of the exponential. 


COUNTER RESPONSE TO RADIATION 
AND APPLICATIONS 


In discussing the response of a liquid scintilla- 
tion counter to various types of radiation, two 
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properties are of interest: efficiency and energy 
response. Charged particles will be counted as 
long as they are of sufficient energy to give rise 
to more than a certain number of photons. 
This number will depend on the light collection 
efficiency, the photomultiplier noise, and the 
noise level of the rest of the system. In this 
respect, the liquid scintillation counter does not 
offer an advantage over gas counters or crystal 
scintillation counters which can likewise count 
charged particles with 100 percent efficiency. 

In regards to energy response, the liquid scin- 
tillators are similar to the organic crystals in that 
for a given energy loss in the phosphor, a lightly 
ionizing particle such as an electron will give 
rise to a larger pulse than will a more densely 
ionizing particle such as a proton or alpha- 
particle. The reasons for this behavior can be 
understood in terms of the ‘‘excitation’’ theory 
of Bowen ef al. Birkes?® has shown that when 
this theory is actually applied to fluorescent hy- 
drocarbons it can quite accurately predict the 
relative pulse heights produced by various 
ionizing particles at different energies. The in- 
organic scintillators such as thallium activated 
sodium iodide give a pulse proportional to energy 
independent of the specific ionization of the 
particle.*° However for a given type of particle 
the response of a liquid scintillator is at least 
roughly proportional to energy.” 

Gamma-rays and neutrons are counted in- 
directly, the former by the electrons produced in 
photoelectric, Compton, or pair production 
processes, and the latter, by the recoil protons 
they give rise to. The counting efficiencies are 
thus determined by the cross sections for these 
effects. In accordance with the remarks in the 
preceding paragraph, a gamma-ray of a given 
energy will give a larger pulse than will a 
neutron of the same energy. 

In determining gamma-ray energies by the 
pulse heights they give rise to, the inorganic 
crystal phosphors have the advantage of con- 
taining high Z elements compared to the organic 
crystals and the liquids. The photoelectric cross 
section, for example, varies as Z°. Since the full 
energy of the gamma-ray, less the binding energy 

29J. B. Birkes, Proc. Phys. Soc. (London) 64A, 874 


(1951). 
30 Franzen, Peele, and Sherr, Phys. Rev. 79, 742 (1950). 












of the electron, appears in the photoelectron, a 
measurement of its pulse size gives the energy 
of the gamma-ray. The Compton process, pro- 
portional to Z, gives a continuum of electron 
energies and unless the phosphor is large enough 
so that there is a good chance of the degraded 
photon being absorbed by a photoelectric proc- 
ess within the counting volume, the full energy 
of the gamma-ray will not be given up. A meas- 
urement of the total energy of a pair produced 
by a higher energy gamma-ray also enables its 
energy to be determined. Here again the in- 
organic crystals are favored since in addition to 
pair production being proportional to Z?, their 
higher density causes the electron ranges to be 
correspondingly shorter so that there is less likeli- 
hood for them to escape from the phosphor. It is 
for these reasons that liquid or organic crystal 
phosphors have not been used in gamma-ray 
spectroscopy. 

In the case of beta-ray spectroscopy, the 
problem of conversion is not present so that the 
organic crystals and liquids can be used. Much 
work has been done by P. R. Bell and his 
associates on this subject using anthracene 
crystals. Liquids offer an advantage if the 
maximum electron range exceeds the dimensions 
of the crystals available. In addition, they offer 
some further possibilities for improving geom- 
etry. For example, Raben and Bloemberger*! 
have experimented with immersing the beta-ray 
source in the liquid, thus utilizing a 47 solid 
angle. 

For fast neutron counting, greatly increased 
efficiencies have been realized using liquid scin- 
tillators. The high efficiencies result from the 
large usable thicknesses and the high concentra- 
tion of hydrogen nuclei. Gas counters using a 
hydrogenous material as an irradiator or having 
a hydrogenous gas filling have an effective 
number of grams of hydrogen per cm* which is 
less by a factor of several hundred than that 
available in a liquid scintillator. The inorganic 
phosphors lack hydrogen while the organic 
crystals cannot be obtained in great volumes. 
Thus neutron counters have generally a low 
counting efficiency. W. F. Hornyak,” however, 


(1951): S. Raben and N. Bloemberger, Science 114A, 363 
1 


32 W. F. Hornyak, Rev. Sci. Instr. 23, 264 (1952). 
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Fic. 6. Integral bias curves for neutrons and y-rays using 
a 1-inch terphenyl-xylene cell with a glass bottom. The rela- 
tive counting rates of the three curves are not significant. 


has achieved a moderately high neutron counting 
efficiency by molding a powdered form of zinc 
sulfide with Lucite, the recoil protons which 
produce scintillations in the ZnS coming from 
the Lucite. Opacity limits the useful thickness of 
this detector to about a centimeter. Compared 
to all these detectors, the liquid scintillator can 
give higher counting efficiencies, which can be as 
high as 50 percent.”6 

Since recoil protons of all energies up to the 
maximum neutron energy are produced with 
equal probability in a hydrogenous liquid, the 
integral bias curve for neutrons detected in a 
liquid scintillator would be expected to be a 
straight line of negative slope intersecting the 
pulse height axis at a point corresponding to 
the maximum neutron energy. Figures 6 and 7 
show bias curves™ for some different energy 
neutrons and gamma-rays. Straight lines, when 
transformed to the logarithmic scale, can be made 
to fit the neutron bias curves fairly well, as 
shown in Fig. 7. 

The disadvantage of the liquid scintillators for 
fast neutron counting lies in their response to 
gamma-rays. It can be seen from Fig. 6 that the 
maximum pulse height for 14-Mev neutrons is 
only about 2.3 times that for 1.3-Mev gamma- 
rays which in turn give pulses about the same 
size as 4-Mev neutrons. Consequently, if gamma- 
rays comparable in energy to that of the neutrons 
are present, they cannot be biased out and can 
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thus seriously limit the use of liquid scintillators 
in neutron counting applications. Harding* has 
overcome this problem by making an emulsion 
of potassium iodide crystals in a-bromonaphtha- 
lene. The potassium iodide was ground up into 
having a diameter equivalent to the range of the 
recoil protons produced by 2.6-Mev neutrons in 
the organic liquid. The range of electrons pro- 
duced by gamma-rays of the same energy 
(2.6 Mev) in the potassium iodide, was much 
larger than the diameter of the crystals. Con- 
sequently all the “‘neutron” energy was lost in the 
crystal while only part of the “gamma-ray” 
energy was dissipated. As a result 2.6-Mev 
neutrons produced pulse heights slightly larger 
than those produced by 2.6-Mev gamma-rays. 
The serious disadvantage of this type of counter, 
however, isitsover-all neutron counting efficiency 
which is less than 1 percent and the fact that a 
special counter has to be designed for each 
energy. It is only possible to use this method 


THEORETICAL 
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NEUTRONS 
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NEUTRONS 
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° 
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10 LZ 
PULSE HEIGHT-VOLTS 


Fic. 7. Relative pulse-height distribution produced by 
various energy neutrons in a 1-inch terphenyl-xylene cell 
with the scintillating solution being in direct contact with 
the 5819. The relative counting rates of the different curves 
are not significant. The dashed curve is the expected pulse 
distribution assuming only one collision of the neutron 
in the liquid. 


33 G. N. Harding, Nature 167, 437 (1951). 


Cc. E. FALE AND HF. 
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if the neutrons have energies higher or at least as 
high as that of the gamma-rays. The zinc sulfide- 
Lucite detector previously mentioned is not 
subject. to this limitation. However, a reaction 
such as deuterons on tritium produces 14-Mev 
neutrons without the presence of gamma-rays. 
The reaction takes place at low bombarding 
energies so that there is no background of 
neutrons or x-rays from the accelerator itself. 
In this case, a liquid scintillator is very well 
suited as a neutron detector and can be used in 
lengths of several inches to get efficiencies of 
about 20 percent. Such a detector has been used, 
for example, in a measurement of the -p cross 
section at 14 Mev.* The d-d reaction is also free 
from gamma-rays but here the higher bombard- 
ing energy may in some cases produce an un- 
desirable x-ray background from the accelerator. 

The liquid counters described above are not 
suitable for detecting thermal neutrons. How- 
ever, Draper*® has immersed a thin layer of 
powdered B,O; in a terphenyl-toluene solution 
and has achieved thermal neutron counting 
efficiencies of approximately 50 percent. Re- 
cently Muehlhause and Thomas** have shown 
that triethyl-, trimethyl-, or tributyl-borate 
mixed with a terphenyl-phenylcyclohexane solu- 
tion will make a slow neutron detector of almost 
100 percent efficiency. The borate does not 
quench the fluorescence and if enriched in B!° 
the lifetime of the neutron in the solution is 
~0.4 microsecond. 

Liquid counters have been used in cosmic- 
ray research where their large counting volumes 
are desirable. Thus Ney and Thon* sent a liquid 
counter in conjunction with a Geiger counter 
telescope aloft in a balloon, recording pulses due 
to the passage of heavy primaries, and it was 
possible with this device to relate the pulse size 
to the Z of the primary. Harrison and Reynolds** 
utilized the speed of a terphenyl-toluene counter 
as well as its large available volume to determine 
the lifetime of u-mesons in copper and antimony. 


’ x“ _ Salant, Snow, and Yuan, Phys. Rev. 87, 11 
1952). 
% J. E. Draper, Rev. Sci. Instr. 22, 543 (1951). 

36 C, O. Muehlhause and G. B. Thomas, private com- 
munication. 

37 EF, P. Ney and D. M. Thon, Phys. Rev. 81, 1069 (1951). 

38 Harrison, Keuffel, Godfrey, and Reynolds, Phys. Rev. 
86, 616 (1952); Harrison, Keuffel, and Reynolds, Phys. 
Rev. 83, 680 (1951). 


















Experiments performed with meson or proton 
beams from the large synchro-cyclotrons use 
liquid scintillators extensively as detectors be- 
cause of their speed, size, and ease of preparation. 
Figure 8 shows an example of this type of counter 
used in these high energy experiments. It shows 
an eight-inch diameter liquid counter designed 
by Messrs. H. Loar, R. Durbin, and Dr. J. 
Steinberger, which has been used at the Nevis 
cyclotron of Columbia University. This type of 
counter is usually used as the final counter in a 
counter telescope, while the other elements of the 
telescope are either thin liquid counters of smaller 
diameter, or thin stilbene crystals.*:4° 


CONCLUSIONS 


The liquid scintillation counter as a research 
tool has some very distinct advantages over other 
detectors including scintillation crystals. The 
rise time of its pulse of ~2X10-* sec, as com- 
pared to a rise time of ~5X10~-° sec* of the 
best crystal scintillator (diphenyl-acetylene), 
makes it one of the fastest detectors in use and 
thus very useful in fast coincidence work (espe- 
cially as pulse heights can be as high as 60 percent 
the pulse height of anthracene). Its high counting 


39 — Loar, and Steinberger, Phys. Rev. 83, 646 
(1951). 


40 Chedester, Isaacs, Sachs, and Steinberger, Phys. Rev. 
82, 958 (1951). 


41 National Radiac Incorporated, Bulletin S-2 (1951). 





The last point that must be discussed . . 


LIQUID SCINTILLATION COUNTERS 


. ts what is meant by the “universe.” Accepting the 
postulate that the velocity of light is the maximum velocity at which influences can be propagated, 
then it is evident that all events on or within our past light cone may affect us, and they certainly 
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Fic. 8. Large liquid scintillation counter of the type 
used frequently in experiments with particles from large 
accelerators. (Counter made and used at Nevis Laboratory, 
Columbia University.) 


efficiency for fast neutrons makes it invaluable 
as a fast neutron detector, especially when small 
neutron intensities are involved (however, only 
when y-ray intensities are absent or small). The 
counters are easy to handle and easy to prepare. 
Furthermore, it is very easy to make counting 
volumes of much larger sizes than can be 
produced with crystals. The latter feature makes 
the liquid counter very useful in cosmic ray and 
meson research. Thus, there are and can be many 
very distinct applications of this tool in research 
especially with the development of better count- 
ing solutions and faster electronic circuits. 

The authors wish to thank Professors H. 
Kallmann and G. T. Reynolds for their helpful 
discussions and Dr. E. O. Salant for his constant 
advice and aid. 


must be included within our definition of the universe. If we extend this to include all events which may 
affect us at some time in the future and all events which have been or will be affected by us, then we 
arrive at the largest set of events that can be considered to be physically linked to us. It is this set 
that is usually considered to constitute the universe. Some authors, however, consider a different 
set, viz. the largest set to which our physical laws (extrapolated in some manner or other) can be 
applied. This set may include points that have never been and will never be physically linked to 
us and may exclude points which are so linked. The physical significance of this ‘‘universe’’ is 
therefore not very clear, but it should be remembered that even with the definition adopted above the 
actually possible observations concern only a very limited part of the universe—Cosmology by H. 


Bonp1 (Cambridge University Press, 1952), p. 10, by permission, 
















































































































































On the Theory of Dipole Interactions with Metals 


S. J. Czyzax* 
University of Detroit, Detroit, Michigan 
(Received January 18, 1952) 


The classical and quantum-mechanical interactions of molecular dipoles and metallic sur- 
faces are derived and numerical calculations given for several molecule-metal pairs. The 
results represent the permanent dipole contributions to the total van der Waals interactions. 
The bearing of the dipole interactions on the problem of adhesion is discussed. 


HERE are many practical applications!” in 
which interactions of organic molecules 

with metallic surfaces play an important part. 
Adhesion, lubrication, friction, wetting, and 
physical absorption are only a few examples in 
which this type of interaction may be involved. 
The pertinent forces in this interaction are 
most often the well-known van der Waals forces. 
Van der Waals forces may be ascribed to the 
following three more or less distinct effects: 
orientation of permanent dipoles, induction, and 
dispersion ; the last being an interaction between 
instantaneous dipoles. This means that, while the 
time average dipole moment vanishes, in this 
case there is at any instant an interaction which 
turns out to be always attractive. It is generally 
found that the dispersion effect does not vary too 
greatly from one molecule to another if referred 
to individual atoms; whereas, the dipole orienta- 
tion and induction effects are dependent upon the 
dipole moment of the molecules and may in the 
case of nonpolar molecules be entirely absent. 


Fic. 1. Diagram of the position of the dipole and image. 


* On leave; Lt. Col. USAF Flight Research Laboratory, 
Wright Field, Dayton, Ohio. 

a Ig Ferguson and E. F. Stephan, A.E.S. Research 
Repos Be, SOE, SVS bhnenes, J ‘Colloid Sci. 2, 163 
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2W. D. Harkins, AAAS Publication No. 7, 19 (1939); 
AAAS Publication No. 20, 40 (1943). 


Thus, the latter two effects may be considered as 
being superimposed on a more or less constant 
background of dispersion forces. It is observed 
experimentally that such phenomena as adhe- 
sion, lubrication, friction, wetting, etc., depend 
upon the polarity of the molecules. It is the 
purpose of this paper to examine the effect of the 
polarity of the molecule on the interaction energy 
with metallic surfaces. 


CLASSICAL TREATMENT OF THE INTERACTION 
OF A DIPOLE WITH THE METAL SURFACE 
The classical treatment of dipole-metal inter- 
actions was first worked out by Lorentz and 
Landé, Magnus, Jacquet, Lennard-Jones, and 
others*-® in connection with the adsorption of 
vapors on metals. The present development 
follows much along the same lines. The inter- 
action energy between two dipoles of moments 


u; and wu: will be the potential energy of dipole wu. 
in the field of uw: 


ui- 
V.=we- VY 


In the scalar notation this equation takes the 
well-known form 


(1) 


U1- as 3u1-Ru2-R 
—)- 7 “RS 


12 
Ve= ~~“ cos@; cos#.—sin@; sinf, cosg), (2) 


where 6; and 62 are the angles of the dipoles with 


(1920). Lorentz and A. Landé, Z. anorg. Chem. 125, 47 
1922 

4A. Magnus, Z. physik. Chem. A142, 401 (1922). 

5 E. Jacquet, Fortschr. Chem., Physik, u. physik. Chem., 
Series B 18, 147 (1925). 

®F. London and R. Ejisenschitz, Z. Physik 60, 520 
(1930); 63, 245 (1930). J. E. Lennard-Jones, Trans. 
Faraday Soc. 28, 334 (1934). 

7H. Margenau, Revs. Modern Phys. 11, 1 (1939). 

8R. H. Fowler and E. A. Guggenheim, Statistical 


Thermodynamics (Cambridge University Press, Cambridge, 
1939). 
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ON THE THEORY OF DIPOLE 


respect to R, and gis the angle between the planes 
defined by the two dipoles and R as shown in 
Fig. 1. Thus the interaction energy in Eq. (1) is 
that of any two randomly oriented dipoles. Since 
we desire the interaction between a dipole and 
metal surface, Eq. (2) becomes by using Kelvin’s 
method of images 

2 


Ve= =~—fi +o, (3) 
R3 


where 0:= —02=80, wi=u2=u, g=0. This equa- 
tion now represents the interaction between a 
dipole and its image. 

From Eq. (3) it can be seen that the maximum 
interaction occurs when 6@=0, i.e., with the 
dipoles perpendicular to the surface. When many 
dipoles are involved in a system, which is the 
case with actual molecules, the dipoles will have 
a tendency to line up perpendicular to the 
interface. However, this tendency is overcome by 


R3 


Simplification and integration leads to 


— . 1 e* 
V 


=—]14—/———_--1 


2xr 1 
f dx 
0 


h=2/kTR? 


:, (8) 


where 


x=cosé. 


QUANTUM-MECHANICAL TREATMENT OF THE 
INTERACTION OF A DIPOLE WITH THE 
METAL SURFACE 


The metal is pictured as being built up of posi- 
tive and negative charges.®—!* Consider the inter- 
action between the molecule with its coordinates 
at the origin and an element of metal a distance r 
from the origin. By an element is meant a volume 


®H. Margenau and W. F. Pollard, Phys. Rev. 57, 557 
(1940). 


( 10 Oy. Margenau and W. F. Pollard, Phys. Rev. 60, 128 
1940). 

11 J. Bardeen, Phys. Rev. 58, 727 (1940). 

12 Prosen, Sachs, and Teller, Phys. Rev. 58, 1066 (1941). 
13 E, Prosen andjR. G. Sachs, Phys. Rev. 61, 65 (1941). 
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thermal agitation which will tend to keep them 
randomly oriented. The latter effect prevails at 
high temperature and, in the limit, complete 
randomness will prevail. If we average the inter- 
action energy over all orientations, we get 


2 


v.=- 





: [J (1-+c0s%) sinédédg, (4) 
T 


which upon integration gives 


V.= —4y?2/3R3 (5) 
and at very low temperatures approaches 
V.= —2p?/R’. (6) 


At intermediate temperatures the value of V will 
lie somewhere between 2y?/R* and 4y?/3R°. 
If the interacting molecules are assumed to 


obey Boltzmann statistics the following result is 
obtained : 


/ f 1-+cos*e) exp[u?(1+cos?6)/RTR®] sinédédy 
be 


(7) 


f exp[u?(1+cos?0)/kTR®] sinédédg 


of smaller dimensions than the wavelength of 
radiation corresponding to the predominant reso- 
nant frequencies of the molecule, but sufficiently 
large to possess the properties of the metals. The 
total interaction can be obtained by summing up 
the elements, and the well-known equations for 
polarizability may be employed in solving the 
problem. The interaction energy may then be 
written in the well-known form as follows: 
‘ ' 


V= ——(XMXO+4 VO VO—27MZ), (9) 
r3 


where X‘, X®), etc., are the components of the 
vectors for the molecule and metal, respectively. 
Let ¥(1, 2) =¥(1)¥(2) represent the unperturbed 
state function for calculating the interaction 
energy, where y¥(1) represents the state of the 
molecule and ¥(2) represents the state of the 
metal. If the molecule is represented as a rigid 
rotator, its wave equation is 


y(1) _ R(r1)8(¢1) O(41) 


= R(r;)P(lim,; cosé)em™*, (10) 
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In an unperturbed metal the charge is uni- 
formly distributed to a very high order of ap- 
proximation which accounts for the fact that 
there is no permanent dipole moment. Thus the 
first-order expression for the interaction energy 
vanishes, AV,=0, since it depends upon the 
product of the dipole moments of the molecule 
and of the metal. The interaction energy is then 
given by: 


| (Zm1; @2|UO| Aiur; be) |? 


—————— ee, (11 
E(a;) — E(b:) + E(a2) — E(b2) 


where each state 6 depends upon the quantum 
number of the molecule (5;) and metal (2). If 
Eq. (11) is averaged over all the magnetic 


quantum numbers of the molecule, then Eq. (11) 
becomes 


1 
 —Wh+1 


qm 


| (Zim; @2|V| Aiur; be) |? 


ninamapemnmatie,, ite 
m1 A141,b2 E(1,) — E(A1)+ E(a2) — E(b2) 


tl ¢ 


7 214,+1 re my A1H1,b2 


CZYZAK 


Since the molecule is regarded as a linear rigid 
dipole and its wave equation is of the form 
¥(1) = R(r1)6(¢1)0(01), the matrix Ua may be 
expanded and written in spherical coordinates. 
By these operations the dipole moment of the 
molecule turns out to be: 


pae f RorsRedr, (13) 


and now the matrix for V is: 


(Lym; @2|U| Ayr; be) 
= {(Limy| £1| Aiur) (@2| eX | be) 
+ (Lym1|n| Aru) (@2|e Y | be) 
—2(dym1| | Aim) (@2|eZ | b2)}, (14) 


where £é, 7, ¢, represent the respective trigono- 
metric functions. Since the molecule is regarded 
as a rigid dipole and fixed in position, only the 
rotational energy exists, no change in vibrational 
energy occurring. Generally the rotational ener- 
gies are small; especially small are the changes in 
rotational energies of the molecule as compared 
with the changes in energies within the metal.!+1° 
Hence E(a2) — E(b2)>E(l,) —E(A1) so that E(J;) 
— E(d,) may be neglected, and V4, m becomes: 


(Lym, | E| Axper) (@2| eX | be) + (Lyi | | Arr) (@2| e VY | be) 
E(a2) — E(b2) 


ly 1 1/1 2 2°) be $ 
aa my| | drp1) (a2 | eZ | be) (5) 
E(a2) — E(b2) 


Upon squaring the terms it is found that the cross-terms drop out, and if the metal is assumed 


to be isotropic, then Eq. (15) reduces to: 


2 


” mh4+17 


By evaluating the matrix, the expression becomes 


" 2 uw? _ | (a2|eR | be) |? 
V 


ye. ae 


3 18% ©=E(bs)—E(as) _ 


The polarizability for the metal in quantum- 
mechanical terms, where a2 represents the ground 


. 1 4 | | | (Zama | E| Kaw) |?+ | | 
m1 A1H1,62 


xj oe] (16) 
E(a:)—E(b:) } J 


state, is: 
2 _ | (a2|eR | be) |? 


as = . 
3 be E(b2) — E(a:) 
From this, Eq. (17) becomes 


(18) 


2 


- B 
Ves = ————"@ig. 
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1.26 


1.26 


1.26 


1.26 


1.26 


1.26 


1.745 


1.745 


1.745 


1.745 


1.745 


1.745 


1.745 


1.745 





1.54 


1.44 





1.50 


1.84 


1.68 


1.70 


0.26 


1.15 


0.60 


1.84 


1.68 


0.26 


1.45 


0.60 


i 





TABLE I. Metal-dipole interactions. 


Classical energy 


—4y2/3R3 
(ergs) 


Lat 3<i0"> 


2.31 X10" 


122 <i 


1.09 x10-% 


1.08 X10-% 


0.368 X 1078 


0.425 X10-" 


0.408 X 107" 


0.846 X 1078 


1.40 x107% 


0.790 X 10-8 


0.723 X10-* 


0.788 X 10-% 


0.253 X 10-8 


0.318 X 10-4 


0.302 X 10-8 








—2p?/R3 
(ergs) 


2.05 107% 





3.46 X10-* 
1.83 107% 
1.64 X10-* 
1.63 X10-% 
0.553 X 1078 


0.638 X 10-4 


0.612 < 10-4 


27 xiG = 
2:10 X10 
1.18 K107% 
1.08 X10~* 
1.18 Ki0* 
0.379 X 1078 


0.476 X10-" 


0.453 X 1078 


Q.M. energy 
—p?/12r8 
(ergs) 


0.685 X 10-8 





1.16 X10-% 


0.611 10% 


0.545 x 1078 


0.54 <X10-% 


0.184 X 1078 


0.212 X10-" 


0.204 x 10-4 


0.423 X 10-8 


0.700 X 10-8 


0.395 X 10-8 


0.362 X 10-8 


0.394 X 10-% 


0.127 X 10-8 


0.159 X 10-"4 


0.151 10-8 


INTERACTIONS WITH METALS 
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Classical 
energy 
—p?/R3Ba 
(ergs) 


x 107% 
x 1078 
x<10-%* 
x<10-8 
x10-% 
10-8 
xX 107% 
x<10-8 
x10-* 
10-3 
x 10-18 
10-8 
x10-% 
x 10-% 
xi¢-* 
04x 10-8 
0.386 X 10% 
0.374 X 107 
0.430 X 10-" 





te ee et ND Gn Ga BUD ON DOO 
BSON’W WON ONKHMANTON NINN 


oO 


0.428 X 10-4 
0.426 X 10-4 
0.453 X 10-4 


0.429 X 10-4 
0.418 X 10-" 
1.04 X10-* 
0.943 X 1078 
0.898 X 10-8 
1.86 X10-* 
1.67 X10-" 
1.54 <X10-" 
0.982 X 10-8 
0.870 x 10-18 
0.830 X 10-8 
0.857 X 10-8 
0.792 X 10-8 
0.755 X 10-8 
0.802 X 10-% 
0.798 X 10-1 
0.790 X 107 
0.270 X 10-8 
0.268 X 1078 
0.258 X 10-8 
0.321 10-" 


' 


0.320 X 10-™ 
0.319 X 10-# 
0.325 X 10-8 


0.313 X10-" 
0.305 X 10-8 





® Data obtained from Chemical Physics (McGraw-Hill Book Company, Inc., New York, 1939), by J. C. Slater, and Modern Theory of Solids 
(McGraw-Hill Book Company, Inc., New York, 1940), by F. Seitz. 
> v5 calculated from second virial coefficient. Data obtained from Introduction to Chemical Physics (McGraw-Hill Book Company, Inc., New York, 
1939), by J. C. Slater for Indene and Vinyl Acetate rs calculated from molecular volume. 
¢ Data obtained from Dipole Moments (Methuen, London, 1948), second edition, by R. J. W. LeFevre except in case of Indene and Vinyl! Acetate 
which was obtained J. Soc. Chem. Ind. (Japan) 43, Suppl. binding 190-1 (1940), by S. Lee. 
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Since by definition the static polarizability of a 
metal is its volume divided by 27, the average 


interaction energy for the element of metal con- 
sidered is 


Vom = — pu 2dv/2nr’. (19) 


Therefore, the expression for the interaction 
energy between the molecules and the entire 
metal becomes upon integration, over infinity in 


the x and y directions and from r to ~ along the z 
direction : 


. uw ¢ dv pe? 
Von= — —=S 


3 (20) 
2rJ, 1273 


DISCUSSION AND RESULTS 


From a practical standpoint the closest dis- 
tance of approach to the surface will be one of the 
prime variables in determining the absolute mag- 
nitude of the interaction energy. This distance is, 
of course, governed by the relative magnitudes of 
the attractive and repulsive forces. Since the 
latter are at present still beyond quantitative 
theoretical treatment, '*~"* it is necessary to resort 
to an empirical method of estimation. There 
appears to be little choice but to adopt a hard- 
sphere model and to take the distance of closest 
approach as the sum of the radii of molecule and 
metal atom. With this assumption several metal- 
dipole interactions have been calculated (Table I). 
The effects of temperature, radii, and dipole 
moment on the interaction energy are readily 
seen. 

If Eq. (8) is examined for the condition T=0 
and T= «, it is found that for T=0 


V.= — 2p?/R?, 
and for T= 


= 4 py? 
V 


3 R? 


so that in the lower limit the Boltzmann sta- 
tistics reduce to the required simple expressions 
(5) and (6). 

For the case of dipoles the total interaction 


4H. Margenau and V. W. Meyers, Phys. Rev. 66, 307 
(1944). 


1 R. Heller, J. Chem. Phys. 9, 154 (1941). 
16 W. H. Stockmeyer, J. Chem. Phys. 9, 398 (1941). 
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energy is 


2 “ i 
Oe 


i=1 Ki; 


(1+-cos?6;) 


nm 1 
+ > — [2 cos6,; cos6;; 
a Rat 5] 3 


—sin6,; sin6;; cos(¢.;— ¢3i) | ’ (21) 


where 7 represents the dipole and j represents the 
image of a different dipole. If the orientations are 
assumed to be random with respect to one another 
the second term of this equation may be neglected 
for in that case one obtains a repulsive force as 
often as an attractive force. This then leads to a 
simple extension of Eq. (3): 


n 1 
Via=—p? > —(1+c0s76;). 


(22) 
i=1 R,8 
While it is questionable whether one is justified 
in assuming complete randomness for a many- 
molecule problem, the net attractive force that 
might be expected as a result of dipole-metal 
interaction can be calculated, provided one is 
willing to accept complete randomness. By 
differentiating Eqs. (8) and (20) with respect to 
R and r the attractive force obtained is 
dV. dV. 
c= = , (Boltzmann statistics) 
dR 2dr 


dV om 
dr — 


(23) 


qm = 


(Quantum mechanical) (24) 


Here F represents the average force exerted by 
each dipole on the metal. To obtain the sum of all 
these forces over an element of unit area, we may 
integrate over all the 7’s since in a statistical 
treatment it is permissible to integrate over a 
continuous interval because the number of 
dipoles is sufficiently large to justify the assump- 
tion of continuity. The total attractive force 
per unit area can then be written as follows: 


Classical: 


oN 
5.=— 
MJ, 


a EP (@)—Ple)}, (25) 
dr =a m ~— 
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300 





Dipole 
force> 


(kg/cm?) 
A 


5380 


200 4280 

HO 150 9900 
300 10,240 

C.H,OH 150 1500 
300 1350 

CH;:COOCH:CH:2* 150 302 
300 286 

Pb NHs 150 3080 
200 2500 

H,O0 150 5740 
300 5780 

C.H,OH 150 926 
300 856 

CH;COOCH:CH:2* 150 218 


210 











TABLE IJ. Adhesion forces for dipole interaction. 


Dipole Total Total 
force® Dispersion adhesion adhesion 
(kg/cm?) force (classical) (Q.M.) 
B (kg/cm?) A B 
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3110 8490 5110 
1805 2785 7065 4590 
3560 2530 12,430 6090 
3960 2810 13,050 6770 
595 1470 2970 2065 
595 1470 2820 2065 
136 
136 
1260 1860 4940 3120 
1120 1670 4170 2790 
2150 1470 7210 3620 
2440 1630 7410 4070 
396 930 1856 1326 
396 930 1786 1326 
101 
100 











b Symbol A represents classical result. _ 
¢ Symbol B represents quantum-mechanical result. 


where §, is the attractive force per unit area, a is 
the minimum distance between metal and mole- 
cule, ¢ is the molecular density, M is the molecular 
weight, and N is Avogadro’s number. Then: 


26 


F.=7.54X 1022 





Ma?’ 


1 e* | 
M1115] 1 


2r 1 
f dx 
0 


Quantum Mechanical: 








woN ¢’ dr poN 
Fom = f = ’ (27) 
4M Y¥, r* 12Ma’ 
or: 
Foam = 5.02 XK 10?2(u2a/ Ma’). (28) 


In Table II the attractive or adhesion forces 
for dipole interaction of some of the molecules 
used in Table I are tabulated. Since the necessary 
constants for ammonia, water, and alcohol are 
known, the dispersion force could be roughly 
estimated and was therefore included in the 
computation of the total estimated adhesion 
force. The total attractive force per unit area 
from the dispersion effect is 


TQ, IIe 
2 


Ma® I;+I2 








Fa=7.50 X10? 


(29) 





® Insufficient data to compute dispersion effect for this molecule-metal interaction. 









where J; and J:=ionization energies for the 
molecule and metal, respectively, and a, and 
a@2=static polarizability of the molecule. 

It is not difficult to understand why at one 
time the van der Waals forces had been ascribed 
solely to the polarity of the molecules when one 
considers the contributions of the dipoles in such 
polar molecules as water and ammonia. However, 
today the contribution of the dispersion forces to 
the total van der Waals force is well known for 
many compounds. This effect is also observed 
very strongly in adhesion studies. For example, 
Kraus and Manson” measured the adhesion of 
polystyrene and polyethylene to metal by ex- 
trapolating to zero thickness the tensile test 
results and obtained adhesion values of the order 
of 500 kg/cm?. Czyzak!$ in calculating the adhe- 
sion of polystyrene to metal by considering the 
effect of permanent dipoles alone (u=0.26D) 
obtained an adhesion value of the order of 1 
kg/cm?*. If we use the latest value of the dipole 
moment of polystyrene (u~0) then the perma- 
nent dipole contribution would obviously be zero. 
This clearly indicates the importance of the other 
effects on the interaction energy and it shows that 
in the case of molecules with negligible perma- 
nent dipole moment a negligible adhesion value 


should certainly not be expected. Kraus and 
17G. Kraus and J. E. Manson, J. Polymer Sci. 6, 625 
(1951). 


18S. J. Czyzak, D. Sc. thesis, University of Cincinnati 
(1948). 
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Manson in discussing the theory of specific adhe- 
sion (total molecular attraction per unit area of 
interface) assume the attractive potential to be 
due to van der Waals forces alone. For a hypo- 
thetical molecule of »=2D and molecular weight 
of 100 they obtain 500 cal/mole for the contribu- 
tion of the permanent dipoles to the total 
interaction energy, whereas the contribution from 
all effects was estimated to be of the order of 8000 
cal/mole. 

Any experimental adhesion force measurements 
are in general sufficiently low that they can be 
accounted for by dispersion forces alone since 
with nonhydrogen bonding dipoles the dispersion 
effect is approximately ninety percent of all the 
effects. G. Kraus,'® using his previous method,” 
obtained an experimental adhesion value of the 
order of 600 kg/cm? for polyvinylacetate while 
the calculated data in Table II for dipole-metal 
(Fe) interaction gives an adhesion force of the 
order of 100 kg/cm? for this compound. Here the 
dipole moment is considerable (u=1.71D) in 
comparison to polystyrene (u~0.0D), and yet the 
dispersion effect clearly outweighs the dipole 
contribution. 

An experimental corroboration of theoretically 
calculated adhesion forces is not possible for the 


19G. Kraus, private communication, University of 
Cincinnati. 
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adhesion force must be measured as the force 
necessary to produce rupture. Rupture phe- 
nomena, in general, are however mainly de- 
pendent on nondominant structural features? 
rather than on the mean molecular attraction. 
The exact relation between molecular forces and 
adhesive rupture is as yet not known; experi- 
mental adhesion forces are, however, always 
much smaller than the total theoretical van der 
Waals attraction.'*® 

A recent contribution to adhesion testing is the 
electromagnetic ultracentrifuge developed by 
Soller, Liang and others including the author. !*-2! 
In this instrument adhesive rupture is produced 
under an extreme gravitational field. The method 
is especially useful in the study of the adhesion of 
coatings. 
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The really remarkable thing about the result stated above is that so interesting and significant 
a conclusion can be reached from such simple observations as the darkness of the night sky when 
compared to the surface brightness of a star, together with the cosmological principle. . 

The reason for the cosmological significance of such a simple fact as the darkness of the night 
sky is that this is one of the phenomena that depend critically on circumstances very far away. 
In laboratory experiments we always deal with quantities that we can influence, so as to examine 
the dependence of one measurement on another. Thereby we automatically exclude phenomena 
that depend only, or at least largely, on very distant matter (which is out of our control). In cosmology 
these are just the observations from which we can derive most information. Difficulty arises only 
because owing to their constancy and uniformity these effects of long-range influences do not attract 
our attention very much. We tend to ascribe an obvious and “‘absolute’’ character to such effects 
(when we notice them at all). The consistent interpretation of such effects forms one of the most 
important and significant parts of cosmology—Cosmology by H. Bonp1 (Cambridge University 
Press, 1952), p. 24, by permission. 
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Using reliable data on thermal conductivities K, viscosities 7, and specific heats cy, values of 
the important ratio f= K/nc, have been determined over a wide range of temperature, 80-380°K, 


for eight gases. For some gases the ratio shows a decrease with increasing temperature and for 
others an increase. The results and existing theories do not agree on all points. For a monatomic 
gas, helium, f has a constant value up to 250°K, agreeing with most theories, but it drops at 
higher temperatures. One diatomic gas, hydrogen, shows a sharp drop, agreeing very well with 
Eucken’s theory. Another diatomic gas, nitric oxide, shows a rise much sharper than predicted 
by the theory. Other diatomic and polyatomic gases show large deviations from Eucken’s 
theory, mainly in that experimental values rise while the theory predicts a decrease with in- 


creasing temperature. 


INTRODUCTION 


N the kinetic theory of gases there is an im- 

portant dimensionless ratio between thermal 
conductivity K, viscosity », and specific heat at 
constant volume per gram ¢,. The ratio f= K/n¢, 
is much greater than unity mainly because 
thermal conduction depends so much more on 
velocities than does viscosity. The actual magni- 
tude depends on the type of molecule concerned, 
as well as temperature. There are at least three 
good reasons for determining and considering the 
values of f and how they vary with several 
parameters. First, is the insight one might gain 
into the subject of inter- and intramolecular 
transfer of energy. Then, of somewhat more 
direct use, is the possibility of determining the 
specific heat of a gas whose properties are such 
as to make direct measurement difficult, for 
example, too low a vapor pressure. Finally, 
there is the great practical use for the engineer 
in heat-exchange calculations. Here the Prandtl 
number Pr=nc,/K=c,/fc, (where c,=specific 
heat at constant pressure) is directly used. Thus, 
it was thought worthwhile to examine the varia- 
tion of f over a wide range of temperatures for 
a variety of gases. 

The value of f has had various theoretical 
treatments, the most complete being made by 
Chapman and Cowling.! They find that a gas 
with only translational molecular energy should 
have f=2.5 within 1 percent. The value may 


*This paper is based on work performed unde- Uni- 
versity of California contract with the AEC. 

1S. Chapman and T. G. Cowling, Mathematical Theory 
of Non-Uniform Gases (Cambridge University Press, 
Cambridge, 1939). 
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vary from 2.500 for molecules with inverse fifth 
power repulsion to 2.522 for hard, elastic spheres. 
Available data on the inert gases at 273°K 
showed f=2.5 within 3 percent. However, in the 
case of more complex gases the addition of 
internal energy causes a proportional rise in 
specific heat but not in heat conductivity, since 
the latter property depends upon the molecular 
velocity to a great extent. Thus the internal 
energy ought to be treated separately from the 
translational and, in general, the over-all f ought 
to be broken into several components, one for 
each type of energy. Thus, 


f= (fieet fitrt fuCut fele)/Co, (1) 


where ¢, 7, u, and e denote translational, rota- 
tional, vibrational, and electronic, respectively. 
A simplified approach to this method was made 
by Eucken.? 

Assuming all forms of internal energy to have 
a rate of propagation equal to the rate of mo- 
mentum propagation, he set f,=f.=f.=fi:=1, 
where the subscript 7 denotes internal. Since, on 
a molar basis, C;=$R cal/mole and C;=C,—3R, 
he obtained a simple expression for fg the 
Eucken function, 


fr=1+4.47/C,=(9y—5)/4, (2) 


where y=C,/C,. Eucken’s formula has been 
used a great deal, for its results have been in 
good agreement with certain experimental data,® 


2 A, Eucken, Physik. Z. 14, 324 (1913). 


3 See, for example, Chapman and Cowling (see reference 
1), Table 22, p. 241, or E. H. Kennard, Kinetic Theory of 
Gases (McGraw-Hill Book Company, Inc., New York, 
1938), p. 180. 
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TEMPERATURE (°K) 


Fic. 1. The ratio f=K/nc, as a function of 
temperature for several gases. 


which, unfortunately, have been limited mainly 
to 273°K. 


PRESENT CALCULATIONS 


Recently there have become available sufficient 
reliable data to give experimental f values over 
a wide temperature range and for a variety 
of molecules—monatomic, diatomic, and poly- 
atomic. For the present study the following 
eight gases were chosen: helium, hydrogen, 
oxygen, carbon monoxide, nitric oxide, carbon 
dioxide, nitrous oxide, and methane. The thermal 
conductivities had been measured over the range 
80°-380°K,* while the viscosities had been 
measured from 80° to 300°K®5 and were, for the 
sake of consistency,® extrapolated by the author 


4H. L. Johnston and E. R. Grilly, J. Chem. Phys. 14, 
233 (1946). 

5a. H. L. Johnston and K. E. ana J. Phys. 
. 44, 1038 (1940); Hz, O2, NO, COs, N20, 

b. H. L. Johnston and E. R. Grilly, J. Phys. Chem. 46, 
948 (1942); He, CO. 

Wobser and Fr. Miiller, Kolloid-Beihefte 52, 165 
(1941, have given Sutherland constants between 293° and 
373°K and viscosities at 293°. The latter agreed with those 
in reference 5 within a few tenths percent, except for NO, 
whose viscosity was lower by 0.8 percent. At 380° their 
data deviated from the extrapolated values used in the 
present computations by the following percentages: He, 
—2.0; Ho, —0.4; O., —0.4; co, —1.3; NO, —1.4; CO:, 
—0.7; N.O, +0.2; CH,, +0.6. 


to 380°K. The heat capacities used had been 
calculated from spectral data.’ In all these cases 
the heat capacities had been calculated for the 
gases in the ideal state. Essentially, the viscosity 
and thermal conductivity data used were also for 
the ideal state since the measurements had been 
made at low pressures. 

The values of f= K/ne, and of fe =1+4.47/C, 
are plotted as functions of temperature in Figs. 
1 and 2, respectively. The Prandtl numbers 
Pr=nc,/K are plotted in Fig. 3. The accuracy of 
f is limited by K and 7; if each has a maximum 
error of 0.5 percent, f should certainly be deter- 
mined to 1 percent. The fg values should be 
good to about 0.1 percent, as the C, values are 
accurate to at least 0.2 percent. For hydrogen f 


f,21+4.47/C, 


TEMPERATURE (°K) 


Fic. 2. The Eucken function fg=1+4.47/C, asa 
function of temperature for different gases. 


7a. W. F. Giauque, J. Am. Chem. Soc. 52, 4816 (1930): 
H:2 below 250°K. b. H. L. Johnston and C. O. Davis, J. Am. 
Chem. Soc. 56, 1045 (1934): He above 250°K. c. H. L. 
Johnston and M. K. Walker, J. Am. Chem. Soc. 55, 172 


(1933): Oz. d. H. L. Johnston and C. O. Davis, J. Am. 
Chem. Soc. 56, 271 (1934): CO. e. H. L. Johnston and A. T. 
Chapman, J. Am. Chem. Soc. 55, 153 (1933): NO. f. E. A. 
Long, Ph.D. dissertation, Ohio State University, 1934: CO 
and N,O. g. D. P. MacDougall, Phys. Rev. 38, 2296 (1931): 
CH, rotational contributions. h. R. D. Vold, ‘ Am. Chem. 
Soc. 57, 1192 (1935): CH, vibrational contributions above 
273°K. i. The CH, vibrational contributions below 273°K 
were calculated with the aid of Einstein functions, using 
the same frequencies as Vold. 












and fg agree within 1 percent over the entire 
range. For helium the agreement is good up to 
250°K, then f drops steadily until it is 7 percent 
below fz at 380°K.® 

Nitric oxide, the only gas considered whose 
specific heat decreases with increasing tempera- 
ture (due to a maximum at 75°K caused by an 
electronic transition) has both f and fg increasing 
with temperature, although at different rates. 
For all the other gases little correlation between 
f and fz is observed over much of the tempera- 
ture range studied, but fair agreement is ob- 
tained for all the gases near 300°K. 


DISCUSSION 


One of the most glaring faults of the simple 
Eucken formula is setting f;, the coefficient for 
internal energy, equal to unity. This indicates 
that the mean free path for transport of internal 
energy is taken to be the same as for transport 
of momentum. A more logical step involves 
taking the mean free path appropriate to self- 
diffusion.°® 

The result is f;=Dj,:d/n, where Dy, is the self- 
diffusion coefficient and d is the density. Since 
experimental D,, values are generally unavail- 
able, we must resort to an indirect method for 
testing this approach. Theoretically, all values 
of D,,d/n seem to be above unity. Chapman and 
Cowling’ give it as a function of m for molecules 
repelling as the inverse mth power of the distance: 
1.204 for m= and 1.543 for n=5. Amdur"™ 
gives equations which yield 1.06<D,;d/n<1.18 
for He, Ne, A, Kr, and Xe, and which also result 
in a slight positive temperature coefficient for 
Di,d/n. However, we can calculate f; empirically 
from a rearrangement of Eq. (1): 


fi=(fC,—7.45)/Ci. (1a) 


The results for the gases concerned here give the 
same pattern, an increase from about 0.7 to 1.3 
with increasing temperature. The values lower 
than unity are hard to explain. The higher ones 


8 Even if the viscosity given by Wobser and Miiiller® 
were accepted, the difference would still amount to 5 
percent. 

® Reference 1, pp. 238-240. 

10 Reference 1, p. 238. 


1], Amdur, J. Chem. Phys. 15, 482 (1947); 16, 190 
(1948). 
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Fic. 3. The Prandtl number Pr=nc,/K as a function 


of temperature for the 8 gases for which data are presented 
in Figs. 1 and 2. 


are compatible with the diffusion mechanism for 
propagation of internal energy. 

An interesting theoretical approach to the 
evaluation of the over-all f has been given by 
Chapman and Hainsworth.” Their molecular 
model was the rough sphere of variable radius, 
where “rough’”’ means nonslipping at the point 
of contact in the collision of two spheres; thus 
the relative velocities are reversed by the colli- 
sion, and the interconversion of translational 
and rotational energies occurs. Their results 
show that f depends on how much the molecular 
radius varies with the velocity of approach: the 
greater the variation (or softer the molecule), the 
smaller is f. Since there is no way to choose 
correct parameters except empirically, the 
theory’s value is mainly limited to determining 
the variation of f with temperature. The positive 
temperature coefficient produced is more in line 
with experimental results than is the usual nega- 
tive coefficient given by Eucken’s theory. 


12S. Chapman and W. Hainsworth, Phil. Mag. 48, 593 
(1924). 
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CONCLUSION 


Looking at these various approaches to inter- 
relating transport properties of gases, we must 
conclude that each has severe limitations. The 
simple Eucken relation is probably still the most 
useful. The disadvantage of its incorrect de- 


E.. R. GRIELY 


pendence on temperature might often be out- 
weighed by its nearly correct average result. To 
test some of the other hypotheses, we need data 
from two fields: (1) self-diffusion in gases and 
(2) exchange between translational and internal 
forms of energy among molecules. 


On the Physical Interpretation of Solutions of the Dirac Equation for a Free Particle 


S. M. NEAMTAN 
University of Manitoba, Winnipeg, Canada 


(Received December 11, 1951) 


The usual choice of an orthogonal set of four plane-wave solutions of the free-particle Dirac 
equation does not lend itself readily to direct and complete physical interpretation except in low 
energy approximation. A different choice of solutions can be made which yields a direct physical 
interpretation at all energies. Besides the separation of positive and negative energy states 
there is a further separation of states for which the spin is respectively parallel or antiparallel 


to the direction of the momentum vector. 


LANE-WAVE solutions of the Dirac equa- 
tion are well known, and it would appear 
to be necessary to justify a reconsideration of 
this problem. The justification lies in the fact 
that the solutions found in textbooks and else- 
where do not lend themselves readily to a direct 
physical interpretation. From a_ pedagogical 
point of view such an interpretation is much to 
be desired. In the usual treatment of the prob- 
lem,! four mutually orthogonal solutions are ob- 
tained for a given value of the momentum p, two 
of which correspond to positive energy and two 
to negative energy. However, the physical dis- 
tinction between a pair of solutions belonging to 
the same energy does not appear until one con- 
siders the low energy approximation in which 
the small components of the state function are 
neglected and the two solutions reduce to eigen- 
functions of o, belonging to the eigenvalues +1. 
In the usual notation, the Dirac Hamiltonian 
for a free particle is 


H=-—ca-p—mc’s. (1) 


This Hamiltonian commutes with the momentum 
vector p, and the usual procedure is to seek 
simultaneous eigenfunctions of H and p. These 
eigenfunctions are, however, not uniquely de- 


1See for example, L. I. Schiff, Quantum Mechanics 
(McGraw-Hill Book Company, Inc., New York, 1949). 


termined, and for given eigenvalues of H and p, 
there remains a twofold degeneracy. In order to 
resolve this degeneracy we seek a dynamical 
variable which commutes with both H and p. 
Such a variable is o-p. It is obvious that this 
variable commutes with p. To verify that it also 
commutes with H, we write 


a=pit= 90/1, 
and recalling that 8 commutes with o, we have 
o- pH —He-p= (cpis-p)(o-p) 
—(e-p)(cpio-p)=0, (2) 
since p,; commutes with oe. 

We now proceed to find simultaneous eigen- 
functions of the commuting variables H, p and 
o-p. We have, since the components of p 
commute, 

(o-p)*=?", (3) 


where p is the magnitude of the momentum 
vector. Thus for a simultaneous eigenstate of p 
and o-p, the value of o-p will be + or —?, 
corresponding to states for which the spin is 
parallel or antiparallel, respectively, to the 
momentum vector. 

A simultaneous eigenfunction of H and p will 
have the form of a plane wave 


vj;=u;exp[i(p-r—Et)/h], j=1,2,3,4, (4) 
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where the y; are the four components of the state 
function and u; four numbers to be determined. 
In the argument of the exponential function, p 
represents the eigenvalues of the components of 
the momentum for this state and E the corre- 
sponding eigenvalue of H. Then E can have 
either of the two values 


E=+e=+(m'ct+cp’)!. (5) 
We now demand that y; be also an eigenfunction 
of o-p belonging to one of the eigenvalues 7, say, 


where r= +p. Employing the usual matrix repre- 
sentation for ¢, we have 


bz pz—thy 
sil piettpy —?P 0 0 
0 0 be = Pz tpy | 
0 0 Ppettpy —DPs 
In the above matrix, p., py, and p, are operators, 
but since this matrix is to operate on an eigen- 
function of p, the operators can be replaced by 
their eigenvalues. We shall, without risk of con- 
fusion, use the same symbols for the eigenvalues 
as for the corresponding operators. 
The eigenvalue equation 


o-py=ry, (7) 


yields the following four equations: 


(6) 


Pelit (p2—tpy)u2=TM1, 
(petipy)ui—ple= Tue, 
Pits t (p2z—tpy)Us= THs, 
(petipy)Us— pPel4= TU. 


These equations are satisfied by 


Ui=pz—tpy, uU2=17— pz, 
U3=X(pzr—ipy), us= (4 — pz), 


where A is a number still to be determined. (The 
consistency of the equations is readily verified, 
recalling that r= +p.) 

We now make use of the requirement that y 
be an eigenfunction of H belonging to the eigen- 
value E. Then 


(9) 


(E—H)y=0; (10) 


that is, 


(E+ cpio: p+mcs)y =0. (11) 


In this Eq. (11) we replace o-p by its eigenvalue 
to obtain 


(E+cmrp,+mc’p)y =0. 


On expansion, we obtain 


(12) 


0=(E+mc)u;+crusg, 
0=(E+mc*)uetcrus, 
0=(E—me*)uzt+cruy, 
0=(E—me*)uy+crue. 


(13) 


The condition E=-+te ensures the consistency 
of these equations. Comparing Eq. (13) with 
Eq. (9), we find 


E+mc cr 
— =— . (14) 


cr E-—mc 


=> 


We can now write the components of yw as 
follows: 


¥j;=u; exp[t(p-r—Et)/h], 


where the uw; are given by 


(15) 


U1=pz—tpy, Uu2=1— pz, 
us= —(E+mce*)(p.—ip,)/cr, 
ug= —(E+mce*)(x—p.)/er. 


Since E can be given either of the two values 
+e and 7, the two values +), we have found for 
given p four linearly independent plane wave 
solutions.? It is easily verified that they are 
mutually orthogonal. 

The physical interpretation of the solutions is 
now clear. Each solution represents a homo- 
geneous beam of particles of definite momentum 
p, of definite energy, either +e, and with the 
spin polarized either parallel or antiparallel to 
the direction of propagation. 


(16) 


2 It will be noted that the solution (16) appears to break 
down for the special case 
pbz=py=0, w=pz. 
It is easily verified that for this case the solution can be 
written 
u,=2p, u2=0 (17) 
ug=—2(E+me*)/c uy=0. 
This form corresponds to that to which Eq. (16) reduces 
for the case pz = py=0, r= — pz, viz: 
uy =0 u2= —2p; 


us=0 ug= —2(E+ me?) /c. (18) 
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The Multiple-Slit Diffraction Pattern 


D. H. RANK 
The Pennsylvania State College, State College, Pennsylvania 


EXTBOOKS on physical optics almost universally 

discuss the characteristics of the multiple-slit dif- 
fraction pattern and then show that this theory can be 
applied directly to the diffraction grating. One of the main 
features of the diffraction pattern is that the secondary 
maxima decrease in intensity extremely rapidly, and that 
there are N—2 of these secondary maxima between orders 
of the grating, where N is the number of slits. The case 
that is most universally treated in theory is that in which a 
is the width of any slit and d is the constant separation of 
adjacent slits. Under these conditions the intensity distri- 
bution J for normal incidence will be given by the following 
expression! 

nae = 
sin*y 

where B=(masin@)/A and y=(zdsin@)/A, @ being the 
angle of diffraction and \ the wavelength. Although any 
real diffraction grating does not quite fulfill the above 
theoretical conditions it is to be expected that the above 
expression will describe the qualitative features of the 
diffraction pattern exactly, and the intensity distribution 
will be modified only in a very minor way. 

To the best of my knowledge there are no examples in 
the literature showing the fine detail of the diffraction 
pattern produced by a diffraction grating. Some years ago, 
when I was testing the figure of a large plane grating 
illuminated by the ‘“‘white”’ light of a mercury arc, I noticed 
a maze of extremely sharp weak fringes in the neighbor- 
hood of the central image. Making use of monochromatic 
light, e.g., 45461 of Hg, it was immediately apparent that 
these fine sharp fringes were the maxima of the multiple- 
slit diffraction pattern produced by the reflection of light 
from the rulings. I have photographed this diffraction 
pattern in the following manner. 

The large plane grating was mounted between colli- 
mating and focusing mirrors each of ten meters focal 
length. The grating was set so that the central image 
produced by the rulings fell in the center of the field of 
the ten-meter camera. Illumination was by means of the 
Hg green line isolated by means of a suitable filter. In order 
to make the detail of the diffraction pattern large compared 
to the grain size of the photographic plate a strip of rulings 
only one-half inch in width, that is, 7500 lines, was used 
in making the photographs. Under these conditions a 
beautiful set of fringes could easily be observed with an 
eyepiece. Schuster? has shown that when a ‘“‘normal”’ slit 
width is used one can expect the observed diffraction 
pattern to be practically indistinguishable from that pro- 
duced by an infinitely narrow slit. The ‘“normal’’ slit 
width is defined as (A/4) F, where F is the ratio of the focal 
length of the collimator to the apparent width of the 
grating viewed from the position of the slit. In the present 


Fic. 1. Enlargements of diffraction patterns produced by a reflection 
grating of 7500 lines. The diffraction pattern was produced by the zero 
order of the grating. The slit width used was the “‘normal” slit width, 
i.e., 0.11 mm. Illumination by 5461 of Hg. A, B, and C are progressively 
heavier exposures made to show detail of the diffraction pattern. The 
detail in the center of the pattern of 1C was brought out by “dodging”’ 
the edges during enlargement, i.e., shading the exposure of the edges 
of the pattern to the printing paper by moving one’s hands sym- 
metrically about the center in such fashion so as to achieve the appear- 
ance of all the detail on the original negative on the positive print. 


case this normal slit width was 0.11 mm. It was noted 
that the diffraction pattern deteriorated rapidly when the 
slit width was increased to several times the normal value. 

In Fig. 1 enlargements of some of the plates obtained are 
reproduced. Figures 1A, B, and C are progressively heavier 
exposures made to show the detail of the multiple-slit 
diffraction pattern. The highest order fringe reproduced in 
Fig. 1C is many thousands of times weaker than the central 
fringe. I have not been able to observe this diffraction 
pattern in the spectra other than the spectrum of zero 
order because of complications produced by various and 
sundry ghosts present even in the best diffraction grating 
spectra. 

1Jenkins and White, Fundamentals of Optics (McGraw-Hill Book 


Company, Inc., New York, 1950), p. 323. 
2 Arthur Schuster, Astrophys. J. 21, 197 (1905). 


A New Vector Identity and Physical Application 


EARLE B. MULLEN* 
University of Scranton, Scranton, Pennsylvania 


HE vector identity for grad(A-B) that one commonly 
finds in textbooks! is 


V(A:B) =AX(VXB)+BX(VXA)+(A-V)B+(B-V)A. 


But another identity of equal validity which a cursory 
search of the literature has not revealed is 


V(A-B) =(AXV) XB+(BXV) XA+A(V-B)+B(V-A). (1) 


This equation can be verified as follows: The ordinary 
vector triple product 


(AX C)XB=C(A-B) —A(C-B) 
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leads us to suspect that 


(AXV) XB=Va(A-B)—A(V-B), (2) 


where the subscript A means that the vector A is being 
held constant. That Eq. (2) is actually true can be verified 
by expanding both sides: 

t j k 
Az Ay As 
d/ax a/ay a/az 


(AXV)XB= xB 








a j 
=|A,a/az—A,0/dy, A~/dx—A,0/d2, 
B, By 


k 
A,0/dy—A,0/dx}, 
Bz 


the x component of which is 
AB,/dx—A,0B,/02—A,0B,/dy+A,0B,/dx. (3) 


The x component of the right-hand side of Eq. (2) is 
given by 


[va(A-B) —A(V-B)]. 
= (a/ax)4(AzB,+A,B,+A.B,) 
—A,(dB,/ax+0B,/dy+0B,/d2) 
=(A,0B,/dx+A,0B,/ax+A,dB,/ax 
—A,0B,/ax—A,0B,/ay—A,0B,/d2). (4) 


When the first and fourth terms are canceled, Eq. (4) is 
seen to be identical with Eq. (3) thus verifying Eq. (2). 
Interchanging A and B in Eq. (2) gives 


(BXV)XA=VB(A-B)—B(V-A). (5) 
Adding Eqs. (2) and (5) we find that 


(AX V) XB+(BXV) XA=Va(A-B) 
+VsB(A-B) —A(V-B)—B(V-A). 


But the first two terms on the right combine to give 
V(A-B). Transposing the other two terms, we obtain 
Eq. (1). 

An application of this identity is a direct integral proof, 
without making use of the concept of the vector potential, 
of the theorem that a closed circuit carrying a current is 
equivalent to a magnetic shell of a certain strength. If m is 
the strength, i.e., the magnetic dipole moment per unit 
area of the shell, then the potential at an external point P 
is given by 


V= - (m-1)r-%dS = — J. m-V(1/r)d5, (6) 


where r is the vector from an element of area dS of the 
shell S to the field point P. The magnetic field at P due to 
the shell is 


ee ae J V(m-Vs(1/r)d5S, (7) 


where Vs indicates that the differentiation is done with 
respect to the coordinates of the surface S and use has 
been made of the fact that V= —Vsz. 





DISCUSSION 






If we apply Eq. (1), we can write B in the form 
B= -f {—m(V-V(1/r) +Vs(1/r)(V-m) 


+(mXV) XVs(1/r) +(Vs(1/r) XV) Xmjd5S. 


The first term in the integral vanishes since V?(1/r) =0. 
The second and fourth terms are zero since m is inde- 
pendent of the coordinates at P. Therefore we are left with 


in J (mXV) XVs(1/r)dS. (8) 


Since mdS=mdS and m is a constant, this can be ex- 
pressed as 


Se J (dSXVs) XVs(m/r). (9) 


Now using a corollary of Stokes’s theorem (see Lass, 
reference 1, p. 112) which states that 


ee (@SxVs)xA=f, dlXA, 
where C is the bounding rim of S, Eq. (9) becomes 
B= if dlXVs(m/r) = I (mdl x r)r-3, 


which, if m is taken to be a certain constant times the 
current, is the Biot-Savart law for the magnetic field due 
to a closed electrical circuit. 


(10) 


* Nowat The Electronics Laboratury, The General Electric Company, 
Electronics Park, Syracuse, New York. 

1E, q Lass, Vector and on _— (McGraw-Hill Book 
Company, Inc., New York, 1950), p. 





A Simple Class Experiment on Plastic Flow 


H. Mykura 
University College of the West Indies, Jamaica 


LTHOUGH most laboratory courses in physics in- 

clude a number of experiments on the determination 
of elastic moduli, there are usually none on plastic flow 
and fracture. 

The reason for this is that the usual type of clamp intro- 
duces such stress concentrations in the specimen at the 
point of clamping that fracture occurs there and the stress 
cannot be measured; so the specimen must be specially 
shaped (as in the engineering tensile test), which makes it 
unsuitable for a physics teaching laboratory. 

A simple clamp will support a wire specimen with suffi- 
cient absence of stress concentration to enable measure- 
ments of the elongation to be made, in loading it to frac- 
ture. The clamp consists essentially of a wooden cylinder, 
about 1 in. in diameter; the wire is wrapped round the 
cylinder for three-quarters of a turn and then clamped to 
the cylinder by two screws and a piece of brass. Friction 
between the wire and cylinder supports most of the tension 
in the wire, so that stresses at the clamp are insufficient to 
deform the wire plastically during loading. For convenience 
in changing the specimen a second brass piece is fitted to 
the cylinder to hold the clamping piece. The other end of 
the wire is fixed to a similar clamp which has a rack for 
slotted weights (or scale pan) attached. 

Figure 1 shows the apparatus in use. Two paper flags 
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Fic. 1. The appa- 
ratus in use. A, speci- 
men; B, paper flags to 
define gauge length; C, 
meter rule; D, cylindri- 
cal clamps. 













are fixed to the specimen to mark off a “gauge length.” 
A meter rule fixed behind the wire enables the length to 
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On Demonstrating a Classical Problem 
in Analytical Mechanics 


HE usual course in analytical mechanics includes no 
laboratory work and very little, if any, demon- 
stration. The problems encountered are, as I choose to 
call them, ‘‘fictional.”” They have application, to be sure, 
but the applications are pretty obscure to beginners, for 
they are usually found in intricate and involved engineering 
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be measured to 0.5 mm quite easily. The wire is loaded, 
the length being measured at each load, till fracture occurs. 
The specimen will not always break in the gauge length; 
but even if it breaks near the clamp the necking at the 
point of fracture is quite symmetrical. 

From the results a load-extension graph can be drawn. 
Measurement of the wire diameter and length enables one 
to convert this into a stress-strain curve. The diameter of 
the broken specimen at the point of fracture can be 
measured by traveling microscope, so the maximum stress 
at the point of fracture can be calculated. 


Q 


























2 ye é 6 
ExTENSION (cm ) 


Fic. 2. Load-Extension curves. Copper specimens 35 cm long, 0.053 cm 
diameter initially. A, as received; B, annealed. 


F Soft iron or copper wires of 22 to 26 S.W.G. are suitable 
materials for this experiment. The wire, as supplied, is 
usually in a work-hardened condition. The copper wire 
can be easily annealed by heating in a Bunsen flame, so 
the tensile properties of the work-hardened and annealed 
material can be compared. With copper an appreciable 
amount of creep occurs, so the results depend to some 
extent on the rate of loading. 

The graphs (Fig. 2) show load-extension curves for 
copper wire in the ‘‘as received”’ and ‘‘annealed”’ condition. 
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devices and machines. It seems appropriate, therefore, to 
demonstrate our findings in the purely ideal problems, if the 
physical set-up can be approximated. In addition, results 
are often come to analytically, that is, mathematically, 
which are not easy to see. They are true enough, we trust, 
since who will gainsay the rigorous analysis, but it is often 
difficult for the beginner to get the feel of the thing. I have 
often come upon a case where I believed the result to be 
true but have in the next breath said “I wish I could see 
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that.’’ Toward this end I have been trying to demonstrate 
certain classical problems in mechanics. Many, of course, 
assume ideal conditions, and these, more often than not, 
cannot be even approximated, let alone realized. The one 
I report here, however, far exceeded my expectations, 
and it does exactly what I want it to do. 

The problem.—A uniform rod of mass m and length / is 
free to rotate about a horizontal axis through one end. 
It is supported horizontally by a vertical string at the 
other end. The string is cut. Find the reaction of the axis 
at every instant of its fall. The analysis shows that the 
vertical hinge reaction is always 5/2 mg cos@, and the 
horizontal hinge reaction is given by 1/4 mg sin@, where 6 
is the angle with the vertical. It follows that: 


(1) The reaction is independent of the length of the rod. 

(2) At the instant of release, when the rod is horizontal, 
the hinge reaction is mg/4. 

(3) When the rod passes through the vertical, the hinge 
reaction is 5/2 mg. 


Obviously, when the rod is he/d in a horizontal position 
by supporting the free end, the hinge reaction is mg/2. 
At the instant this support is removed the hinge reaction 
goes to mg/4. This is very interesting! Can we show these 
things? 

Demonstration—A meter stick with a hole on one end 
is supported by the hook of a dial-type spring balance. 
In the vertical position at rest the scale reads (in my 
experiment) 150 grams; when held horizontally the scale 
reads 75 grams. At the instant the supporting string is cut 
(done best by burning it through) the needle flicks back to 
about 40 grams. As the stick passes through the vertical 
the scale reads very nearly 370 grams. Repeated trials 
yield a pretty good confirmation, after skill is acquired in 
catching the reading. 

Now cut the meter stick in half and fasten the halves 
together along their length. Elastic bands do it well. We 
have the same mass as before but half the length. Shall we 
expect the demonstration now to confirm our analytical 
results? It should. A logical extension is to cut the halves 
into quarters, which leads us to ask about the limiting case. 

I discover the demonstration to lead to much more 
penetrating inquiry than any blackboard analysis provides. 


Jutius SUMNER MILLER 
Dillard University 
New Orleans 22, Louisiana 





Discussion of Units and Dimensions in Physics 


AY I add something to the discussion of the wel- 

come article by C. H. Page! on this subject which 
appeared recently in this journal. I say that the article is 
welcome because even in our texts of best repute the 
author will occasionally yield to the temptation to con- 
sider the equations expressing physical laws as being 
numerical equations. Such a numerical equation contains 
our old friend k, which has a certain value corresponding 
to each choice of unit, and which may be made unity by 
the proper choice of units. 
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The equation most mistreated in this manner is probably 
the one chosen by Page, namely, the statement of Newton’s 
second law, F=ma. For example, in that excellent text by 
Sears? equations are generally treated as physical equa- 
tions, but F=ma is, in the beginning of the discussion, at 
least, treated as a numerical equation. 

A possible reason for this might be found if we examine 
the procedure by which we obtain the formulation of 
Newton’s second law. We usually describe experiments 
whose results are expressed by the proportion, F «ma, 
where for the purpose of measuring F in this experiment 
we must use some preconceived system of units, which we 
are at liberty to discard later if we wish, once their func- 
tion in obtaining the law has been fulfilled. If, now, we are 
using a system of units based on units of mass, length, 
and time, then the units of m and a are determined, but 
the units of F remain to be defined. The question arises, 
how are we to define F? The answer is that we define F by 
the equation F EQUALS ma. Thus, if we are using the mks 
system, our force will have the units kg-m/sec?, a composite 
unit which we have chosen, purely for the sake of con- 
venience, to call the newton. Similarly, in the cgs system 
we get force in g-cm/sec?, a composite unit which we 
call the dyne. If for some reason we should choose to 
measure distances in feet, and mass in kilograms, we would 
express our force in kg-ft/sec?, a perfectly respectable, but 
unconventional force which no one has bothered to name, 
since it would seldom be used. If we do come out with 
such a unit for a force we hasten to convert the result to 
newtons, dynes, or pounds, which everyone understands. 

The fundamental objection to & in this case is that it 
conceals the fact the force 7s mass times acceleration. As an 
example of a place where the use of & is appropriate, 
consider the equation F=kx, relating, say, the tensile 
force on an elastic bar to the resulting extension. Here 
force is not extension, but its magnitude is proportional to 
the magnitude of the extension (F«x), and this propor- 
tionality is made an equality by inserting the constant k, 
giving it the appropriate value and units. Thus, to write 
F=kma is to say that force is mot mass times acceleration, 
whereas it actually is. 

Having made this strong statement, may I retreat a bit 
and admit that a case might be made for the use of & in 
the example which I am about to give. Suppose that we 
had arrived at the proportional formulation of Newton's 
second law with a unit of force (say, the ugh) already 
defined, in addition to the units of mass, length, and time, 
and the ugh did not happen to have been chosen just right. 
Then indeed, we might write F=kma, where k would have 
such units as ugh-sec?/kg-m, or, ugh/newton. 

There are many places where this question of k does not 
enter. For example, momentum is defined by the equation, 
M=mv. Typical units for momentum are g-cm/sec. 
Since no one has thought it worthwhile to coin a name for 
this set of units, the question of whether or not to use k 
does not arise. 

A further step in the direction of the numerical equation 
is taken when a statement is made such as ‘‘In other words, 
the weight of a body, when expressed in terms of the force 
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unit of any system, is numerically equal to the mass of the 
body, in the mass unit of that system, multiplied by the 
corresponding value of the acceleration of gravity’? (italics 
mine). This is an example of the explanation of an equation, 
an explanation which should be unnecessary. The use of 
the explanation hides from the student the fundamental 
nature of the equation, and gives him the idea that he 
has to memorize a lot of arbitrary ‘‘formulas’”’ to get along. 
If he understands the equation he can figure out the con- 
sequences. 

The method of conversion of units used by Page in his 
article is similar to the one used by Worthing,‘ based on 
the idea that all conversion factors are unity. This well- 
known idea, applied to the example in the article, would 
give, 

mi 1 hr 


mi 
=60— =60— x ————_- 
" - 3600 sec 


5280 ft hi ft 
hr hr 7 : 


88 


1 mi sec 


In this connection may I remark that the reference to an 
article by Worthing® leaves the impression that Worthing 
was one of those who regarded the equations as numerical 
ones. Actually, he was a vigorous proponent of the idea of 
physical equations. 

In a final section of his article Page has discussed the 
question of whether to use three or four fundamental 
dimensions—whether /, m, and ¢ are sufficient, or whether 
some electrical dimension, such as g, should be added to 
the above. His discussion of this is rather long, and seems 
to the present writer to miss the fundamental point of 
distinction, which is a philosophical one. This point comes 
up regularly at the beginning of a course, Electricity and 
Magnetism, taught by the writer, when the different 
systems of units are discussed in connection with Coulomb’s 
law for the attraction of two electrical charges. The 
writer presents to the class the following viewpoint: the 
fundamental distinction between the cgs and the mks 
systems is that the former is based on three dimensions, 
m, 1, and t, whereas the latter has a fourth dimension, q. 
Do you believe that an electrical phenomenon can be 
expressed by an equation containing only mass, length, 
and time? If you do, there is no philosophical reason for 
preferring one system over the other. If you believe, on 
the other hand (as the writer does), that there is something 
different about electricity which must be expressed by an 
additional electrical dimension, then the mks system must 
be your choice. Any difference of magnitudes between units 
in the two systems is incidental, and the choice of the 
particular electrical dimension is, of course, arbitrary 
The writer feels that this is the basis on which the dis- 
cussion of the fourth dimension should be based. 


; Francis T. WORRELL 
Rensselaer Polytechnic Institute 


Troy, New Yor. 


; (mae Page, ‘Units and dimensions in physics,’"” Am. J. Phys. 20, 
2F. W. Sears, Principles of Physics (Addison-Wesley Press, Inc., 
Cambridge, Massachusetts, 1947), pp. 75, 76. 
3 Reference 1, page 79. 
4A. G. Worthing, Am. J. Phys. 8, 199 (1940). 
5 A. G. Worthing, Am. J. Phys. 14, 359 (1946). 
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AGREE with Professor Worrell that electrical phenom- 

ena are different from purely mechanical phenomena, 
or at least are conveniently separated in our present 
theories. It is difficult to understand how a logical proof 
of the necessity of another dimension can follow, for our 
present knowledge of electromagnetism was developed in 
the m, 1, system. A comparable situation arises in thermo- 
dynamics. The subject is usually taught with the conven- 
tion that temperature is dimensionless, but assigning a 
new dimension to TJ certainly simplifies dimensional 
analysis. 

It may be argued that assigning additional dimensions 
for reasons of convenience only is arbitrary. My viewpoint 
is that amy use of dimensions in physical theories is purely 
for convenience; the criterion of consistent behavior there- 
fore calls for using the most convenient set of dimensions. 


CHESTER H. PAGE 
15400 Norwood Pike, Rt. 1 


Silver Spring, Maryland 


A Demonstration on Relationships between 
ac and dc Voltage 


HE equality of the heating effects of ac and dc power 
when applied across equal resistances can be used to 
define the equality of ac and dc voltages. This follows from 
the definition of the ac voltage in terms of the root mean 
square of the variable potential applied across the re- 
sistance. 

However, the root mean square concept is often not 
grasped by elementary students, whereas a demonstration 
of the equality of the heating effects of ac and dc presents 
tangible evidence. The ordinary methods of calorimetry 
are too involved for such a demonstration. We describe 
below a photometric technique for comparison of heating 
effects simple enough to be used as a classroom demonstra- 
tion. Before describing it, we mention a qualitative demon- 
stration designed to arouse the 8 o’clock student who has 
come to class without breakfast. Two hot dogs are heated, 
one on ac and one on dc, by inserting nails as electrodes 
into their ends and sending the ac or dc through them. 
Very quickly they show visible signs of heat. They both 
smell alike and cause similar salivary disturbances. They 
can supply physical, if not mental food, but the demonstra- 
tion lacks quantitative precision. 

To supplement the above with intellectual meat we 
perform the following demonstration. A paraffin photom- 
eter! is used to compare the luminosity of two lamps. 
The student may be cautioned against thinking that a 
greater luminosity always implies a greater heating effect. 
In this experiment, however, we can validly assume that 
equal luminosities imply equal heating effects. 

Figure 1 shows the essential parts and connections. Two 
similar 40-watt lamps have been found adequate. If both 
lamps are first connected to the same variable ac source 
and the photometer is placed so that photometric balance 
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RHEOSTAT AUTOTRANSFCRMER 


PARAFFIN BLOCK 


Fic. 1. The equality of the heating effects of similar ac and dc voltage 
may be judged photometrically and used to calibrate an ac voltmeter. 
Since the two lamps will have the same color when balanced, color 
differences will not introduce error. 


is achieved, it may be left there regardless of geometrical 
asymmetry. To run the demonstration, one lamp is con- 
nected through a rheostat to dc and the other to an auto- 
transformer and ac as shown. In practice, the dc voltage 
across lamp No. 1 is set arbitrarily by adjusting the 
rheostat. Then, looking at the photometer but not at the 
ac voltmeter, the autotransformer is turned until lamp 
No. 2 is just as bright as lamp No. 1 as judged photo- 
metrically. The reading on the ac meter may then be com- 
pared with that of the meter. We have found that it is 
easy to get agreement within 3 percent. 

The above is adequate for an experiment, but to exploit 
the dramatic possibilities, a voltmeter whose readings can 
be projected is very effective. For example, let us assume 
that the ac meter readings can be projected without 
seriously affecting the photometer balance. A good pro- 
cedure is to have the instructor read and record the dc 
voltage. He then announces to the class that he will set 
the ac meter at the same value, looking at the photometer. 
The class watches intently to see how closely he can 
“call his shots.’’ A suggested alternative is to project the 
needle movement but not the scale. Numbers correspond- 
ing to the dc values obtained at different balance points 
are drawn in chalk on the screen, which may be a white 
board. Later, the actual ac scale values are exposed and on 
projection may be compared with the chalked-in values. 

A supplement to the above demonstration shows experi- 
mentally that Vpeak/ Vrms = V2. Figure 2 shows how a neon 
bulb can be connected alternately to sources of variable 
ac and dc. The demonstration consists in connecting the 
bulb to dc and raising the voltage from zero until the bulb 
first lights. This voltage is recorded as V;. Then the neon 
bulb is connected to ac, and the voltage is raised from zero 
until the bulb first lights. The ac voltmeter reading, which 
gives the rms value, is recorded as V2. Having demon- 
strated several times the fact that the bulb does not light 
until the potential across it is at least Vi, the student may 
be convinced that Vi=Vyeax and V2=Vrms. In several 
trials we have obtained values of about 70 volts con- 
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Fic. 2. The minimum striking voltage of a neon bulb may 
be used to demonstrate that Vpeak =VrmsV2. 


sistently for V; and about 50 volts for V2. This yields an 
experimental value of 1.4 for the theoretically predicted v2. 


; ALBERT V. BAEz 
UNESCO Technical Assistance Mission 


University College 
Baghdad, Iraq 


1 Joly diffusion photometer. Sutton, Demonstration Experiments in 
~~ (McGraw-Hill Book Company, Inc., New York, 1938), L-12, 
p. b 


A Laboratory Experiment on Relationships 
between ac and dc Voltage 


HE following instructions have been used successfully 

for an experiment based on the letter entitled: 

“A demonstration on relationships between ac and dc 
voltage.” 

A. 1. Arrange parts and make connections according to 
Fig. 1 (see preceding letter). 

2. Before proceeding, however, connect both lamps 
temporarily to the same autotransformer output. With 50 
volts ac on both lamps, balance the paraffin photometer by 
moving it between the lamps. Repeat this procedure at 70, 
90, and 110 volts ac. Once you have found the best balance 
position, leave the paraffin block fixed. 

3. Now connect one lamp to dc and the other to ac. 
Start at 110 volts dc and adjust the autotransformer on 
the ac side until both lamps are equally bright, as judged 
photometrically. Record the readings on dc and ac volt- 
meters. 

4. Repeat (3) at dc voltages of 100, 90, 80, 70, and 60 
volts. Record the ac voltmeter reading. 

5. Calculate the percent difference in each case and the 
average percent difference. 

B. 1. Make connections to determine the minimum dc 
potential which will cause the neon bulb (see Fig. 2 of 
preceding letter) to light. 

2. Determine and record several times the minimum dc 
voltage needed to light the neon bulb. 

3. Repeat (2) with ac voltage. 

4. Determine the ratio of Vpeak/ Vrms from the above as 
suggested by the preceding discussion. 

ALBERT V. BAEZ 
UNESCO Technical Assistance Mission 


University College 
Baghdad, Iraq 
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ANNOUNCEMENTS AND NEWS 


Book Reviews 


Radiations from Radioactive Substances. Sir ERNEST 
RUTHERFORD, JAMES CHADWICK, AND C. D. ELLIs. 
Reissue of the edition of 1930. Pp. 588+xi, 1624 cm. 
The Cambridge University Press, New York, 1951. 
Price $11.00. 


This classical treatise on the early work of radioactivity 
has been out of print for many years. The current reissue 
has been corrected for the previous few typographical 
errors. 

The presentation of the theory of successive transforma- 
tions in Chapter I has not yet been improved upon by 
current authors of nuclear books. Although some of the 
numerical values associated with the branch products of 
various chain reactions have been changed in recent years, 
the fundamental points of view are most clearly presented. 
All alpha-particle phenomena are thoroughly treated, such 
as scintillations, ~ange, specific ionization, stopping power, 
straggling, scattering, radioactive haloes, and even physio- 
logical actions. 

The presentation of the general collision relations of 
nuclei, leading to the concept and evaluation of a potential 
field about a nucleus, can hardly be improved upon. 
Chapters 8 and 9 are, therefore, very sound and important 
reading for nuclear students today. 

The chapters on artificial disintegration by alpha- 
particles and on beta-decay are primarily of historical 
interest, since no high energy accelerators existed in 1930. 
By using beta-ray spectrographs, Geiger counters, gamma- 
ray spectrometers, etc., although somewhat crude accord- 
ing to modern standards, nuclear energy levels were, never- 
theless, arrived at for numerous radioactive nuclei. 

The latter portion of the 588 page book is devoted to 
more recent (1930) discoveries such as radioactivity of 
ordinary matter, new scintillation materials, recoil method 
in chemical separations, and the structure of the nucleus. 

Radiations from Radioactive Substances cannot, of course, 
be used as a text in a beginning course in nuclear physics, 
but the contents of this book should formulate the basis 
of much of the material presented in such a course. How- 
ever, for nuclear research students, nuclear scientists, and 
teachers of nucleonics this book should be among those im- 
portant ones that are within easy reach. 


M. L. PooL 
The Ohio State University 


Eclipses of the Sun. SAMUEL ALFRED MITCHELL. Fifth 
Edition. Pp. 445+-xv, Figs. 48, 15.523 cm. Columbia 
University Press, New York, 1951. Price $6.50. 


This interesting and instructive book is now in its fifth 
edition and includes the major solar researches of recent 
years. Dr. Mitchell is a veteran of eclipse expeditions, 
having taken part in ten of these. He has made photo- 
graphic and spectrographic records during twenty-six 


minutes of totally eclipsed sun, and has published extensive 
data dealing with the wavelengths and intensities of 
spectral features of the solar chromosphere. In recent years 
this body of data has been re-examined by Wildt and has 
yielded important results. 

A number of chapters are devoted to a discussion of 
ancient eclipses, and several deal with the author’s work 
and experiences at recent eclipses. The material on early 
eclipses is most interesting and contains some authentic 
Biblical history. The reader will acquire a fair knowledge 
of the problems of an eclipse expedition from the discussion 
of the more recent eclipses. Dr. Mitchell has managed to 
convey something of the excitement and strain which 
exist during the period of totality. 

Several chapters are devoted to a discussion of spectro- 
scopic equipment and elementary spectroscopic theory. 
The chapter on spectroscopes could well have been ex- 
panded to include material on modern gratings and their 
possibilities. One small error may be noted. On page 77 
the word extinction rather than absorption should be used. 
There is an important physical distinction between these 
processes. In chapters 10 and 11 more emphasis is required 
on the processes of excitation and ionization since the 
“theart”’ of the book deals with the results of these processes 
operating in different regions of the solar atmosphere. 

Those interested in the physics of the solar atmosphere 
will consider chapter 13-18 the important part of the book. 
Dr. Mitchell has reviewed carefully the large amount of 
literature dealing with problems of the reversing layer, the 
chromosphere, and the corona, and the careful reader will 
be brought up to date on this subject. The reviewer felt 
that the discussion would be strengthened by the inclusion 
of formulas which illustrate the calculation of various 
physical conditions. For example, the equations which lead 
to the calculation of excitation temperatures from the 
intensities of emission lines would add weight to several 
sections of the material. 

Dr. Mitchell’s book will be found useful by both the 
research worker and those who have only a general interest 
in the field. 

W. PETRIE 
University of Saskatchewan 


Principles of Radio. Keira HENNEY AND GLENN A. 
RICHARDSON. Sixth Edition. Pp. 655+-vii, Figs. 405, 
14.522 cm. John Wiley and Sons, Inc., New York, 
1952. Price $5.50. 


This is the sixth edition of a well-known book on the 
principles of radio. Although it is quite technical, it assumes 
little previous knowledge on the part of its reader. The first 
chapter starts with such elementary topics as the plotting 
of one variable against another, and the use of powers of ten 
to express large and small numbers. About one-third of the 
whole book is used to cover the fundamental theory of dc 
circuits, electrical measurements, electromagnetism, and 
ac circuits. The remainder of the book covers all the 
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common applications of electronic devices in radio, in- 
cluding television and radar. 

As stated in the author’s preface, the book is written 
from a viewpoint useful to those who must learn radio 
without the help of a teacher. Considering the great amount 
of material included in the book’s six hundred and fifty- 
five pages, the authors are to be commended on a truly 
remarkable achievement of detailed and clear explanation. 
The readers for whom it is intended will find it an excellent 
text. It would make an excellent textbook for use in a trade 
school or technical institute. Readers with a background of 
college algebra, trigonometry, and physics would find it less 
satisfactory as a textbook. For them the simpler concepts 
are over-explained, and the more difficult topics are treated 
too briefly. However, such people will find it valuable as a 
quick reference book which may give them all the informa- 
tion required, for a given purpose. For such reference 
purposes, the book is up-to-date and remarkably complete. 

Solved examples of problems are included throughout the 
text along with a large number of problems for which the 
answers are not given. In one of the earlier chapters, brief 
directions for a number of experiments on resonant circuits 
are given, following the style of a laboratory manual. There 
did not appear to be any other place in the book where such 
instructions were included, and no reference to these in- 
structions as such was made in the table of contents, in the 
preface, or in the index. 

In addition to the chapters dealing with the application 
of electronic tubes in amplifiers, rectifiers, and oscillators, 
the book contains chapters treating transmission lines, 
antennas, electromagnetic radiation, frequency modulation 
and detection, ultra-high-frequency phenomena, electronic 
instruments, transients and wave-shaping circuits, tele- 
vision and radar. Listed in the complete index are found 
items of recent interest such as klystrons, lighthouse 
triodes, television, ghosts, pulse modulation, etc. Inside the 
front cover is a table of the characteristics of the more 
common types of receiving tubes. Inside the back cover 
is a table of wavelength, frequency, and L—C combina- 
tions, and a table of trignonometric functions. 


JoserpH H. Howey 
Georgia Institute of Technology 


Principles of Quantum Mechanics. WILLIAM V. Houston. 
Pp. 283+vii, Figs. 21, 23.5X16cm. McGraw-Hill 
Book Company, Inc., New York, 1951. Price $6.00. 


“Quantum mechanics appears as a difficult and some- 
what mysterious subject to a beginner. The mystery 
gradually decreases as more examples are tried out, but one 
never quite loses the feeling that there is something peculiar 
about the subject.’’ This frank remark made recently by 
an eminent quantum mechanician emphasizes the special 
challenge which quantum mechanics presents to teachers, 
students, and textbook writers. 

The difficulties which arise in a first encounter with 
quantum mechanics are of twofold origin: namely those 
which, in effect, the student brings with him, and those 
which may be properly attributed to specifically new 
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quantum concepts. This reviewer feels that the elimination 
of the difficulties of the first type should be attempted 
much earlier in a physics student’s studies than at the 
beginning of his quantum mechanics course. Indeed most 
of these difficulties have to do with questions which an 
inquiring student might easily have raised in any branch 
of classical physics were he encouraged to do so, e.g., what 
should one expect of a physical theory? What assumptions 
are implicit in the processes of measurement? What is the 
status of such concepts as energy, field, inertia, etc., which 
he is taught to manipulate so freely? Such philosophically 
tainted questions are usually ignored as being of no use in 
most undergraduate physics instruction only to appear, 
seemingly of necessity, in the introductory quantum theory 
course where the student meets them for the first time 
amidst a rapid succession of new physical phenomena, 
mathematical techniques, and concepts which hardly 
afford him the opportunity to find his bearings. 

In this textbook Professor Houston meets this problem 
on a rather modest scale. His words on the philosophical 
matters are generally few in number, clear, and carefully 
stated. If they do not thoroughly explore all aspects of the 
matter at hand, neither do they leave the reader with any 
glib or dogmatic conclusions, e.g., the reader is told that 
other interpretations of quantum mechanics might be 
possible. 

The reference to ‘‘Principles” in the title of the book is 
not entirely appropriate since most of the book is devoted 
to special methods and applications of quaantum-mechanical 
calculations, rather than to any systematic elucidation of 
general principles. The author states that his object has 
been “‘to acquaint the student with a formulation of 
quantum mechanics such that he can make applications 
for himself and can understand published work.” The for- 
mulation he has chosen is the wave mechanics of the 
Schrédinger equation. Sometimes matrix elements come 
sufficiently to the fore, as in the chapters on spectroscopy, 
to make one wish that the author had bothered to show 
the relationship between the wave and matrix forms of 
quantum mechanics. It is questionable whether such an 
omission will help the student in ‘“‘understanding published 
work.”’ Similarly one might doubt the wisdom of avoiding 
a systematic presentation of some of the more abstract 
methods when these same methods are pulled out of the 
hat for special problems and thereby appear as isolated 
artifices rather than as a part of a general method. An 
example of this is the use of symmetry operators in dis- 
cussing lattice wave functions, or the parity operator, 
which is first introduced in connection with two-electron 
spectra. How much better it would be if the concept of 
invariance under a group of symmetry operators were in- 
troduced at the beginning and applied systematically to 
deduce such consequences as the conservation laws, selec- 
tion rules, degeneracy of wave functions, and constants 
of motion which could then be illustrated with each of 
the special systems considered later! The absence of a 
systematic exposition of linear algebra or transformation 
theory invites similar objections. 













































The book provides a generous selection of topics to 
acquaint the student with the power and breadth of 
quantum mechanics. Part I, General Outline of Quantum 
Mechanics, contains most of the standard topics of non- 
relativistic Schrédinger theory, including the Pauli theory 
of electron spin. Part II, Applications to Spectroscopy, has 
considerable detail not to be found in other texts at this 
level on two-electron spectra. Part III, Applications to 
Collision Problems, gives a summary of scattering theory. 
Part IV, Applications to Electrons in Solids, is of unusual 
interest. The chapters on the motion of the electrons in a 
periodic potential includes a treatment of Brillouin zones, 
and will be welcomed by many. On the other hand, a 
chapter on quantum-statistical mechanics leaves much 
room for improvement, primarily because it is all based on 
the use of a quantum-mechanical density operator whose 
presence is justified only by a footnote referring to Tolman’s 
monumental book which no student is likely to have 
mastered. In Part V, Application to Electromagnetic 
Radiation, Professor Houston has taken care to indicate 
the classical analogies which help make the quantum treat- 
ment more amenable to the student’s classical intuition, 
and the student will be able to learn from this discussion 
(which follows Fermi) much more about some properties of 
the electromagnetic field than can be found in the standard 

classical texts on electromagnetic theory. 

One of the rewarding features of the author’s presenta- 
tion is the interlacing of problems with text. These prob- 
lems, as in his book on mathematical physics, are never 
trivial, and not only give the student a chance to test his 
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Kentucky Section 


The spring meeting of the Kentucky Section, American 
Association of Physics Teachers, was held on May 17, 1952, 
in the Physics Building, University of Kentucky, Lexington, 
Kentucky. The meeting was attended by 35 members and 
guests. Dr. PAuL C. OVERSTREET, Morehead State College, 
presided. 

After the presentation of the contributed papers, the 
laboratories and shop facilities were open for inspection. 
The meeting concluded with a luncheon in the Student 
Union Cafeteria. The contributed papers were as follows: 


1. Energy spectra of slow neutron-capture gamma-rays. 
R. J. REITHEL AND B. D. KERN, University of Kentucky.— 
A photographic technique suggested by Pringle and Isford! 
is being used to measure the energy spectra of slow neutron- 
capture gamma-rays in manganese, lead, and magnesium. 
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understanding of the text, but also serve to complete and 
extend the discussion. 


D. L. FALKOFF 
Massachusetts Institute of Technology 


New Members of the Association 


The following persons have been made members or junior members 
(J) of the American Association of Physics Teachers since the pub- 
lication of the preceding list [Am. J. Phys. 20, 386 (1952)]. 


Adelson, Harold E. (J), 2077 Vyse Ave., Bronx, New York 
60, N. Y. 

Banner, Richard Duer (J), 470 E. Mill St., Liberty, Mo. 

Bird, John B., 35 Ronan St., Hillcrest, Binghamton, N. Y. 

Colodzin, Martin (J), c/o Brecher, 700 Ocean Ave., 
Brooklyn, N. Y. 

Conway, Robert Lawrence (J), 3751 College Ave., San 
Diego, Calif. 

Forrester, Robert Louis (J), Rt. 4, Box 2084, Vista, Calif. 

Francis, Sister Marguerite, College of St. Elizabeth, Con- 
vent Station, N. J. 

Hall, Lewis Summerlin (J), 160 Barber St., Athens, Ga. 

McCoy, Dallas Perry (J), 1203 Berly, Pacific Beach, San 
Diego 9, Calif. 

Saenger, Alfred (J), 2162 University Ave., Bronx, New 
York, N. Y. 

Stow, Richard Wolcott, Mayowood, Route 2, Rochester, 
Minn. 

Ulloa, Fernando Duran (J), 2859 Broadway, San Diego 2, 
Calif. 
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A 50-mc Ra-Be source in a slow neutron howitzer? supplies 
thermal neutrons. The scintillation spectrometer consists 
of a NaI-Tl crystal, an RCA 5819 photomultiplier tube, a 
linear amplifier, and a Tetronix Model 511A oscilloscope. 
The 1.17-Mev and 1.33-Mev gamma-rays from Co® are 
used for calibration purposes. 

The face of the oscilloscope is photographed by a 
Dumont oscillograph record camera using Ansco Supreme 
film. The light transmission of the negative is measured on 
a Hilger comparator-microphotometer which has been 
converted into a recording comparator microphotometer.® 
This technique is checked by comparison of the microphoto- 
meter trace with a graph obtained by use of a single-channel 
differential discriminator. The width of the Cs!*7 gamma- 
ray line is the same by both methods. 

1R. W. Pringle and G. Isford, Phys. Rev. 83, 467 (1951). 


2 Dana P. Mitchell, Phys. Rev. 49, 453 (1936). 
3 Reithel, Kern, and Hanau, Am. J. Phys. 20, 388 (1952). 
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2. A study of the spectrum of the cathode spot of the 
mercury arc. ROBERT E. WEAVER AND S. V. GALGINAITIS, 
University of Louisville—A discharge tube was designed 
and built in such a way that the cathode spot of a mercury 
arc could be viewed from above and focused on the slit of 
a spectrograph for as long a time as desired. The spectro- 
graph used was a Hilger medium quartz spectrograph. The 
spectrograms verified the presence of a continuous spec- 
trum. By comparison with the spectra from other parts of 
the arc, it was definitely determined that the continuum 
originates at the cathode spot. 

The most intense continuum occurred in two principal 
regions. The maximum intensities in these regions occurred 
at approximately 4850 and 3420A. The positions of these 
maxima did not seem to depend on the current passing 
through the arc. It will be noted that the wavelengths 
mentioned above are approximately those of the heads of 
two of the molecular bands of mercury. 

Further study is planned, using a grating spectrograph. 
It is hoped that with higher dispersion, the exact nature 
of the continuum can be better determined. The ultimate 
hope is to gain more knowledge concerning the nature of 
electron emission at the mercury cathode. 


3. Spectrographic analysis of archaeological artifacts. 
RicHARD Hanau, University of Kentucky.—A preliminary 
qualitative spectrographic analysis was made of six copper 
bracelets found in eastern Kentucky Indian mounds. The 
samples, coated with green and black encrustations, were 
cleaned in HCl. Assuming the coating to be 80 percent 
malachite and 20 percent CuO, the average amount of 
original Cu oxidized is about 18 percent. There appears 
no significant difference in this value between samples 
from different mounds. The cleaned samples show crevices, 
etched in black CuO, indicating that the thicker bracelets 
were probably made by rolling thin sheets, or hammering 
thin strips together, to form a rod. 

The samples, ranging in size from 3%; in.-diameter elec- 
trodes to small slivers imbedded into graphite electrodes, 
were burned in a 2300-volt, 1-ampere, ac arc, gap width 
? mm. Ten-minute exposures indicate seven noncopper 
lines between 2150 and 4450A. Two are the very sensitive 
Ag lines, 3280 and 3382A. Five are air lines or not positively 
identified. The spectra show no qualitative difference be- 
tween samples. The method of forming the bracelets, the 
constitutional similarity between samples, and the high 
purity of copper, agree with results of metallographic 
studies on these same bracelets.! 


1 James E. Fox (unpublished material). 


4. Recent modifications of the University of Kentucky 
1-Mev Van de Graaff generator. H. H. Givin, J. L. Ryan, 
L. W. Cocuran, T. M. Haun, ann B. D. KERN, University 
of Kentucky.—Recent modifications of the University of 
Kentucky Van de Graaff generator include: (1) employ- 
ment of a mixture of carbon dioxide and nitrogen in the 
pressure system; (2) replacement of the corona ladder with 
a voltage divider consisting of S.S. White 200-megohm 
resistors; (3) collimation of the beam at the ion source; (4) 
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installation of a corona feedback circuit, capable of 
voltage stabilization within 0.25 percent; (5) installation 
of a 90° magnetic analyzer. 

These modifications have resulted in greater stability, 
improved focusing, and increased ease of operation. 


5. Calibration of the University of Kentucky 1-Mev 
Van de Graaff generator. J. L. Ryan, H. H. Givin, L. W. 
Cocuran, B. D. KERN, AND T. M. Haun, University of 
Kentucky.—Calibration of the Kentucky Van de Graaff 
generator was effected using the well-known Li’(p, ~)Be® 
resonance! at 440 kev and the resonances below 1 Mev 
seen in the F!9(p, ~)Ne®® reaction.? The calibration was 
made in terms of the magnetic analyzer current. 


sue Lauritsen, Morrison, and Fowler, Revs. Modern Phys. 22, 
1950). 
2 Chao, Tollestrup, Fowler, and Lauritsen, Phys. Rev. 79, 108 (1950). 


6. Measurement of the gamma-ray yield from proton 
bombardment of boron. L. W. Cocuran, J. L. Ryan, 
H. H. Givin, T. M. Hawn, Ann B. D. KERN, University of 
Kentucky.—Thick and thin target gamma-ray yields from 
the proton bombardment of natural boron have been 
measured in the energy region between 130 kev and 1 Mev 
using the University of Kentucky electrostatic generator. 
From the thick target data an upper limit of 7107 
gamma-ray per proton was obtained for the reported! 
resonances at 850 kev and 950 kev. Analysis of the thin 
target data indicates the presence of a broad resonance 
at a resonant energy of 670+15 kev. with a width of 390 
kev. 


1 Kerb, Kerst, and McKibben, Phys. Rev. 51, 691 (1937); Curren, 
Dee, and Petrzilka, Proc. Roy. Soc. (London) 169, 269 (1939). 


7. Discussion of the gamma-ray yield from proton 
bombardment of boron. L. W. Cocuran, H. H. Givin, 
J. L. Ryan, B. D. KERN, AND T. M. Hawn, University of 
Kentucky.—The broad 390-kev resonance observed at a 
resonant energy of 670+15 kev would seem to indicate 
that the observed gamma-rays are from an excited state 
of Be® following alpha-particle emission from an excited 
state of C2. However Walker! has observed the energy 
spectrum of gamma-rays from the proton bombardment 
of boron and did not observe gamma-rays of energies other 
than those associated with the excited states of C. 


1R. L. Walker, Phys. Rev. 79, 172 (1950). 


Officers of the Section are: President, Paut C. OvER- 
STREET, Morehead State Teachers College; Vice-president, 
Cart Apams, University of Louisville; Secretary-Treasurer, 
RICHARD HANAu, University of Kentucky; Representative 
to the Council of AAPT, R. A. Lorinc, University of 
Louisville. 

RicHARD HANAU, Secretary-Treasurer 


Michigan Teachers of College Physics 


Kalamazoo College and Western Michigan College were 
joint hosts for the spring meeting of the Michigan Teachers 
of College Physics on May 17, 1952. Registration included 
83 members and wives. 
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The morning session in McCracken Hall, the new 
science building of Western Michigan College, opened with 
an account by Dr. E. R. PHEeLps of Wayne University of 
his activities as Detroit’s ‘‘TV Weatherman.” Dr. E. R. 
GaRRETT of the Upjohn Company presented a lecture on 
some of the many uses of physics in the pharmaceutical 
industry. After a discussion of physics in general education, 
led by Dr. R. SCHLEGEL and Dr. V. W. MILEs, the group 
was joined for lunch by the wives, who had been touring 
the campus and following their own program on music 
and art. 

Dr. WILLARD LisBy of the University of Chicago was the 
principal speaker of the afternoon session at Kalamazoo 
College. After reviewing techniques for measuring car- 
bon-14 and the theory of its production in the atmosphere 
by cosmic-ray neutrons, he presented some of the findings 
in dating and other applications. Dr. H. R. CRANE then 
gave a review of some of the highlights of the Washington 
A.P.S. meeting, notably work with “V-particles” and the 
status of the Brookhaven “‘cosmotron.” A concluding tea 
was held at the Barbours’ home on the campus. 

Tan G. BARBOUR 
Kalamazoo College 


Oregon Section 


The sixtieth meeting of the Oregon Section of AAPT was 
held at the University of Oregon on April 26, 1952. The 
President of the Section, Dr. WILLIAM L. PARKER, presided. 
Invitations were accepted from the following institutions 
for meetings during 1952-53: Fall, University of Portland; 
Winter, Linfield College; Spring, Oregon State College. 

Officers for the coming year were elected as follows: 
President, ROBERT L. PurRBRICK, Willamette University; 
Secretary, KENNETH E. Davis, Reed College; Historian, 
B. GODFREY VASSALLO, University of Portland; AAPT 
Council Member, Witt V. Norris, University of Oregon. 


Program 


The following papers were presented at the general 
sessions held in connection with the dedication of the new 
Science Building at the University of Oregon: 


Problems in modern biology. G. W. BEADLE, California 
Institute of Technology. 

Chemical institutions and scientific personnel in Europe. 
G. Ross RoBErtson, University of California, Los Angeles, 

Current research in nuclear physics. S. K. ALLISON. 
University of Chicago. 
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Education for science, panel discussion. S. K. ALLISON, 
G. W. BEADLE, ALAN T. WATERMAN, DONALD STOTHER, 


CLARENCE DIEBEL, P. J. VAN RyYSSELBERGHE, E. G. 
EBBIGHAUSEN. 


The following contributed papers were presented: 


1. An electrical method for solving the lens equation. 
R. L. Pursricx, Willamette University—An electrical 
method for solving a reciprocal equation of the lens- 
formula type was described. The operation of the device 
is based upon Ohm’s law and employs three sliding resis- 
tors and nomographic principles. 


2. An aberration produced by crossed cylindrical lenses. 
J. J. Faris, Pacific University—A spherical wave front 
refracted by a toric lens is transformed to a toric front. 
After this wave front has proceeded a distance ¢ beyond 
the lens, the wave front is no longer toric. If cF, and cF, 
are both small compared to one, the difference between the 
curvature of this wave front and a toric surface can be 
expressed approximately as 


A(6) = —c(F.— F,)? sin*26 


where F, and F, are the powers of the principal meridians 
of the lens and @ is measured from the y meridian. 

Since this wave front is not toric, a toric lens cannot 
convert it to a spherical wave front. 


3. V-particles. LAWRENCE S. GERMAIN, Reed College. 


4. Electric breakdown of solid dielectrics. HENRY R. 
KAISER, Oregon State College. 


5. Field emission and the vacuum arc. W. W. Doan, 
J. K. TROLAN, and W. P. Dyke, Linfield College. 


6. New absorption bands of alkali atoms due to the 
presence of foreign gases. SHANG-Y1 Cn’EN, University 
of Oregon.—A resumé of the data, about the narrow and 
diffuse absorption bands of alkali atoms due to the presence 
of foreign gases, appearing on the short wavelength side of 
certain members of their principal series, was given. The 
data included in addition to the case for Rb and Cs as 
already published also those for Li, Na, and: K in the 
presence of He, Ne, A, Nz, and He as recently observed. 
Theoretical interpretations of the existence of these bands 
were discussed. 


Frep W. DECKER, Secretary 


To spin words and make pretty sentences is my one talent, and I must make the best of it. And I 
am fallen on an age when not one man in 20,000 knows good English from bad, and not one in 
50,000 thinks the difference is of any importance-—D’ARcy WENTWORTH THOMPSON, quoted 
by JoHN TYLER BONNER, Scientific American, August 1952, p. 63. 
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HE 1952 summer meeting of AAPT met 
with the Colloquium of College Physicists 
at the University of Iowa, June 11, 12, 13, and 14. 
This was the 14th annual meeting of the Collo- 
quium,—an event sponsored by the University 
of Iowa under the inspired direction of Professor 
G. W. Stewart. Approximately 200 physics 
teachers from 100 colleges and universities, from 
Massachusetts Institute of Technology on the 
east coast to Stanford on the west, attended the 
meetings. 
Some 30 short contributed papers, whose ab- 
stracts appear below, and 3 invited papers com- 
prised the AAPT portion of the program. 


Invited Papers 


The Raydist system. VERNE F. Swarm, Naval Aviation 
Ordnance Test Station. 

A survey of general physics laboratories in the United 
States. SANBORN C. Brown, Massachusetts Institute of 
Technology. 

Movies of magnetic domains. R. M. Bozortu, Bell 
Telephone Laboratories. 

The program of the Colloquium portion of the meeting 
included the following lectures: 

Operations research. W. L. Wuitson, Johns Hopkins 
University. 

Edge tones and whistles. H. K. ScHILiinG, Pennsylvania 
State College. 

Recent developments in geophysical prospecting. H. B. 
PEACOCK, Geophysical Service, Incorporated. 

Researches in magnetism. R. M. Bozortu, Bell Tele- 
phone Laboratories. 

Demonstration experiment round table—participated in 
by R. H. Cook, South Dakota School of Mines; H. C. 
JENSEN, Lake Forest College; M. J. Pryor, New York State 
College for Teachers; Paul Roop, Western Michigan College 
of Education. Each of the speakers described a method of 
presenting a demonstration experiment, emphasizing the 
reasons for the methods used. 

Probably the most unique feature of the Colloquium 
was the exhibit of new devices. More than 35 items were 
brought by as many members of the Colloquium. An 
anonymous donor provided cash prizes for the best exhibits 
in each class of devices,—experimental and nonexperi- 
mental. These were determined by ballot and were awarded 
with suitable public recognition at the Friday luncheon. 
The prize winners were as follows: 


Experimental 


First prize, PROFESSOR FRANK VERBRUGGE, Carleton 
College. 
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Second prize, PROFESSOR RICHARD Howe, Denison 
University. 

Third prize, PROFESSOR HAYM KRUGLAK, University of 
Minnesota. 

Honorable mention: 
PROFESSOR PAUL Roop, Western Michigan College. 
ProFeEssor J. G. Winans, University of Wisconsin. 
P. W. Wituiams, South Dakota State College. 
Dr. RosE A. Carney, St. Procopius College. 
Dr. R. R. Hancox, Midwest Research Institute. 
PRoFEssOR Harry Mason, Jamestown College. 


Nonexperimental 


Mr. Paut J. Haicu, University of Illinois. 
PROFESSOR I. WALERSTEIN, Purdue University. 


The Colloquium also included a series of lectures spon- 
sored by the Research Corporation and presented by 
GeEorGE E. UHLENBECK of the University of Michigan. 
These were entitled The so-called elementary particles 
and Some famous unsolved problems in statistical physics. 

The usual informal social evening at the home of 
Professor Stewart,—an event which requires the use of his 
lawn—had to be called off because of rain. It was trans- 
ferred with considerable success to one of the lounges in 
the Iowa Union. On this occasion, when called upon for 
some informal remarks, Dr. Uhlenbeck responded with a 
brief anecdotal history of the discovery of the spin of the 
electron, recounting how he and S. A. Goudsmit came to 
propose the idea of spin while still students of Professor 
Ehrenfest at Leyden. 


Contributed Papers 


1. Catalysis of students for careers in physics. LEONARD 
O. OLSEN, Case Institute of Technology.—This paper con- 
sidered the type of curriculum best suited for the education 
of students at the bachelor level, who are primarily inter- 
ested in entering employment after securing this initial 
degree. The opportunities for such employment were also 
considered. 


2. Triangular inches as a bias detergent for students of 
physics. T. D. PurLiips, Marietta College——The measure- 
ment of area by the use of ‘triangular inches’’ as units was 
suggested as a recreational project. A ‘‘triangular inch’’ is 
defined as an equilateral triangle one inch on a side. On 
this basis formulas for computing areas of a number of 
simple figures were developed to illustrate the possibilities. 
As a companion project an unconventional trigonometry 
was described. In this example the elementary functions 
were defined relative to a sixty-degree triangle rather than 
to the right triangle. The ulterior motive of the project was 
to let the student demonstrate for himself that ‘‘we can do 
things differently if we want to.” 















3. Quaternions in undergraduate physics. J. G. WINANS, 

University of Wisconsin—A quaternion or scalar plus 
vector (Q=S+V) results from the multiplication or divi- 
sion of arbitrary vectors. One vector a is replaced by the 
sum of two components, one parallel and one perpendicular 
to the other vector b. Thus a=a;+a2, and ab=a,b+a»b, 
or ab=S-+V, since a; is parallel to b and a: is perpendicular 
to b. Now a,:b=ab cos@= (a-b), and a2b=ab sind=[aXb]; 
so ab=(a-b)+[aXb]. Also a/b=a,/b+a:/b=S2+V2. 
Distributiveness has been assumed. Cos#=adjacent side/ 
hypotenuse is a quaternion as is also sin but tané =oppo- 
site side/adjacent side is a vector. Since sin*@+-cos*@ = 1, we 
have Q,?+ Q.?=1, and the square of a quaternion is a scalar. 
(S+V)?= Q?, or S2+SV+VS+V?= Q?. Separating scalars 
and vectors leaves S?+V?=Q?, and SV+VS=0. Thus 
multiplication of scalar and vector is not commutative. 


4. A simple method for obtaining numerical solutions 
to problems governed by the one-dimensional wave equa- 
tion. WALTER P. Rep, Naval Ordnance Test Station, China 
Lake, California.—If y(x, t) is the general solution to the 
one-dimensional wave equation, then 


y(x+vh, t-+H) =y(x+0H, t+h) 

+y(x—vH, t—h) —y(x—vh,t—H). 
This formula may be new, although it is so simple that one 
would expect it to be well-known. At any rate, it can be 
shown that its use will frequently enable one to calculate 


very easily numerical results of problems governed by the 
one-dimensional wave equation. 


5. Film-loops for physics teaching. JoHn J. HEILMAN, 
Ursinus College, Collegeville, Pennsylvania.—Continuous 
film-loops with running times of 20 to 120 seconds were 
shown illustrating phenomena in vibratory motion, wave 
motion, and impact. These loops may be used either as 
animated lantern slides in lecturing or by the students to 
make quantitative studies of the phenomena. 


6. A simple high frequency demonstration oscillator. 
RIcHARD H. Howe, Denison University.—A high frequency 
demonstration oscillator of simple construction was ex- 
hibited and briefly described. The oscillator was used to 
demonstrate several of the fundamental uses of electro- 
magnetic radiation with particular application to high 
frequency arrays. 

The power output of the oscillator is sufficient to provide 
ample visual display of the effects to a large group. It 
operates at a frequency which permits the demonstration 
equipment to be of dimensions convenient for class- 
room use. 


7. Echo ranging with audiofrequencies. GEORGE BRAD- 
LEY, Western Michigan College of Education.—A device has 
been constructed which detects objects in a room by means 
of a reflected sound signal. The sound pulse is made up of 
very few cycles of audiofrequencies, and therefore the 
amplifiers and detectors operate in the audio range. The 
device, in using only a few cycles in each pulse, gives 
graphic demonstration to the student that not all sound 
waves are periodic. 
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This apparatus was primarily designed for a demonstra- 
tion piece for the classroom. It has been utilized success- 
fully (a) to measure the speed of sound in free air, (b) to 
illustrate the reflection of sound waves, (c) to measure the 
reverberation time of enclosures, and (d) to illustrate the 
principle of sonar and radar detection in the classroom. 
Aside from a standard power supply, the apparatus 
utilizes 8 tubes in the electronic circuits, mounted on a 
7 in.X7 in. chassis. A small permanent magnet speaker 
acts as the transducer and a dynamic microphone as the 
detector. Data were presented on a standard oscilloscope. 


8. Preparation of pictorial material for classroom use. 
Z. V. Harvaik, IST, University of Arkansas.—Photo- 
graphic methods are employed to assemble pictorial 
material. Emphasis is given to the use of 35-mm film 
because of the easy handling, low costs, and the small 
storage space required for the slides. 

Material was demonstrated which dealt with tables, 
graphs, apparatus layouts, experimental records, and 
microfilming. 

A simple apparatus was described and shown which 
enabled almost routine-type speed in taking photographs of 
material useful in lectures, demonstrations, and reference. 


9. Color demonstration with a small projection lantern. 
V. E. Eaton, Wesleyan University—A Model MC De- 
lineascope and a transmission grating were used to project 
the continuous spectrum of a lamp source and the band 
spectra of various filters. By substituting a rutile prism for 
the grating, two large brilliant spectra were produced. 
The same lantern was used to demonstrate the subtractive 
method of color photography. 

For other color experiments the light from the lantern 
was reflected onto the screen from three mirrors, mounted 
on universal joints and separately controlled. When three 
primary filters were properly mounted in the lantern the 
arrangement could be used as a color mixer and the in- 
tensity of each component could be controlled by means of 
Polaroids. This arrangement could be used also to demon- 
strate the additive method of color photography and to 
project the Maxwell color triangle. 


10. Individual projects on the intermediate level as a 
training aid and as an economy measure. O. B. Youn, 
Southern Illinois University—An outline of advantages of 
individual projects performed by undergraduate physics 
students on the intermediate level was presented. The 
method is highly adaptable and flexible, meeting different 
types of student needs and departmental circumstances. 
Twenty-three different points were considered. Special 
emphasis was given to student motivation, development 
of interest in physics, economy of teacher time, of depart- 
mental equipment, and application to students of varying 
abilities. 


11. A most unusual solar phenomenon witnessed at 
Cedar Falls, Iowa, February 27, 1952. R. A. ROGERs, 
Iowa State Teachers College—From about 11:00 A.M. 
until 2:00 P.M. on_February 27, 1952, conditions ap- 
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parently were well-nigh perfect in the Cedar Falls—Water- 
loo area of Iowa for an extremely fine and unusual display 
of solar halos. Portions of the phenomenon were visible 
also, in less intensity, throughout a large portion of Iowa 
and as far north and east as Madison and Beloit, Wis- 
consin. The display included both the 22° and the 46° halos, 
the parhelic circle, brilliant parhelia or sundogs, tangent 
arcs of the 22° halo, lateral tangent arcs of the 46° halo, 
parhelia of 120°, and other less certainly identified phe- 
nomena. The phenomenon was unusual not only because 
of its complexity but also because of its long duration and 
spectacular intensity. 


12. Light does not weaken because of its spreading. 


F. M. StEapMAN, Los Angeles, California. (To be read by 
title only). 


13. Series of lectures for high school students. J. W. 
Bucuta, University of Minnesota.—On many occasions 
college and university teachers blame the high schools for 
failure to give adequate background in the sciences, 
especially the physical sciences. Perhaps some of this 
censure is deserved, but often the college and university 
staff member does not appreciate the task of a high school 
teacher; and too often the high school teacher is given 
little aid and support by the scientists in higher education. 

In an attempt to do something definite and concrete in 
aiding the high school group, we at the University of 
Minnesota have for a number of years offered a series of 
lectures designed especially for high school students. The 
famous lectures for ‘a juvenile auditory’ of the Royal 
Institution were our guide and inspiration. The lectures 
have been well received by the students and teachers. 
Attendance averages 150. The lectures combine some enter- 
tainment along with solid, rigorous physics. 

In this effort we believe we have contributed toward 
relieving the manpower shortage in science. The efforts are 
directed at the grass-roots level, at the time when the 
careers of many of our youth are determined. 


14. A professional course in physics for teachers of 
elementary grades. W. B. Miner, Northern Illinois State 
Teachers College-—This course was organized because of our 
realization that (1) the breadth of training required of an 
elementary teacher makes it practically impossible for him 
to take any of the courses required of a physics major, (2) 
experienced teachers are crying for courses which give them 
ideas and materials they can use in their own classrooms 
—they want something besides mere educational theory. 

In taking up a new topic each student first makes his 
own equipment from materials which are usually readily 
available to him in an actual teaching situation. The lec- 
tures and discussions follow as the student feels the need 
for help. After completing the project the student takes the 
equipment to use in the school in which he is teaching. 

The comments from experienced teachers have been very 
favorable, as have also been the comments from youngsters 
taught by these teachers. Our colleagues in the department 
of education have been well pleased with the results. They 
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have not only permitted, but often encouraged, students to 
substitute this course for one in education. 


15. The atom: A unifying concept for a general educa- 
tion physics course. ALBERT J. Hatcu, New Mexico College 
of Agriculture and Mechanic Arts.—A basic requirement for 
a general education physics course is a unifying concept 
which will unite the traditional topics of physics into a 
cohesive logical structure. Such a concept can be provided 
by an elementary study of the atom. 

There are three major aspects of the atom which can be 
utilized in a course of this type. First, an analysis of the 
orbital motion of the electron in the planetary model of the 
atom can be used as a basis for the study of the classical 
topics of mechanics and electrostatics. Second, the Bohr 
theory of atomic structure and behavior serves to unite 
several additional topics of classical physics with quantum 
physics. Finally, the natural extension of the study of the 
atom to nuclear phenomena leads to such important topics 
as mass-energy equivalence, nuclear fission, and nuclear 
forces. 

A one-semester general education physics course using 
the atom as the central theme has been presented to Arts 
and Science students at New Mexico A. and M. A. College 
for the past eight semesters with very satisfactory results. 


16. The physics library. WALLACE A. HILTON and OPAL 
R. Carin, William Jewell College-—A study of library 
facilities for students and staff members in physics depart- 
ments of colleges and universities offering at least a 
bachelor’s degree in physics was presented. The study was 
based on returns from a questionnaire which was sent to the 
physics department chairmen of 249 selected colleges and 
universities. 


17. A simplified cloud chamber for the physics labora- 
tory. ELMER Nusssaum., Taylor University—This paper 
dealt briefly with the history of cloud chambers, followed 
by a demonstration of a simple diffusion-type cloud 
chamber. The latter chamber, as demonstrated, was quite 
similar to that described in recent literature except that 
it was constructed of materials more readily available to a 
small laboratory, and its design overcame a number of 
disadvantages inherent in the earlier models. 


18. Undergraduate determination of alpha-particle 
range. WALTER C. MICHELS and Eva WIENER, Bryn 
Mawr College-—Simple apparatus, using a scintillation 
counter with a 931-A photomultiplier, was described. The 
range of alpha-particles from Po was determined with 
this apparatus by varying the pressure in a shielded con- 
tainer, keeping the geometry constant. Data were pre- 
sented on a 0.6 millicurie source, yielding a range of 
(3.810.03) cm at 760 mm Hg and 15°C. This value was 
obtained after correction for finite source thickness. The 
experiment offers particular value in giving students an 
opportunity to consider statistical fluctuations and to 
analyze the effects of finite source thickness, straggling, etc. 
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19. Quantitative experiment with the forced mechanical 
oscillator. JOHN R. PRUETT and WALTER C. MICHELS, 
Bryn Mawr College—The forced mechanical oscillator 
offers pedagogical advantages as compared with the elec- 
trical oscillating circuit, because it can be introduced at an 
earlier stage of the student’s career and because its charac- 
teristics are directly observable. In the experiment that 
was described, a mass of 9.06 K is suspended from a spring 
with a constant of 670 N/M, giving a natural frequency 
of 1.4/S. The driving force is supplied through a lower 
spring with a constant of about 100 N/M, and guides are 
used to suppress unwanted modes. The driving spring is 
connected to a variable speed drive, which may be adjusted 
to give frequencies between 0.5/S and 2.0/S. Observations 
of amplitude and amplitude and phase angle (measured 
stroboscopically) reveal the characteristics of a coupled 
oscillator. A straightforward analysis allows the force 
applied to the simple oscillator to be obtained. The re- 
sulting amplitude vs frequency curve agrees quite well 
with those derived by Hartog for forced oscillators with 
mixed viscous and Coulombic friction. 


20. Teaching a philosophy of experimentation in a course 
in electrical measurements. FRANCIS T. WORRELL, Rens- 
selaer Polytechnic Institute-—A course in electrical measure- 
ments is an excellent place in which to teach a philosophy 
of experimentation; the basic ideas are easy to grasp, so the 
student can devote more effort to understanding the experi- 
mental process; the equipment is quite well standardized, 
and generally works, and the experiment is not usually 
confused by the intangible effects that are so annoying in, 
say, thermal measurements. 

In the laboratory under discussion, emphasis was 
placed on the experimentation. The student was required, 
therefore, to keep a notebook of the type he would use in a 
research laboratory, and no reports were asked for. Experi- 
ments were chosen not on the basis of specific electrical 
measuring techniques but to illustrate certain general ex- 
perimental procedures. Ample time was allowed for the 
performance of each experiment, so only a few were done 
during the term. 

Certain of the experiments were discussed by way of 
illustration. 


21. An attempt to make the elementary laboratory a 
more effective scientific experience. GLEN F. PIPPERT 
and DuANE ROLLER, Wabash College—A laboratory pro- 
cedure has been tried that is similar in many respects to 
those described by F. K. Richtmeyer and others.! It is 
designed to present experimental problems to the student 
much as he might be faced with them when working in an 
actual laboratory. As much as possible, the preparation for 
the laboratory period is left to the student; this is done in 
part to encourage preliminary library study. The labora- 
tory instructor has a position similar to a director of a 
laboratory. The student hands in a progress report at the 
end of each 3-hour laboratory period; when a particular 
problem is finished, he prepares a report somewhat like one 
he might submit to a periodical. The work to be done is 








467 


outlined by the instructor, but the initiative for selecting 
the particular experimental methods is left to the student, 
commensurate with the equipment available. Acquaintance 
with apparatus, such as measuring instruments, is acquired 
as an incidental part of the main problem not as a finger 
exercise in a laboratory period by itself. The work is 
arranged so that the student uses a number of instruments 
repeatedly throughout the course, thus gaining the 
familiarity with them needed to promote self-confidence. 

1 Committee Report. Am. J. Phys. (Am. Phys. Teacher) 2, 104 (1934); 


G. B. Welch, Am. J. Phys. 3, 69 (1935); F. D. Cruikshank, Am. J. Phys. 
17, 15 (1949). 


22. The effect of multiple-response quizzes on the lin- 
guistic skill of physics students; a conjecture. JuLius 
SUMNER MILLER, Dillard University.—Multiple-response 
quizzes and examinations have taken hold quite firmly in 
our physics teaching. A student can now go through an 
entire course without ever writing a connected account of 
a topic or without ever giving an old-fashioned classroom 
recitation. He is rarely called upon to organize his informa- 
tion and to discourse on it fluently. 

This argument can be extended to students in all 
quarters. Vast numbers now come up from the grades and 
continue even into graduate school without ever having 
written a connected account of anything. Every subject 
has its formalized, standard, multiple-response questions. 
They are used even in English and in Literature! 

These testing devices have efficacy for some purposes, 
but close observation over a dozen years or more points up 
this: the student is deprived of an exercise which has neither 
peer nor superior. The discipline of an essay examination 
cannot be met by any means save writing an essay. The 
student who writes nothing but multiple-choice quizzes 
and examinations develops a natural immunity for exposi- 
sition, both oral and written, and is shamefully lacking in 
value-judgment, thought, and composition. 

I believe it time to at least abate this practice, and there 


may be no real hurt in returning altogether to classical 
methods. 


23. Varieties of objective examinations in general 
physics. Pau KirKPATRICK, Stanford University.—Sample 
examinations with original and perhaps novel features 
were described, and exhibited in part by lantern slides. 
Virtues and limitations of the several types were discussed. 
Proposals for the interchange of examination questions 
among physics teachers were considered. 


24. What happens to the negative ions in the electrolysis 
of copper? Francis E. ToRow, Cornell College.—The usual 
presentation of electrolysis in physics texts leaves the im- 
pression that only those ions which enter or leave the 
solution participate in the charge transfer. But the electric 
field causes all ions to migrate, so that the concentration 
of the electrolyte will generally increasingly depart from 
uniformity—and may greatly reduce the flow of current 
—if convection currents are prevented. These effects are 
easy to demonstrate experimentally. 
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25. Potentials and currents in dc circuits. J. L. BowMAN, 
McPherson College.—Potential is a concept which may be 
used very effectively in teaching beginning students the 
fundamentals of electric currents. The potential, ¢ at a 
point may be defined as the work done in carrying one 
coulomb from the earth to the point, if the potential of the 
earth is taken as zero. Ohm’s law may be stated in two 
forms, v=ri and E=RI. Potential may be calculated using 
v=ri in which 7 is the resistance of any resistor, 7 is the 
current through it, and v is the potential difference. 
E=RI can be applied to a complete circuit, in which E 
is the combined emf, R is the combined resistance, and I 
is the main current. By using these laws in a definite order 


it is very easy to calculate all the currents in complicated 
circuits. 


26. A laboratory method of analyzing wave forms in ac 
networks. RosE A. CARNEY and DONALD SCHEIBER, St. 
Procopius College.—A description was given of a laboratory 
experiment for which the equipment is easily constructed 
by the students. The measurements are carried out in the 
low frequency range so that the student can view the wave 
forms on the oscilloscope. A number of wave forms were 
analyzed such as the half-wave rectified form, the wave 
form obtained from two generators in series operating in a 
circuit in which there is a nonlinear resistance. 


27. Transients in L-C networks. WILLIAM J. SHONKA 
and RoBerT Lari, St. Procopius College.—An experimnet 
was described in which quantitative measurements were 
made from photographic records of the event with the 
60-cycle ac being used as a reference standard. 


28. Quick analysis of musical tones. NEwTon GAINEs, 
Texas Christian University—Brightness of the flashing 
lamp in a Stroboconn (precision frequency meter described 
by Young and Loomis, J. Opt. Soc. Am., October, 1938; 
Kent, Electronics, September, 1943) closely follows the 
wave form of the sound entering the microphone; hence 
approximate analysis of a musical tone may be accom- 
plished (in less than five minutes) by ‘‘stopping’’ the 
fundamental and harmonics at the appropriate windows. 
Thus, if the tone is middle C, the C window will indicate 
the fundamental, second, fourth, and eighth harmonics on 
successive circles of ‘“‘spokes”; the G window, after slight 
adjustment of the tuning fork, will reveal third and sixth 
harmonics; the E window the fifth harmonic; and the B-flat 
window, the seventh. An appropriate portion of a gray 
scale is used to estimate relative intensities on a basis of 
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ten degrees of intensity for the harmonic indicated by the 
darkest, most distinct circle of spokes. 


29. Manpower resources in physics. MArsH W. WHITE, 
The Pennsylvania State College—Data were presented to 
show the number of physicists-in-training graduating at 
various levels during the academic years from 1946-1947 to 
the present. Following World War II there was a vast in- 
crease in the numbers of physics majors. The number of 
bachelor’s degrees awarded by 525 institutions reached a 
maximum of 3770 in 1949-1950 and declined each year 
thereafter, with 3000 estimated for 1951-1952. The numbers 
of master’s degrees awarded by 158 departments was highest 
at 1065 in 1949-1950 and has fallen off to an estimated 950 
for the current year. At the doctorate level the production 
is continuing to increase, although at a reduced rate, with 
475 estimated for 1951-1952. The total enrollment of under- 
graduate physics majors reached a maximum of 12,670 
in 1949-1950 and has now declined to 10,200. The total 
number of physics graduate students has been nearly 
constant for the last two years at 5750. Data concerning 
professional physicists were presented, these statistics 
having been taken from a study of the physicists enrolled 
in the National Scientific Register. These data showed the 
fields of specialization, age, citizenship status, military 
status, education, employment and income of 6600 physi- 
cists and 1300 graduate students. 


30. The current science activities of UNESCO. Ira M. 
FREEMAN, Rutgers University—Some of the major present 
and projected activities of UNESCO in science on an in- 
ternational basis were described. These include the estab- 
lishment of international research centers and the dissemi- 
nation of scientific information through regional Field 
Science Offices, through publications and exhibitions. The 
speaker served as a member of the Department of Natural 
Sciences of UNESCO last year and was engaged in activities 
connected with science education and popularization. 


31. A small cyclotron for educational purposes. THOMAS 
A. BENHAM, Haverford College. (read by Herbert Williams 
of Stitz Laboratories)—A small cyclotron which will 
produce particles of 1 Mev energy was described. It is 
believed that if several schools are interested in obtaining 
such an instrument it should be possible to construct them 
for about $10,000 each. Interested schools should contact 
Professor Benham. 

R. R. PALMER 
Beloit College 





